約 789,339 件
https://w.atwiki.jp/toho/pages/6582.html
Ark dark half サークル:鶫 Number Track Name Arranger Lyrics Vocal Original Works Original Tune Length 01 青い星、そのスピード NSY リボン リボン 東方神霊廟 リジッドパラダイス [04 16] 02 Silent シュリンプ Cutie Pai まゆちゃん Cutie Pai まゆちゃん 東方神霊廟 死霊の夜桜 [05 29] 03 ひらひらと nice_s ぽらぽら。 ぽらぽら。 東方神霊廟 夢殿大祀廟 [05 38] 04 E・A・R ワニ リボン リボン 東方神霊廟 小さな欲望の星空 [04 58] 05 St. SHINE nice_s リボン メリーベル 東方神霊廟 大神神話伝 [04 43] 06 千夜の月が綺麗でも、もう願いは叶わない。 ワニ Cutie Pai まゆちゃん Cutie Pai まゆちゃん 東方神霊廟 妖怪寺へようこそ [04 17] 07 8 -eight- NSY リボン リボン 東方神霊廟 ゴーストリード [03 33] 詳細 コミックマーケット83(2012/12/30)にて初頒布 イベント価格:1000円 ショップ価格:1300円(税込:1365円) レビュー 名前 コメント
https://w.atwiki.jp/ce00582/pages/5037.html
class tax{ int s,t,j,n1,n2,q,qx,nx1,nx2,nsc,nsy,e; double th[]= new double[101]; double c[]= new double[101]; double y[]= new double[101]; double ux[][][]= new double[101][3][3]; double wx[][][]= new double[101][3][3]; double v[][][][]= new double[100][3][3][21]; int gotoc[][][][]= new int[100][3][3][21]; int gotoy[][][][]= new int[100][3][3][21]; int gotoq[][][][]= new int[100][3][3][21]; double endv[][]=new double[3][3]; int endc[][]=new int[3][3]; int endy[][]=new int[3][3]; int endq[][]=new int[3][3]; int opc[]= new int[101]; int opy[]= new int[101]; int opq[]= new int[101]; double tl,tr,w1,u1,c1,l1,vs,v1; double h,dc,dy; int pp; public static void main(String args[]){ tax test=new tax(); } tax(){ for (s=1;s 101;s++){ th[s]=0.02*s; } tl=tls(th); tr=trs(tl,th); for (s=1;s 101;s++){ y[s]=th[s]*lx(s,tl,tr,th); c[s]=(1-tl)*y[s]+tr; } h=0.01; t=0; while(t 1000){ for (s=1;s 101;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ c1=c[s]+n1*h; l1=(y[s]+n2*h)/th[s]; u1=u(c1,1-l1); ux[s][n1+1][n2+1]=u1; } } } for (s=1;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ c1=c[s]+n1*h; l1=(y[s]+n2*h)/th[s+1]; u1=u(c1,1-l1); wx[s][n1+1][n2+1]=u1; } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ v[1][n1+1][n2+1][q+10]=-999; } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ q=n2-n1; v[1][n1+1][n2+1][q+10]=ux[1][n1+1][n2+1]; } } seekv(); seekendv(); seekop(); e=0; for (s=1;s 101;s++){ e=e+opc[s]*opc[s]+opy[s]*opy[s]; } for (s=1;s 101;s++){ c[s]=c[s]+opc[s]*h; y[s]=y[s]+opy[s]*h; } System.out.println(e); if (e 2)h=h/2; if (h 0.0001)t=10000; t=t+1; } for (s=1;s 100;s++){ dy=y[s+1]-y[s]; dc=c[s+1]-c[s]; if(dy 0)System.out.println(1-dc/dy); } } static double lx(int s,double tl,double tr,double th[]){ double w1,l1; w1=(1-tl)*th[s]; l1=(w1-tr)/(2*+w1); if (l1 0)l1=0; return l1; } static double bud(double tl,double tr,double th[]){ double c1,y1,b1; int s; b1=0; for (s=1;s 101;s++){ y1=th[s]*lx(s,tl,tr,th); c1=(1-tl)*y1+tr; b1=b1+y1-c1; } return b1; } static double trs(double tl,double th[]){ double tr1,tr2,tr3,b1,b2; int t; tr1=0.1; tr2=0.2; b1=bud(tl,tr1,th); t=0; while(t 100){ b2=bud(tl,tr2,th); tr3=tr2-b2*(tr2-tr1)/(b2-b1); tr1=tr2; tr2=tr3; b1=b2; if (b2*b2 0.0001)t=1000; t=t+1; } return tr2; } static double tls(double th[]){ int m; double tl,tr,w1,maxw,maxtl; maxw=-999; maxtl=0; for (m=10;m 50;m++){ tl=0.01*m; tr=trs(tl,th); w1=wel(tl,tr,th); if (w1 maxw)maxtl=tl; if (w1 maxw)maxw=w1; } return maxtl; } static double wel(double tl,double tr,double th[]){ double c1,l1,w1; int s; w1=0; for (s=1;s 101;s++){ l1=lx(s,tl,tr,th); c1=(1-tl)*th[s]*l1+tr; w1=w1+u(c1,1-l1);; } return w1; } static 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; } void seekv(){ double sy1,sy2; for (s=2;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ u1=ux[s][n1+1][n2+1]; qx=q-n1+n2; pp=0; if (qx 10)pp=100; if (qx -10)pp=100; if (pp 50)qx=0; vs=-999; nsc=0; nsy=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[s-1][nx1+1][nx2+1][qx+10]; if (wx[s-1][nx1+1][nx2+1] u1)v1=-999; sy1=y[s]+n2*h; sy2=y[s-1]+nx2*h; if(sy2 sy1)v1=-999; if (v1 vs)nsc=nx1; if (v1 vs)nsy=nx2; if (v1 vs)vs=v1; } } if (pp 50)vs=-999; gotoc[s][n1+1][n2+1][q+10]=nsc; gotoy[s][n1+1][n2+1][q+10]=nsy; gotoq[s][n1+1][n2+1][q+10]=qx; v[s][n1+1][n2+1][q+10]=vs; } } } } } void seekendv(){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ u1=ux[100][n1+1][n2+1]; qx=n2-n1; vs=-999; nsc=0; nsy=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[99][nx1+1][nx2+1][qx+10]; if (wx[99][nx1+1][nx2+1] u1)v1=-999; if (v1 vs)nsc=nx1; if (v1 vs)nsy=nx2; if (v1 vs)vs=v1; } } endc[n1+1][n2+1]=nsc; endy[n1+1][n2+1]=nsy; endq[n1+1][n2+1]=qx; endv[n1+1][n2+1]=vs; } } } void seekop(){ vs=-999; nsc=0; nsy=0; for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ if (endv[n1+1][n2+1] vs)nsc=n1; if (endv[n1+1][n2+1] vs)nsy=n2; if (endv[n1+1][n2+1] vs)vs=endv[n1+1][n2+1]; } } opc[100]=nsc; opy[100]=nsy; opc[99]=endc[opc[100]+1][opy[100]+1]; opy[99]=endy[opc[100]+1][opy[100]+1]; opq[99]=endq[opc[100]+1][opy[100]+1]; for (j=1;j 99;j++){ s=99-j; opc[s]=gotoc[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opy[s]=gotoy[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opq[s]=gotoq[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; } } }
https://w.atwiki.jp/ce00582/pages/2016.html
class tax103{ public static void main(String args[]){ int s,t,j,n1,n2,q,qx,nx1,nx2,nsc,nsy,e; double th[]= new double[101]; double c[]= new double[101]; double y[]= new double[101]; double ux[][][]= new double[101][3][3]; double wx[][][]= new double[101][3][3]; double v[][][][]= new double[100][3][3][21]; int gotoc[][][][]= new int[100][3][3][21]; int gotoy[][][][]= new int[100][3][3][21]; int gotoq[][][][]= new int[100][3][3][21]; double endv[][]=new double[3][3]; int endc[][]=new int[3][3]; int endy[][]=new int[3][3]; int endq[][]=new int[3][3]; int opc[]= new int[101]; int opy[]= new int[101]; int opq[]= new int[101]; double tl,tr,w1,u1,c1,l1,vs,v1; double h,dc,dy; int pp; for (s=1;s 101;s++){ th[s]=0.02*s; } tl=tls(th); tr=trs(tl,th); for (s=1;s 101;s++){ y[s]=th[s]*lx(s,tl,tr,th); c[s]=(1-tl)*y[s]+tr; } h=0.01; t=0; while(t 1000){ for (s=1;s 101;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ c1=c[s]+n1*h; l1=(y[s]+n2*h)/th[s]; u1=u(c1,1-l1); ux[s][n1+1][n2+1]=u1; } } } for (s=1;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ c1=c[s]+n1*h; l1=(y[s]+n2*h)/th[s+1]; u1=u(c1,1-l1); wx[s][n1+1][n2+1]=u1; } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ v[1][n1+1][n2+1][q+10]=-999; } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ q=n2-n1; v[1][n1+1][n2+1][q+10]=ux[1][n1+1][n2+1]; } } for (s=2;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ u1=ux[s][n1+1][n2+1]; qx=q-n1+n2; pp=0; if (qx 10)pp=100; if (qx -10)pp=100; if (pp 50)qx=0; vs=-999; nsc=0; nsy=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[s-1][nx1+1][nx2+1][qx+10]; if (wx[s-1][nx1+1][nx2+1] u1)v1=-999; if (v1 vs)nsc=nx1; if (v1 vs)nsy=nx2; if (v1 vs)vs=v1; } } if (pp 50)vs=-999; gotoc[s][n1+1][n2+1][q+10]=nsc; gotoy[s][n1+1][n2+1][q+10]=nsy; gotoq[s][n1+1][n2+1][q+10]=qx; v[s][n1+1][n2+1][q+10]=vs; } } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ u1=ux[100][n1+1][n2+1]; qx=n2-n1; vs=-999; nsc=0; nsy=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[99][nx1+1][nx2+1][qx+10]; if (wx[99][nx1+1][nx2+1] u1)v1=-999; if (v1 vs)nsc=nx1; if (v1 vs)nsy=nx2; if (v1 vs)vs=v1; } } endc[n1+1][n2+1]=nsc; endy[n1+1][n2+1]=nsy; endq[n1+1][n2+1]=qx; endv[n1+1][n2+1]=vs; } } vs=-999; nsc=0; nsy=0; for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ if (endv[n1+1][n2+1] vs)nsc=n1; if (endv[n1+1][n2+1] vs)nsy=n2; if (endv[n1+1][n2+1] vs)vs=endv[n1+1][n2+1]; } } opc[100]=nsc; opy[100]=nsy; opc[99]=endc[opc[100]+1][opy[100]+1]; opy[99]=endy[opc[100]+1][opy[100]+1]; opq[99]=endq[opc[100]+1][opy[100]+1]; for (j=1;j 99;j++){ s=99-j; opc[s]=gotoc[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opy[s]=gotoy[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opq[s]=gotoq[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; } e=0; for (s=1;s 101;s++){ e=e+opc[s]*opc[s]+opy[s]*opy[s]; } for (s=1;s 101;s++){ c[s]=c[s]+opc[s]*h; y[s]=y[s]+opy[s]*h; } System.out.println(e); if (e 2)h=h/2; if (h 0.0001)t=10000; t=t+1; } for (s=1;s 100;s++){ dy=y[s+1]-y[s]; dc=c[s+1]-c[s]; if (dy 0)System.out.println(1-dc/dy); } } static double lx(int s,double tl,double tr,double th[]){ double w1,l1; w1=(1-tl)*th[s]; l1=(w1-tr)/(2*+w1); if (l1 0)l1=0; return l1; } static double bud(double tl,double tr,double th[]){ double c1,y1,b1; int s; b1=0; for (s=1;s 101;s++){ y1=th[s]*lx(s,tl,tr,th); c1=(1-tl)*y1+tr; b1=b1+y1-c1; } return b1; } static double trs(double tl,double th[]){ double tr1,tr2,tr3,b1,b2; int t; tr1=0.1; tr2=0.2; b1=bud(tl,tr1,th); t=0; while(t 100){ b2=bud(tl,tr2,th); tr3=tr2-b2*(tr2-tr1)/(b2-b1); tr1=tr2; tr2=tr3; b1=b2; if (b2*b2 0.0001)t=1000; t=t+1; } return tr2; } static double tls(double th[]){ int m; double tl,tr,w1,maxw,maxtl; maxw=-999; maxtl=0; for (m=10;m 50;m++){ tl=0.01*m; tr=trs(tl,th); w1=wel(tl,tr,th); if (w1 maxw)maxtl=tl; if (w1 maxw)maxw=w1; } return maxtl; } static double wel(double tl,double tr,double th[]){ double c1,l1,w1; int s; w1=0; for (s=1;s 101;s++){ l1=lx(s,tl,tr,th); c1=(1-tl)*th[s]*l1+tr; w1=w1+u(c1,1-l1);; } return w1; } static 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; } }
https://w.atwiki.jp/ce00582/pages/2210.html
class tax{ public static void main(String args[]){ int s,t,j,n1,n2,q,qx,nx1,nx2,nsc,nsy,e; double th[]= new double[101]; double c[]= new double[101]; double y[]= new double[101]; double ux[][][]= new double[101][3][3]; double wx[][][]= new double[101][3][3]; double v[][][][]= new double[100][3][3][21]; int gotoc[][][][]= new int[100][3][3][21]; int gotoy[][][][]= new int[100][3][3][21]; int gotoq[][][][]= new int[100][3][3][21]; double endv[][]=new double[3][3]; int endc[][]=new int[3][3]; int endy[][]=new int[3][3]; int endq[][]=new int[3][3]; int opc[]= new int[101]; int opy[]= new int[101]; int opq[]= new int[101]; double tl,tr,w1,u1,c1,l1,vs,v1; double h,dc,dy; int pp; for (s=1;s 101;s++){ th[s]=0.02*s; } tl=tls(th); tr=trs(tl,th); for (s=1;s 101;s++){ y[s]=th[s]*lx(s,tl,tr,th); c[s]=(1-tl)*y[s]+tr; } h=0.01; t=0; while(t 1000){ for (s=1;s 101;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ c1=c[s]+n1*h; l1=(y[s]+n2*h)/th[s]; u1=u(c1,1-l1); ux[s][n1+1][n2+1]=u1; } } } for (s=1;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ c1=c[s]+n1*h; l1=(y[s]+n2*h)/th[s+1]; u1=u(c1,1-l1); wx[s][n1+1][n2+1]=u1; } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ v[1][n1+1][n2+1][q+10]=-999; } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ q=n2-n1; v[1][n1+1][n2+1][q+10]=ux[1][n1+1][n2+1]; } } for (s=2;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ u1=ux[s][n1+1][n2+1]; qx=q-n1+n2; pp=0; if (qx 10)pp=100; if (qx -10)pp=100; if (pp 50)qx=0; vs=-999; nsc=0; nsy=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[s-1][nx1+1][nx2+1][qx+10]; if (wx[s-1][nx1+1][nx2+1] u1)v1=-999; if (v1 vs)nsc=nx1; if (v1 vs)nsy=nx2; if (v1 vs)vs=v1; } } if (pp 50)vs=-999; gotoc[s][n1+1][n2+1][q+10]=nsc; gotoy[s][n1+1][n2+1][q+10]=nsy; gotoq[s][n1+1][n2+1][q+10]=qx; v[s][n1+1][n2+1][q+10]=vs; } } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ u1=ux[100][n1+1][n2+1]; qx=n2-n1; vs=-999; nsc=0; nsy=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[99][nx1+1][nx2+1][qx+10]; if (wx[99][nx1+1][nx2+1] u1)v1=-999; if (v1 vs)nsc=nx1; if (v1 vs)nsy=nx2; if (v1 vs)vs=v1; } } endc[n1+1][n2+1]=nsc; endy[n1+1][n2+1]=nsy; endq[n1+1][n2+1]=qx; endv[n1+1][n2+1]=vs; } } vs=-999; nsc=0; nsy=0; for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ if (endv[n1+1][n2+1] vs)nsc=n1; if (endv[n1+1][n2+1] vs)nsy=n2; if (endv[n1+1][n2+1] vs)vs=endv[n1+1][n2+1]; } } opc[100]=nsc; opy[100]=nsy; opc[99]=endc[opc[100]+1][opy[100]+1]; opy[99]=endy[opc[100]+1][opy[100]+1]; opq[99]=endq[opc[100]+1][opy[100]+1]; for (j=1;j 99;j++){ s=99-j; opc[s]=gotoc[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opy[s]=gotoy[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opq[s]=gotoq[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; } e=0; for (s=1;s 101;s++){ e=e+opc[s]*opc[s]+opy[s]*opy[s]; } for (s=1;s 101;s++){ c[s]=c[s]+opc[s]*h; y[s]=y[s]+opy[s]*h; } System.out.println(e); if (e 2)h=h/2; if (h 0.0001)t=10000; t=t+1; } for (s=1;s 100;s++){ dy=y[s+1]-y[s]; dc=c[s+1]-c[s]; if (dy 0)System.out.println(1-dc/dy); } } static double lx(int s,double tl,double tr,double th[]){ double w1,l1; w1=(1-tl)*th[s]; l1=(w1-tr)/(2*+w1); if (l1 0)l1=0; return l1; } static double bud(double tl,double tr,double th[]){ double c1,y1,b1; int s; b1=0; for (s=1;s 101;s++){ y1=th[s]*lx(s,tl,tr,th); c1=(1-tl)*y1+tr; b1=b1+y1-c1; } return b1; } static double trs(double tl,double th[]){ double tr1,tr2,tr3,b1,b2; int t; tr1=0.1; tr2=0.2; b1=bud(tl,tr1,th); t=0; while(t 100){ b2=bud(tl,tr2,th); tr3=tr2-b2*(tr2-tr1)/(b2-b1); tr1=tr2; tr2=tr3; b1=b2; if (b2*b2 0.0001)t=1000; t=t+1; } return tr2; } static double tls(double th[]){ int m; double tl,tr,w1,maxw,maxtl; maxw=-999; maxtl=0; for (m=10;m 50;m++){ tl=0.01*m; tr=trs(tl,th); w1=wel(tl,tr,th); if (w1 maxw)maxtl=tl; if (w1 maxw)maxw=w1; } return maxtl; } static double wel(double tl,double tr,double th[]){ double c1,l1,w1; int s; w1=0; for (s=1;s 101;s++){ l1=lx(s,tl,tr,th); c1=(1-tl)*th[s]*l1+tr; w1=w1+u(c1,1-l1);; } return w1; } static 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; } }
https://w.atwiki.jp/ce00582/pages/2211.html
class tax{ public static void main(String args[]){ int s,t,j,n1,n2,q,qx,nx1,nx2,nsc,nsy,e; double th[]= new double[101]; double c[]= new double[101]; double y[]= new double[101]; double ux[][][]= new double[101][3][3]; double wx[][][]= new double[101][3][3]; double v[][][][]= new double[100][3][3][21]; int gotoc[][][][]= new int[100][3][3][21]; int gotoy[][][][]= new int[100][3][3][21]; int gotoq[][][][]= new int[100][3][3][21]; double endv[][]=new double[3][3]; int endc[][]=new int[3][3]; int endy[][]=new int[3][3]; int endq[][]=new int[3][3]; int opc[]= new int[101]; int opy[]= new int[101]; int opq[]= new int[101]; double tl,tr,w1,u1,c1,l1,vs,v1; double h,dc,dy; int pp; for (s=1;s 101;s++){ th[s]=seekth(s); } tl=tls(th); tr=trs(tl,th); for (s=1;s 101;s++){ y[s]=th[s]*lx(s,tl,tr,th); c[s]=(1-tl)*y[s]+tr; } h=0.01; t=0; while(t 1000){ for (s=1;s 101;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ c1=c[s]+n1*h; l1=(y[s]+n2*h)/th[s]; u1=u(c1,1-l1); ux[s][n1+1][n2+1]=u1; } } } for (s=1;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ c1=c[s]+n1*h; l1=(y[s]+n2*h)/th[s+1]; u1=u(c1,1-l1); wx[s][n1+1][n2+1]=u1; } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ v[1][n1+1][n2+1][q+10]=-999; } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ q=n2-n1; v[1][n1+1][n2+1][q+10]=ux[1][n1+1][n2+1]; } } for (s=2;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ u1=ux[s][n1+1][n2+1]; qx=q-n1+n2; pp=0; if (qx 10)pp=100; if (qx -10)pp=100; if (pp 50)qx=0; vs=-999; nsc=0; nsy=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[s-1][nx1+1][nx2+1][qx+10]; if (wx[s-1][nx1+1][nx2+1] u1)v1=-999; if (v1 vs)nsc=nx1; if (v1 vs)nsy=nx2; if (v1 vs)vs=v1; } } if (pp 50)vs=-999; gotoc[s][n1+1][n2+1][q+10]=nsc; gotoy[s][n1+1][n2+1][q+10]=nsy; gotoq[s][n1+1][n2+1][q+10]=qx; v[s][n1+1][n2+1][q+10]=vs; } } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ u1=ux[100][n1+1][n2+1]; qx=n2-n1; vs=-999; nsc=0; nsy=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[99][nx1+1][nx2+1][qx+10]; if (wx[99][nx1+1][nx2+1] u1)v1=-999; if (v1 vs)nsc=nx1; if (v1 vs)nsy=nx2; if (v1 vs)vs=v1; } } endc[n1+1][n2+1]=nsc; endy[n1+1][n2+1]=nsy; endq[n1+1][n2+1]=qx; endv[n1+1][n2+1]=vs; } } vs=-999; nsc=0; nsy=0; for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ if (endv[n1+1][n2+1] vs)nsc=n1; if (endv[n1+1][n2+1] vs)nsy=n2; if (endv[n1+1][n2+1] vs)vs=endv[n1+1][n2+1]; } } opc[100]=nsc; opy[100]=nsy; opc[99]=endc[opc[100]+1][opy[100]+1]; opy[99]=endy[opc[100]+1][opy[100]+1]; opq[99]=endq[opc[100]+1][opy[100]+1]; for (j=1;j 99;j++){ s=99-j; opc[s]=gotoc[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opy[s]=gotoy[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opq[s]=gotoq[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; } e=0; for (s=1;s 101;s++){ e=e+opc[s]*opc[s]+opy[s]*opy[s]; } for (s=1;s 101;s++){ c[s]=c[s]+opc[s]*h; y[s]=y[s]+opy[s]*h; } System.out.println(e); if (e 2)h=h/2; if (h 0.0001)t=10000; t=t+1; } for (s=1;s 100;s++){ dy=y[s+1]-y[s]; dc=c[s+1]-c[s]; if (dy 0)System.out.println(1-dc/dy); } } static double lx(int s,double tl,double tr,double th[]){ double w1,l1; w1=(1-tl)*th[s]; l1=(w1-tr)/(2*+w1); if (l1 0)l1=0; return l1; } static double bud(double tl,double tr,double th[]){ double c1,y1,b1; int s; b1=0; for (s=1;s 101;s++){ y1=th[s]*lx(s,tl,tr,th); c1=(1-tl)*y1+tr; b1=b1+y1-c1; } return b1; } static double trs(double tl,double th[]){ double tr1,tr2,tr3,b1,b2; int t; tr1=0.1; tr2=0.2; b1=bud(tl,tr1,th); t=0; while(t 100){ b2=bud(tl,tr2,th); tr3=tr2-b2*(tr2-tr1)/(b2-b1); tr1=tr2; tr2=tr3; b1=b2; if (b2*b2 0.0001)t=1000; t=t+1; } return tr2; } static double tls(double th[]){ int m; double tl,tr,w1,maxw,maxtl; maxw=-999; maxtl=0; for (m=10;m 50;m++){ tl=0.01*m; tr=trs(tl,th); w1=wel(tl,tr,th); if (w1 maxw)maxtl=tl; if (w1 maxw)maxw=w1; } return maxtl; } static double wel(double tl,double tr,double th[]){ double c1,l1,w1; int s; w1=0; for (s=1;s 101;s++){ l1=lx(s,tl,tr,th); c1=(1-tl)*th[s]*l1+tr; w1=w1+u(c1,1-l1);; } return w1; } static 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; } static double seekth(int s){ double p; double mu; double sig; double yy; double ths; mu=0; sig=0.39; p=0.01*s-0.005; yy=seeky(p,mu,sig); ths=Math.exp(yy); return ths; } static double f(double x,double mu,double sig){ double pi,x1,x2,x3,fx; pi = 3.1415; x1=-Math.pow(x - mu,2) / (2*Math.pow(sig,2)); x2=Math.exp(x1); x3=sig*Math.pow(2*pi,0.5); fx=x2/x3; return fx; } static double g(double y,double mu, double sig){ double gx,h,x; int n,t; gx=0; h=0.001; t=(int)(y/h); for (n=-2000;n t;n++){ x=n*h; gx=gx+f(x,mu,sig)*h; } return gx; } static double seeky(double p,double mu,double sig){ double g1,g2,y1,y2,y3; int t; y1=0.4; y2=-0.2; g1=g(y1,mu,sig); t=0; while(t 100){ g2=g(y2,mu,sig); y3=y2+(p-g2)*(y2-y1)/(g2-g1); y1=y2; y2=y3; g1=g2; if (Math.pow(p-g2,2) 0.0001)t=1000; t=t+1; } return y2; } }
https://w.atwiki.jp/ce00582/pages/4174.html
import java.awt.*; import java.awt.event.*; class tax0313 extends Frame{ int s,t,j,n1,n2,q,qx,nx1,nx2,nsc,nsy,e; double th[]= new double[101]; double c[]= new double[101]; double y[]= new double[101]; double ux[][][]= new double[101][3][3]; double wx[][][]= new double[101][3][3]; double v[][][][]= new double[100][3][3][21]; int gotoc[][][][]= new int[100][3][3][21]; int gotoy[][][][]= new int[100][3][3][21]; int gotoq[][][][]= new int[100][3][3][21]; double endv[][]=new double[3][3]; int endc[][]=new int[3][3]; int endy[][]=new int[3][3]; int endq[][]=new int[3][3]; int opc[]= new int[101]; int opy[]= new int[101]; int opq[]= new int[101]; double tl,tr,w1,u1,c1,l1,vs,v1; double h,dc,dy; int pp; public static void main(String args[]){ tax0313 f=new tax0313(); f.setTitle("tax0313"); f.setSize(700,700); f.setBackground(Color.yellow); f.setVisible(true); } tax0313(){ for (s=1;s 101;s++){ th[s]=0.02*s; } tl=tls(th); tr=trs(tl,th); for (s=1;s 101;s++){ y[s]=th[s]*lx(s,tl,tr,th); c[s]=(1-tl)*y[s]+tr; } h=0.01; t=0; while(t 1000){ step(); System.out.println(e); if (e 2)h=h/2; if (h 0.0001)t=10000; t=t+1; } for (s=1;s 100;s++){ dy=y[s+1]-y[s]; dc=c[s+1]-c[s]; if (dy 0)System.out.println(1-dc/dy); } } static double lx(int s,double tl,double tr,double th[]){ double w1,l1; w1=(1-tl)*th[s]; l1=(w1-tr)/(2*+w1); if (l1 0)l1=0; return l1; } static double bud(double tl,double tr,double th[]){ double c1,y1,b1; int s; b1=0; for (s=1;s 101;s++){ y1=th[s]*lx(s,tl,tr,th); c1=(1-tl)*y1+tr; b1=b1+y1-c1; } return b1; } static double trs(double tl,double th[]){ double tr1,tr2,tr3,b1,b2; int t; tr1=0.1; tr2=0.2; b1=bud(tl,tr1,th); t=0; while(t 100){ b2=bud(tl,tr2,th); tr3=tr2-b2*(tr2-tr1)/(b2-b1); tr1=tr2; tr2=tr3; b1=b2; if (b2*b2 0.0001)t=1000; t=t+1; } return tr2; } static double tls(double th[]){ int m; double tl,tr,w1,maxw,maxtl; maxw=-999; maxtl=0; for (m=10;m 50;m++){ tl=0.01*m; tr=trs(tl,th); w1=wel(tl,tr,th); if (w1 maxw)maxtl=tl; if (w1 maxw)maxw=w1; } return maxtl; } static double wel(double tl,double tr,double th[]){ double c1,l1,w1; int s; w1=0; for (s=1;s 101;s++){ l1=lx(s,tl,tr,th); c1=(1-tl)*th[s]*l1+tr; w1=w1+u(c1,1-l1);; } return w1; } static 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; } void step(){ for (s=1;s 101;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ c1=c[s]+n1*h; l1=(y[s]+n2*h)/th[s]; u1=u(c1,1-l1); ux[s][n1+1][n2+1]=u1; } } } for (s=1;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ c1=c[s]+n1*h; l1=(y[s]+n2*h)/th[s+1]; u1=u(c1,1-l1); wx[s][n1+1][n2+1]=u1; } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ v[1][n1+1][n2+1][q+10]=-999; } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ q=n2-n1; v[1][n1+1][n2+1][q+10]=ux[1][n1+1][n2+1]; } } for (s=2;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ u1=ux[s][n1+1][n2+1]; qx=q-n1+n2; pp=0; if (qx 10)pp=100; if (qx -10)pp=100; if (pp 50)qx=0; vs=-999; nsc=0; nsy=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[s-1][nx1+1][nx2+1][qx+10]; if (wx[s-1][nx1+1][nx2+1] u1)v1=-999; if (v1 vs)nsc=nx1; if (v1 vs)nsy=nx2; if (v1 vs)vs=v1; } } if (pp 50)vs=-999; gotoc[s][n1+1][n2+1][q+10]=nsc; gotoy[s][n1+1][n2+1][q+10]=nsy; gotoq[s][n1+1][n2+1][q+10]=qx; v[s][n1+1][n2+1][q+10]=vs; } } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ u1=ux[100][n1+1][n2+1]; qx=n2-n1; vs=-999; nsc=0; nsy=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[99][nx1+1][nx2+1][qx+10]; if (wx[99][nx1+1][nx2+1] u1)v1=-999; if (v1 vs)nsc=nx1; if (v1 vs)nsy=nx2; if (v1 vs)vs=v1; } } endc[n1+1][n2+1]=nsc; endy[n1+1][n2+1]=nsy; endq[n1+1][n2+1]=qx; endv[n1+1][n2+1]=vs; } } vs=-999; nsc=0; nsy=0; for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ if (endv[n1+1][n2+1] vs)nsc=n1; if (endv[n1+1][n2+1] vs)nsy=n2; if (endv[n1+1][n2+1] vs)vs=endv[n1+1][n2+1]; } } opc[100]=nsc; opy[100]=nsy; opc[99]=endc[opc[100]+1][opy[100]+1]; opy[99]=endy[opc[100]+1][opy[100]+1]; opq[99]=endq[opc[100]+1][opy[100]+1]; for (j=1;j 99;j++){ s=99-j; opc[s]=gotoc[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opy[s]=gotoy[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opq[s]=gotoq[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; } e=0; for (s=1;s 101;s++){ e=e+opc[s]*opc[s]+opy[s]*opy[s]; } for (s=1;s 101;s++){ c[s]=c[s]+opc[s]*h; y[s]=y[s]+opy[s]*h; } } public void paint(Graphics g){ int gy; for (s=1;s 101;s++){ gy=(int)(500*c[s]); g.fillRect(5*s+100,600-gy,5,gy); } }
https://w.atwiki.jp/ce00582/pages/2357.html
class pro0902{ public static void main(String args[]){ int s,t,j,n1,n2,q,qx,nx1,nx2,nsc,nsy,e; double th[]= new double[101]; double c[]= new double[101]; double y[]= new double[101]; double ux[][][]= new double[101][3][3]; double wx[][][]= new double[101][3][3]; double v[][][][]= new double[100][3][3][21]; int gotoc[][][][]= new int[100][3][3][21]; int gotoy[][][][]= new int[100][3][3][21]; int gotoq[][][][]= new int[100][3][3][21]; double endv[][]=new double[3][3]; int endc[][]=new int[3][3]; int endy[][]=new int[3][3]; int endq[][]=new int[3][3]; int opc[]= new int[101]; int opy[]= new int[101]; int opq[]= new int[101]; double tl,tr,w1,u1,c1,l1,vs,v1; double h,dc,dy; int pp; for (s=1;s 101;s++){ th[s]=0.02*s; } tl=tls(th); tr=trs(tl,th); for (s=1;s 101;s++){ y[s]=th[s]*lx(s,tl,tr,th); c[s]=(1-tl)*y[s]+tr; } h=0.01; t=0; while(t 1000){ for (s=1;s 101;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ c1=c[s]+n1*h; l1=(y[s]+n2*h)/th[s]; u1=u(c1,1-l1); ux[s][n1+1][n2+1]=u1; } } } for (s=1;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ c1=c[s]+n1*h; l1=(y[s]+n2*h)/th[s+1]; u1=u(c1,1-l1); wx[s][n1+1][n2+1]=u1; } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ v[1][n1+1][n2+1][q+10]=-999; } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ q=n2-n1; v[1][n1+1][n2+1][q+10]=ux[1][n1+1][n2+1]; } } for (s=2;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ u1=ux[s][n1+1][n2+1]; qx=q-n1+n2; pp=0; if (qx 10)pp=100; if (qx -10)pp=100; if (pp 50)qx=0; vs=-999; nsc=0; nsy=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[s-1][nx1+1][nx2+1][qx+10]; if (wx[s-1][nx1+1][nx2+1] u1)v1=-999; if (v1 vs)nsc=nx1; if (v1 vs)nsy=nx2; if (v1 vs)vs=v1; } } if (pp 50)vs=-999; gotoc[s][n1+1][n2+1][q+10]=nsc; gotoy[s][n1+1][n2+1][q+10]=nsy; gotoq[s][n1+1][n2+1][q+10]=qx; v[s][n1+1][n2+1][q+10]=vs; } } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ u1=ux[100][n1+1][n2+1]; qx=n2-n1; vs=-999; nsc=0; nsy=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[99][nx1+1][nx2+1][qx+10]; if (wx[99][nx1+1][nx2+1] u1)v1=-999; if (v1 vs)nsc=nx1; if (v1 vs)nsy=nx2; if (v1 vs)vs=v1; } } endc[n1+1][n2+1]=nsc; endy[n1+1][n2+1]=nsy; endq[n1+1][n2+1]=qx; endv[n1+1][n2+1]=vs; } } vs=-999; nsc=0; nsy=0; for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ if (endv[n1+1][n2+1] vs)nsc=n1; if (endv[n1+1][n2+1] vs)nsy=n2; if (endv[n1+1][n2+1] vs)vs=endv[n1+1][n2+1]; } } opc[100]=nsc; opy[100]=nsy; opc[99]=endc[opc[100]+1][opy[100]+1]; opy[99]=endy[opc[100]+1][opy[100]+1]; opq[99]=endq[opc[100]+1][opy[100]+1]; for (j=1;j 99;j++){ s=99-j; opc[s]=gotoc[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opy[s]=gotoy[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opq[s]=gotoq[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; } e=0; for (s=1;s 101;s++){ e=e+opc[s]*opc[s]+opy[s]*opy[s]; } for (s=1;s 101;s++){ c[s]=c[s]+opc[s]*h; y[s]=y[s]+opy[s]*h; } System.out.println(e); if (e 2)h=h/2; if (h 0.0001)t=10000; t=t+1; } for (s=1;s 100;s++){ dy=y[s+1]-y[s]; dc=c[s+1]-c[s]; if (dy 0)System.out.println(1-dc/dy); } } static double lx(int s,double tl,double tr,double th[]){ double w1,l1; w1=(1-tl)*th[s]; l1=(w1-tr)/(2*+w1); if (l1 0)l1=0; return l1; } static double bud(double tl,double tr,double th[]){ double c1,y1,b1; int s; b1=0; for (s=1;s 101;s++){ y1=th[s]*lx(s,tl,tr,th); c1=(1-tl)*y1+tr; b1=b1+y1-c1; } return b1; } static double trs(double tl,double th[]){ double tr1,tr2,tr3,b1,b2; int t; tr1=0.1; tr2=0.2; b1=bud(tl,tr1,th); t=0; while(t 100){ b2=bud(tl,tr2,th); tr3=tr2-b2*(tr2-tr1)/(b2-b1); tr1=tr2; tr2=tr3; b1=b2; if (b2*b2 0.0001)t=1000; t=t+1; } return tr2; } static double tls(double th[]){ int m; double tl,tr,w1,maxw,maxtl; maxw=-999; maxtl=0; for (m=10;m 50;m++){ tl=0.01*m; tr=trs(tl,th); w1=wel(tl,tr,th); if (w1 maxw)maxtl=tl; if (w1 maxw)maxw=w1; } return maxtl; } static double wel(double tl,double tr,double th[]){ double c1,l1,w1; int s; w1=0; for (s=1;s 101;s++){ l1=lx(s,tl,tr,th); c1=(1-tl)*th[s]*l1+tr; w1=w1+u(c1,1-l1);; } return w1; } static 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; } }
https://w.atwiki.jp/ce00582/pages/1960.html
class tax{ public static void main(String args[]){ int s,t,j,n1,n2,q,qx,nx1,nx2,nsc,nsy,e; double th[]= new double[101]; double c[]= new double[101]; double y[]= new double[101]; double ux[][][]= new double[101][3][3]; double wx[][][]= new double[101][3][3]; double v[][][][]= new double[100][3][3][21]; int gotoc[][][][]= new int[100][3][3][21]; int gotoy[][][][]= new int[100][3][3][21]; int gotoq[][][][]= new int[100][3][3][21]; double endv[][]=new double[3][3]; int endc[][]=new int[3][3]; int endy[][]=new int[3][3]; int endq[][]=new int[3][3]; int opc[]= new int[101]; int opy[]= new int[101]; int opq[]= new int[101]; double tl,tr,w1,u1,c1,l1,vs,v1; double h,dc,dy; int pp; for (s=1;s 101;s++){ th[s]=0.02*s; } tl=tls(th); tr=trs(tl,th); for (s=1;s 101;s++){ y[s]=th[s]*lx(s,tl,tr,th); c[s]=(1-tl)*y[s]+tr; } h=0.01; t=0; while(t 1000){ for (s=1;s 101;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ c1=c[s]+n1*h; l1=(y[s]+n2*h)/th[s]; u1=u(c1,1-l1); ux[s][n1+1][n2+1]=u1; } } } for (s=1;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ c1=c[s]+n1*h; l1=(y[s]+n2*h)/th[s+1]; u1=u(c1,1-l1); wx[s][n1+1][n2+1]=u1; } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ v[1][n1+1][n2+1][q+10]=-999; } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ q=n2-n1; v[1][n1+1][n2+1][q+10]=ux[1][n1+1][n2+1]; } } for (s=2;s 100;s++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (q=-10;q 11;q++){ u1=ux[s][n1+1][n2+1]; qx=q-n1+n2; pp=0; if (qx 10)pp=100; if (qx -10)pp=100; if (pp 50)qx=0; vs=-999; nsc=0; nsy=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[s-1][nx1+1][nx2+1][qx+10]; if (wx[s-1][nx1+1][nx2+1] u1)v1=-999; if (v1 vs)nsc=nx1; if (v1 vs)nsy=nx2; if (v1 vs)vs=v1; } } if (pp 50)vs=-999; gotoc[s][n1+1][n2+1][q+10]=nsc; gotoy[s][n1+1][n2+1][q+10]=nsy; gotoq[s][n1+1][n2+1][q+10]=qx; v[s][n1+1][n2+1][q+10]=vs; } } } } for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ u1=ux[100][n1+1][n2+1]; qx=n2-n1; vs=-999; nsc=0; nsy=0; for (nx1=-1;nx1 2;nx1++){ for (nx2=-1;nx2 2;nx2++){ v1=u1+v[99][nx1+1][nx2+1][qx+10]; if (wx[99][nx1+1][nx2+1] u1)v1=-999; if (v1 vs)nsc=nx1; if (v1 vs)nsy=nx2; if (v1 vs)vs=v1; } } endc[n1+1][n2+1]=nsc; endy[n1+1][n2+1]=nsy; endq[n1+1][n2+1]=qx; endv[n1+1][n2+1]=vs; } } vs=-999; nsc=0; nsy=0; for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ if (endv[n1+1][n2+1] vs)nsc=n1; if (endv[n1+1][n2+1] vs)nsy=n2; if (endv[n1+1][n2+1] vs)vs=endv[n1+1][n2+1]; } } opc[100]=nsc; opy[100]=nsy; opc[99]=endc[opc[100]+1][opy[100]+1]; opy[99]=endy[opc[100]+1][opy[100]+1]; opq[99]=endq[opc[100]+1][opy[100]+1]; for (j=1;j 99;j++){ s=99-j; opc[s]=gotoc[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opy[s]=gotoy[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; opq[s]=gotoq[s+1][opc[s+1]+1][opy[s+1]+1][opq[s+1]+10]; } e=0; for (s=1;s 101;s++){ e=e+opc[s]*opc[s]+opy[s]*opy[s]; } for (s=1;s 101;s++){ c[s]=c[s]+opc[s]*h; y[s]=y[s]+opy[s]*h; } System.out.println(e); if (e 2)h=h/2; if (h 0.0001)t=10000; t=t+1; } for (s=1;s 100;s++){ dy=y[s+1]-y[s]; dc=c[s+1]-c[s]; if (dy 0)System.out.println(1-dc/dy); } } static double lx(int s,double tl,double tr,double th[]){ double w1,l1; w1=(1-tl)*th[s]; l1=(w1-tr)/(2*+w1); if (l1 0)l1=0; return l1; } static double bud(double tl,double tr,double th[]){ double c1,y1,b1; int s; b1=0; for (s=1;s 101;s++){ y1=th[s]*lx(s,tl,tr,th); c1=(1-tl)*y1+tr; b1=b1+y1-c1; } return b1; } static double trs(double tl,double th[]){ double tr1,tr2,tr3,b1,b2; int t; tr1=0.1; tr2=0.2; b1=bud(tl,tr1,th); t=0; while(t 100){ b2=bud(tl,tr2,th); tr3=tr2-b2*(tr2-tr1)/(b2-b1); tr1=tr2; tr2=tr3; b1=b2; if (b2*b2 0.0001)t=1000; t=t+1; } return tr2; } static double tls(double th[]){ int m; double tl,tr,w1,maxw,maxtl; maxw=-999; maxtl=0; for (m=10;m 50;m++){ tl=0.01*m; tr=trs(tl,th); w1=wel(tl,tr,th); if (w1 maxw)maxtl=tl; if (w1 maxw)maxw=w1; } return maxtl; } static double wel(double tl,double tr,double th[]){ double c1,l1,w1; int s; w1=0; for (s=1;s 101;s++){ l1=lx(s,tl,tr,th); c1=(1-tl)*th[s]*l1+tr; w1=w1+u(c1,1-l1);; } return w1; } static 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; } }
https://w.atwiki.jp/hmiku/pages/44419.html
【検索用 sensyuu 登録タグ 作S 作SE-H 作り手】 + 目次 目次 特徴 リンク 曲 CD 動画 コメント 【ニコニコ動画】VOCALOID%E7%B3%BB%E8%87%AA%E4%BD%9C%E5%8B%95%E7%94%BB 特徴 作り手名:『sensyuu』(せんしゅう) 和風の曲をよく作るボカロP。 2009年10月、「日暮し」にて初めてVOCALOIDオリジナル曲を作曲。 ピアプロコラボ「自信ないコラボ」に所属。 2021年6月、「夢幻桜」にてついに初めての殿堂入りを達成。 使用VOCALOIDは鏡音リン、鏡音レン、KAITO、VY1、VY2。 リンク piapro 曲 とこしえの華 日暮し 夢幻桜 CD まだCDが登録されていません 動画 コメント 名前 コメント
https://w.atwiki.jp/takaradzuka/pages/143.html
来島わたるをお気に入りに追加 来島わたるとは 来島わたるの41%は鉄の意志で出来ています。来島わたるの33%は優雅さで出来ています。来島わたるの12%は気合で出来ています。来島わたるの8%は柳の樹皮で出来ています。来島わたるの2%は純金で出来ています。来島わたるの2%は見栄で出来ています。来島わたるの1%は犠牲で出来ています。来島わたるの1%は大人の都合で出来ています。 来島わたるの報道 熱海・初島で「ピングー」コラボイベント 「ピングー号」運航やコラボメニュー提供(みんなの経済新聞ネットワーク) - Yahoo!ニュース - Yahoo!ニュース 沖縄八重山諸島オリジナルSDGs学習ゲーム完成! ゲームを起点に地域・企業・自治体が団結し、八重山ならではのローカルSDGs推進へ - PR TIMES 商船三井、LNG燃料自動車船4隻発注。NSYと新来島、24年納期 - 日本海事新聞 【Ship of the Year 2020特集】新来島豊橋造船建造船「SAKURA LEADER」受賞、国内初のLNG燃料大型自動車船 - 日本海事新聞 日本郵船、自動車輸送GHG削減加速。LNG焚き12隻、新来島・NSYと覚書 - 日本海事新聞 愛媛沖貨物船沈没事故 潜水での船内の捜索を終了 今治海保 - 毎日新聞 - 毎日新聞 来島わたるのウィキペディア 来島わたる 来島わたるの掲示板 名前(HN) カキコミ すべてのコメントを見る 来島わたるのリンク #blogsearch2 ページ先頭へ 来島わたる 宝塚歌劇団 このページについて このページは来島わたるのインターネット上の情報を時系列に網羅したリンク集のようなものです。ブックマークしておけば、日々更新される来島わたるに関連する最新情報にアクセスすることができます。 情報収集はプログラムで行っているため、名前が同じであるが異なるカテゴリーの情報が掲載される場合があります。ご了承ください。 リンク先の内容を保証するものではありません。ご自身の責任でクリックしてください。