都市い - 経済統計練習帳

package line8;
 
 
public class pro {
 
    int[][] dis=new int[101][101]; 
   int[] x=new int[101];
     int[] y=new int[101];
    int s,sx,number,tx;
    double[] shop=new double[101]; 
    double[] dx=new double[101]; 
     double[] sale=new double[101]; 
     double[] rate=new double[101]; 
 
    public static void main(String[] args) {
        pro trs=new pro();
    }
 
    pro(){
 
        tx=0;
        for(s=1;s<11;s++){
        for(sx=1;sx<11;sx++){ 
        tx=tx+1;
        x[tx]=s;
        y[tx]=sx;
        }}
 
        number=tx;
 
        System.out.println(number);
 
 
         for(s=1;s<101;s++){
        for(sx=1;sx<101;sx++){ 
        int d1=(x[s]-x[sx])*(x[s]-x[sx]);
        d1=d1+(y[s]-y[sx])*(y[s]-y[sx]);
        dis[s][sx]=(int)Math.sqrt(d1);
        }}
 
           for(s=1;s<101;s++){
           shop[s]=100;
           }
 
           int tr;
 
           for(tr=1;tr<1001;tr++){
 
       for(s=1;s<101;s++){
     double d3=0;
      for(sx=1;sx<101;sx++){
 
    double d2=shop[sx]/(1+dis[s][sx]);
     d3=d2+d3;      
         }
    dx[s]=d3;
       }
 
     for(sx=1;sx<101;sx++){
     double d4=0;
     for(s=1;s<101;s++){
    double d5=shop[sx]/(1+dis[s][sx]);
     d4=d4+100*d5/dx[s];      
         }
    sale[sx]=d4;
 
       }
 
            for(sx=1;sx<101;sx++){
            rate[sx]=sale[sx]/shop[sx];
            }       
 
            for(s=1;s<101;s++){
            shop[s]=0.5*shop[s]+0.5*rate[s]*shop[s];
 
            }
 
        System.out.println(tr+","+sale[45]);
 
           }
 
 
 
 
    }
 
 
 
 
 
}