ｓぺえｄ ２ - プログラム倉庫

package w;
 
class nonlinear{ 
 
 
int[] th= new int[101]; 
int[] c= new int[101]; 
int[] y= new int[101]; 
 
int[] opw= new int[101]; 
int[] opb= new int[101]; 
int[][][] opc= new int[101][11][11]; 
int[][][] opy= new int[101][11][11]; 
 
double h;
double z;
 
 
double[][][] ux= new double[101][11][11]; 
double[][][] wx= new double[101][11][11]; 
double[][][] v= new double[101][11][101];
int[][][] gotow= new int[101][11][101];
int[][][] gotob= new int[101][11][101];
 
int[][][] suby= new int[101][11][101];
 
 
int[][][] nearu= new int[101][11][11]; 
int[][][] nearw= new int[101][11][11]; 
double[][][] chu= new double[101][11][11]; 
int[] ws= new int[101]; 
 
int s,n1,n2,m1,m2,mx;
double c1,l1,u1;
 
double maxu;
 
void makedata(){ 
 
 
for (s=1;s<101;s++){ 
for (n1=-5;n1<6;n1++){ 
for (n2=-5;n2<6;n2++){ 
c1=(c[s]+n1)*h; 
l1=(double)(y[s]+n2)/th[s]; 
u1=u(c1,1-l1);   
ux[s][n1+5][n2+5]=u1; 
} 
} 
} 
 
 
 
for (s=1;s<100;s++){ 
for (n1=-5;n1<6;n1++){ 
for (n2=-5;n2<6;n2++){ 
c1=(c[s]+n1)*h; 
l1=(double)(y[s]+n2)/th[s+1]; 
u1=u(c1,1-l1);   
wx[s][n1+5][n2+5]=u1; 
} 
} 
} 
 
 
 
for (s=1;s<101;s++){ 
for (n1=-5;n1<6;n1++){ 
for (n2=-5;n2<6;n2++){
nearu[s][n1+5][n2+5]=search(ux[s][n1+5][n2+5]); 
} 
} 
} 
 
 
for (s=1;s<100;s++){ 
for (n1=-5;n1<6;n1++){ 
for (n2=-5;n2<6;n2++){
nearw[s][n1+5][n2+5]=search(wx[s][n1+5][n2+5]); 
} 
} 
} 
 
for (s=1;s<100;s++){ 
ws[s]=nearw[s][5][5];
}
 
int x1,x2;
 
for (s=1;s<100;s++){ 
for (m1=-5;m1<6;m1++){ 
for (m2=-5;m2<6;m2++){
maxu=-999;
x1=0;
x2=0;
for (n2=-5;n2<6;n2++){
n1=n2-m2;
z=0;
if(n1>5)z=100;
if(n1<-5)z=100;
if(z>50)n1=0;
u1=ux[s][n1+5][n2+5];
mx=nearw[s][n1+5][n2+5];
if(mx>m1+ws[s])z=100;
if(z>50)u1=-999;
if(u1>maxu)x1=n1;
if(u1>maxu)x2=n2;
if(u1>maxu)maxu=u1;
}
chu[s][m1+5][m2+5]=maxu;
opc[s][m1+5][m2+5]=x1+c[s];
opy[s][m1+5][m2+5]=x2+y[s];
}
}
}
 
for (m1=-5;m1<6;m1++){ 
for (m2=-50;m2<51;m2++){
v[1][m1+5][m2+50]=-999;
}
}
 
 
 
for (m1=-5;m1<6;m1++){ 
for (m2=-5;m2<6;m2++){
v[1][m1+5][m2+50]=chu[1][m1+5][m2+5];
suby[1][m1+5][m2+50]=0;
}
}
 
int bx,bxs,p1,bps,ps;
double v1,maxv;
int yy,yys;
 
for (s=2;s<100;s++){ 
for (m1=-5;m1<6;m1++){ 
for (m2=-50;m2<51;m2++){
maxv=-999;
ps=0;
bps=0;
yys=0;
for(bx=-5;bx<5;bx++){
u1=chu[s][m1+5][bx+5];
 
yy=opy[s][m1+5][bx+5];
 
p1=search(u1)-ws[s-1];
bxs=m2-bx;
z=0;
if(bxs>50)z=100;
if(bxs<-50)z=100;
if(p1>5)z=100;
if(p1<-5)z=100;
if(z>50)bxs=0;
if(z>50)p1=0;
 
if(suby[s-1][p1+5][bxs+50]>yy)z=100;
 
v1=u1+v[s-1][p1+5][bxs+50];
if(z>50)v1=-999;
if(v1>maxv)bps=bxs;
if(v1>maxv)ps=p1;
if(v1>maxv)yys=yy;
if(v1>maxv)maxv=v1;
}
v[s][m1+5][m2+50]=maxv;
suby[s][m1+5][m2+50]=yys;
gotob[s][m1+5][m2+50]=bps;
gotow[s][m1+5][m2+50]=ps;
}
}
}
 
 
int nx1,nx2,px;
 
s=100;
maxv=-999;
nx1=0;
nx2=0;
bxs=0;
px=0;
for (n1=-5;n1<6;n1++){ 
for (n2=-5;n2<6;n2++){
u1=ux[s][n1+5][n2+5];
p1=nearu[s][n1+5][n2+5]-ws[s-1];
 
bx=n1-n2;
z=0;
if(p1<-5)z=100;
if(p1>5)z=100;
if(z>50)p1=0;
 
v1=u1+v[s-1][p1+5][bx+50];
 
if(z>50)v1=-999;
if(v1>maxv)nx1=n1;
if(v1>maxv)nx2=n2;
if(v1>maxv)px=p1;
if(v1>maxv)bxs=bx;
 
if(v1>maxv)maxv=v1;
 
 
} 
} 
 
 
opb[99]=bxs;
opw[99]=px;
 
int j;
for(j=1;j<99;j++){
s=100-j;
bx=opb[s];
px=opw[s];
opb[s-1]=gotob[s][px+5][bx+50];
opw[s-1]=gotow[s][px+5][bx+50];
}
 
s=1;
bx=opb[s];
px=opw[s];
c[s]=opc[s][px+5][bx+5];
y[s]=opy[s][px+5][bx+5];
 
for(s=2;s<100;s++){
bx=opb[s]-opb[s-1];
px=opw[s];
c[s]=opc[s][px+5][bx+5];
y[s]=opy[s][px+5][bx+5];
}
 
c[100]=c[100]+nx1;
 
y[100]=y[100]+nx2;
 
System.out.println(maxv);
 
 
 
} 
 
 
 
int search(double u1){ 
double c1;
int n1;
n1=-999;
c1=0;
if(u1>-900)c1=Math.exp(u1);
if(u1>-900)n1=(int)(c1/h);
return n1; 
} 
 
 
double u(double c1,double x1){ 
double u1,c2,x2; 
int pp; 
pp=0; 
c2=c1; 
x2=x1; 
if (c2<0)pp=100; 
if (x2<0)pp=100; 
if (x2>1)pp=100; 
if (x2==0)pp=100; 
if (pp>50)c2=0.5; 
if (pp>50)x2=0.5; 
u1=Math.log(c2)+Math.log(x2); 
if (pp>50)u1=-999; 
return u1; 
} 
 
}