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

package w;
 
class nonlinear{ 
 
int mini=50;
int mouse=100;       
 
 
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][2*mini+1][2*mini+1]; 
int[][][] opy= new int[101][2*mini+1][2*mini+1]; 
 
double h;
double z;
 
 
double[][][] ux= new double[101][2*mini+1][2*mini+1]; 
double[][][] wx= new double[101][2*mini+1][2*mini+1]; 
double[][][] v= new double[101][2*mini+1][2*mouse+1];
int[][][] gotow= new int[101][2*mini+1][2*mouse+1];
int[][][] gotob= new int[101][2*mini+1][2*mouse+1];
 
 
 
int[][][] nearu= new int[101][2*mini+1][2*mini+1]; 
int[][][] nearw= new int[101][2*mini+1][2*mini+1]; 
double[][][] chu= new double[101][2*mini+1][2*mini+1]; 
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=-mini;n1<mini+1;n1++){ 
for (n2=-mini;n2<mini+1;n2++){ 
c1=(c[s]+n1)*h; 
l1=(double)(y[s]+n2)/th[s]; 
u1=u(c1,1-l1);   
ux[s][n1+mini][n2+mini]=u1; 
} 
} 
} 
 
 
 
for (s=1;s<100;s++){ 
for (n1=-mini;n1<mini+1;n1++){ 
for (n2=-mini;n2<mini+1;n2++){ 
c1=(c[s]+n1)*h; 
l1=(double)(y[s]+n2)/th[s+1]; 
u1=u(c1,1-l1);   
wx[s][n1+mini][n2+mini]=u1; 
} 
} 
} 
 
 
 
for (s=1;s<101;s++){ 
for (n1=-mini;n1<mini+1;n1++){ 
for (n2=-mini;n2<mini+1;n2++){
nearu[s][n1+mini][n2+mini]=search(ux[s][n1+mini][n2+mini]); 
} 
} 
} 
 
 
for (s=1;s<100;s++){ 
for (n1=-mini;n1<mini+1;n1++){ 
for (n2=-mini;n2<mini+1;n2++){
nearw[s][n1+mini][n2+mini]=search(wx[s][n1+mini][n2+mini]); 
} 
} 
} 
 
for (s=1;s<100;s++){ 
ws[s]=nearw[s][mini][mini];
}
 
int x1,x2;
 
for (s=1;s<100;s++){ 
for (m1=-mini;m1<mini+1;m1++){ 
for (m2=-mini;m2<mini+1;m2++){
maxu=-999;
x1=0;
x2=0;
for (n2=-mini;n2<mini+1;n2++){
n1=n2-m2;
z=0;
if(n1>mini)z=100;
if(n1<-mini)z=100;
if(z>50)n1=0;
u1=ux[s][n1+mini][n2+mini];
mx=nearw[s][n1+mini][n2+mini];
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+mini][m2+mini]=maxu;
opc[s][m1+mini][m2+mini]=x1+c[s];
opy[s][m1+mini][m2+mini]=x2+y[s];
}
}
}
 
 
 
for (m1=-mini;m1<mini+1;m1++){ 
for (m2=-mouse;m2<mouse+1;m2++){
v[1][m1+mini][m2+mouse]=-999;
}
}
 
 
 
for (m1=-mini;m1<mini+1;m1++){ 
for (m2=-mini;m2<mini+1;m2++){
v[1][m1+mini][m2+mouse]=chu[1][m1+mini][m2+mini];
}
}
 
int bx,bxs,p1,bps,ps;
double v1,maxv;
int yy,yys;
 
for (s=2;s<100;s++){ 
for (m1=-mini;m1<mini;m1++){ 
for (m2=-mouse;m2<mouse+1;m2++){
maxv=-999;
ps=0;
bps=0;
for(bx=-10;bx<11;bx++){
u1=chu[s][m1+mini][bx+mini];
 
 
p1=search(u1)-ws[s-1];
bxs=m2-bx;
z=0;
if(bxs>100)z=100;
if(bxs<-100)z=100;
if(p1>10)z=100;
if(p1<-10)z=100;
if(z>50)bxs=0;
if(z>50)p1=0;
v1=u1+v[s-1][p1+mini][bxs+mouse];
if(z>50)v1=-999;
if(v1>maxv)bps=bxs;
if(v1>maxv)ps=p1;
if(v1>maxv)maxv=v1;
}
v[s][m1+mini][m2+mouse]=maxv;
gotob[s][m1+mini][m2+mouse]=bps;
gotow[s][m1+mini][m2+mouse]=ps;
}
}
}
 
 
int nx1,nx2,px;
 
s=100;
maxv=-999;
nx1=0;
nx2=0;
bxs=0;
px=0;
for (n1=-mini;n1<mini+1;n1++){ 
for (n2=-mini;n2<mini+1;n2++){
u1=ux[s][n1+mini][n2+mini];
p1=nearu[s][n1+mini][n2+mini]-ws[s-1];
 
bx=n1-n2;
z=0;
if(p1<-mini)z=100;
if(p1>mini)z=100;
if(z>50)p1=0;
 
v1=u1+v[s-1][p1+mini][bx+mouse];
 
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+mini][bx+mouse];
opw[s-1]=gotow[s][px+mini][bx+mouse];
}
 
s=1;
bx=opb[s];
px=opw[s];
c[s]=opc[s][px+mini][bx+mini];
y[s]=opy[s][px+mini][bx+mini];
 
for(s=2;s<100;s++){
bx=opb[s]-opb[s-1];
px=opw[s];
c[s]=opc[s][px+mini][bx+mini];
y[s]=opy[s][px+mini][bx+mini];
}
 
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; 
} 
 
}