棋盘分割
終于過了。
題目:http://poj.org/problem?id=1191
dp,水平有限,想不清更新順序,所以打的記憶化,
f[n,x,y,xx,yy]表示還剩下n塊,以及還剩下的這部分的左上角右下角坐標,所存儲的最優值
需要對公式變形為n*σ*σ=Σ(xi-x )^2,維護的即為各項的和
代碼:
program sky; constmaxn = 10000000; vari,j,n,x:longint;ave:extended;a:array[0..8,0..8] of longint;f:array[0..16,0..9,0..9,0..9,0..9] of extended; function min(qq,ww:extended):extended; beginif qq>ww then exit(ww); exit(qq); end; function get(x,y,xx,yy:longint):extended; vartp:longint; begintp:=a[xx,yy]-a[x-1,yy]-a[xx,y-1]+a[x-1,y-1];get:=sqr(tp-ave); end; function dp(n,x,y,xx,yy:longint):extended; vari:longint; beginif f[n,x,y,xx,yy]<>0 then exit(f[n,x,y,xx,yy]);if n=1 thenbeginf[n,x,y,xx,yy]:=get(x,y,xx,yy);exit(f[n,x,y,xx,yy]);end;f[n,x,y,xx,yy]:=maxn;for i:=x to xx-1 do{題目的限制只用遞歸一邊,另一邊不再切割}begin{一開始加了一些限制條件,比如(i-x+1>=n)才搜,后來發現不對}{因為加上會丟狀態,有可能橫著切不了豎著可以切}f[n,x,y,xx,yy]:=min(f[n,x,y,xx,yy],dp(n-1,x,y,i,yy)+dp(1,i+1,y,xx,yy));f[n,x,y,xx,yy]:=min(f[n,x,y,xx,yy],dp(1,x,y,i,yy)+dp(n-1,i+1,y,xx,yy));end;for i:=y to yy-1 dobeginf[n,x,y,xx,yy]:=min(f[n,x,y,xx,yy],dp(n-1,x,y,xx,i)+dp(1,x,i+1,xx,yy));f[n,x,y,xx,yy]:=min(f[n,x,y,xx,yy],dp(1,x,y,xx,i)+dp(n-1,x,i+1,xx,yy));end;exit(f[n,x,y,xx,yy]); end; beginassign(input,'fgqp.in'); reset(input);assign(output,'fgqp.out'); rewrite(output);readln(n);for i:=1 to 8 dofor j:=1 to 8 dobeginread(x);a[i,j]:=a[i-1,j]+a[i,j-1]-a[i-1,j-1]+x;end;ave:=a[8,8]/n;writeln(sqrt(dp(n,1,1,8,8)/n):0:3);close(input); close(output); end.轉載于:https://www.cnblogs.com/skysun/archive/2012/05/24/2516286.html
總結
- 上一篇: PostgreSQL学习手册(角色和权限
- 下一篇: String 类详解