java实现Pareto Optimal算法

//根据对Pareto Optimal算法的理解，用java实现的源代码，贴出来跟大家分享一下
//源代码经测试都是已经可以运行的

public class PO {


public void PO_test(double[][] B){

    double A[][]  = new double[B.length][B[0].length]; //初始的时候将输出矩阵A设为空
    for(int i=0;i<B.length;i++){
    for(int j=0;j<B[i].length;j++)
    A[i][j]=0;
    }
    double cc[][] = new double[B.length][B[0].length]; //下面for循环中利用cc存放每次的满足PO的每一行
    for(int i=0;i<B.length;i++){
    for(int j=0;j<B[i].length;j++)
    cc[i][j]=0;
    }
int sz1 = B.length;//sz1=size(B,1); 取矩阵B的第一列的个数 为8


int jj=0; //jj为b Vector的下标，初始为0
int kk[] = new int[B.length];//kk(8)={0,0,0,0,0,0,0,0}
for(int i=0;i<sz1;i++){
kk[i]=0;
}

double c[][] = new double[B.length][B[0].length]; 
for(int i=0;i<sz1;i++){
    for(int j=0;j<B[i].length;j++)
    c[i][j]=0;
    }

//bb=c;//把矩阵c赋给临时变量bb
double bb[][] = new double[B.length][B[0].length];
for(int i=0;i<sz1;i++){
    for(int j=0;j<B[i].length;j++)
    bb[i][j]=c[i][j];
    }

for(int k=0;k<sz1;k++){ //for k=1:sz1 //k从1到8行

int j=0;
//将第k行的元素放进临时变量ak中
double ak[] = new double[B[0].length];
for(int m=0;m<B[k].length;m++){
ak[m]=B[k][m];
}
//这个for循环主要对1到8行中， 用第k行-第i行的值存入临时的bb的第j行中，j自增1
  for(int i=0;i<sz1;i++){
  if( i!=k){
      j=j+1;
      for(int n=0;n<B[i].length;n++){//  bb(j,:)=ak-B(i,:);
      bb[j-1][n]=ak[n]-B[i][n];
      }
  }
}
  //将bb转置
  double bbt[][] = new double[B[0].length][B.length];
for(int p=0;p<bb.length;p++)
for(int q=0;q<bb[0].length;q++){
    bbt[q][p]=bb[p][q];
    }
//testing
//System.out.println("bb.length="+bb.length);
int s=0;
for(int p=0;p<bb.length-1;p++){
for(int q=0;q<bb[0].length;q++){
if(bbt[q][p]<0){
s++;
break;
}
}
}
if(s>6){
jj=jj+1;
//System.out.println("jj="+jj);
kk[jj-1]=k;
for(int m=0;m<ak.length;m++){
cc[jj-1][m]=ak[m];
}
}
/**   //此时的j=8，对于当前的bb矩阵的1-8行进行转置，如果任意一一列存在小于0的元素，那么将这行元素放入c中，并且记录行号放入kk中
  
**/
}
if(jj>0){
//如果jj不等于0，那么说明c矩阵和kk数组中都有元素存在
System.out.println("OutPut Matrix A:");
for(int p=0;p<jj;p++)
{
for(int q=0;q<cc[0].length;q++)
{
A[p][q]=cc[p][q];

System.out.print(A[p][q]+"  ");
}
System.out.println();
}

}else{
System.out.println("There are no Pareto points. The result is an empty matrix.");
}
System.out.println("The total number of PO points is:"+jj);
System.out.println("Pareto Optimal points：");
for(int i=0;i<jj;i++){
System.out.print(kk[i]+"  ");
}

}


}



//测试用例

public class Run_main {

/**
* @param args
*/
public static void main(String[] args) {
PO po = new PO();
double B[][]={{0 ,-1 ,-2},{1 ,-2 ,-3},{-3, 2, 1},{4 ,0, -2},{2,-2 ,1},{1 ,-1, 2},{2, -1, 1},{0, -2 ,2}};
po.PO_test(B);
}

}