【BZOJ1048】分割矩陣（記憶化搜索，動態規劃）

題面

BZOJ
洛谷

題解

一個很簡單的\(dp\)，寫成記憶化搜索的形式的挺不錯的。

#include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> #include<cmath> #include<algorithm> using namespace std; inline int read() {int x=0;bool t=false;char ch=getchar();while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();if(ch=='-')t=true,ch=getchar();while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();return t?-x:x; } double sqr(double x){return x*x;} int A,B,n; int g[11][11]; double f[11][11][11][11][11],avg; double dfs(int x1,int y1,int x2,int y2,int n) {if(n>(x2-x1+1)*(y2-y1+1))return 1e9;if(f[x1][y1][x2][y2][n]<1e9)return f[x1][y1][x2][y2][n];double ret=1e9;for(int a=1;a<n;++a){for(int i=x1;i<x2;++i)ret=min(ret,dfs(x1,y1,i,y2,a)+dfs(i+1,y1,x2,y2,n-a));for(int i=y1;i<y2;++i)ret=min(ret,dfs(x1,y1,x2,i,a)+dfs(x1,i+1,x2,y2,n-a));}return f[x1][y1][x2][y2][n]=ret; } int main() {A=read();B=read();n=read();for(int i=1;i<=A;++i)for(int j=1;j<=B;++j)avg+=(g[i][j]=read());for(int i=1;i<=A;++i)for(int j=1;j<=B;++j)g[i][j]+=g[i-1][j]+g[i][j-1]-g[i-1][j-1];avg/=n;memset(f,127,sizeof(f));for(int i=1;i<=A;++i)for(int j=1;j<=B;++j)for(int k=i;k<=A;++k)for(int l=j;l<=B;++l)f[i][j][k][l][1]=sqr(g[k][l]-g[i-1][l]-g[k][j-1]+g[i-1][j-1]-avg);printf("%.2lf\n",sqrt(dfs(1,1,A,B,n)/n));return 0; }

轉載于:https://www.cnblogs.com/cjyyb/p/9451575.html

總結

以上是生活随笔為你收集整理的【BZOJ1048】分割矩阵（记忆化搜索，动态规划）的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。