000 - プログラム倉庫

class pro{
public static void main(String args[]){
int s,n1,n2,nx1,nx2;
double th[]=new double[101];
double x[]=new double[101];
double t[]=new double[101];
double rev[][][]=new double[101][3][3];
double u[][][]=new double[101][3][3];
double w[][][]=new double[101][3][3];
double v[][][]=new double[101][3][3];
double beta;
double x1;
double u1,r1,r2,v1,w1;
double p,c,cost;
double px,cx;
double maxr,maxv;
double h;
for (s=1;s<101;s++){
th[s]=s;
}
cost=10;
beta=1;
px=0;
cx=0;
maxr=0;
x1=0;
for (n1=1;n1<61;n1++){
for (n2=1;n2<500;n2++){
p=n1;
c=n2;
r2=0;
for (s=1;s<101;s++){
x1=(th[s]-p)/(2*beta);
if (x1<0)x1=0;
u1=th[s]*x1-beta*x1*x1-p*x1-c;
if (u1<0)x1=0;
r1=p*x1+c-cost*x1;
if (x1==0)r1=0;
r2=r2+r1;
}
if (r2>maxr)px=p;
if (r2>maxr)cx=c;
if (r2>maxr)maxr=r2;
}
}
System.out.println(px);
System.out.println(cx);
System.out.println(maxr);

p=px;
c=cx;
for (s=1;s<101;s++){
x1=(th[s]-p)/(2*beta);
if (x1<0)x1=0;
u1=th[s]*x1-beta*x1*x1-p*x1-c;
if (u1<0)x1=0;
r1=p*x1+c;
if (x1==0)r1=0;
x[s]=x1;
t[s]=r1;
}

h=0.01;
for (s=1;s<101;s++){
for (n1=-1;n1<2;n1++){
for (n2=-1;n2<2;n2++){
x1=x[s]+n1*h;
r1=t[s]+n2*h;
rev[s][n1+1][n2+1]=r1-cost*x1;
}
}
}


for (s=1;s<101;s++){
for (n1=-1;n1<2;n1++){
for (n2=-1;n2<2;n2++){
x1=x[s]+n1*h;
r1=t[s]+n2*h;
u1=th[s]*x1-beta*x1-r1;
if (u1<0)u1=-999;
if (x1<0)u1=-999;
u[s][n1+1][n2+1]=u1;
}
}
}

for (s=1;s<101;s++){
for (n1=-1;n1<2;n1++){
for (n2=-1;n2<2;n2++){
x1=x[s]+n1*h;
r1=t[s]+n2*h;
u[s][n1+1][n2+1]=th[s]*x1-beta*x1-r1;
}
}
}

for (s=1;s<100;s++){
for (n1=-1;n1<2;n1++){
for (n2=-1;n2<2;n2++){
x1=x[s]+n1*h;
r1=t[s]+n2*h;
w[s][n1+1][n2+1]=th[s+1]*x1-beta*x1-r1;
}
}
}

for (n1=-1;n1<2;n1++){
for (n2=-1;n2<2;n2++){
v[1][n1+1][n2+1]=rev[1][n1+1][n2+1];
}
}



for (s=1;s<101;s++){
for (n1=-1;n1<2;n1++){
for (n2=-1;n2<2;n2++){
r1=rev[s][n1+1][n2+1];
u1=u[s][n1+1][n2+1];
maxv=-999;
for (nx1=-1;nx1<2;nx1++){
for (nx2=-1;nx2<2;nx2++){
v1=r1+v[s-1][nx1+1][nx2+1];
w1=w[s-1][nx1+1][nx2+1];
if (w1>u1)v1=-999;
if (v1>maxv)maxv=v1;
}
}
v[s][n1+1][n2+1]=maxv;
}
}
}


}
}