矩形分割（動態規劃）

 矩形分割（cut）
Description
出於某些方面的需求，我們要把一塊N×M的木板切成一個個1×1的小方塊。
對於一塊木板，我們只能從某條橫線或者某條豎線（要在方格線上），而且這木板是不均勻的，從不同的線切割下去要花不同的代價。而且，對於一塊木板，切割一次以後就被分割成兩塊，而且不能把這兩塊木板拼在一起然後一刀切成四塊，只能兩塊分別再進行一次切割。
現在，給出從不同的線切割所要花的代價，求把整塊木板分割成1×1塊小方塊所需要耗費的最小代價。
Input Format
輸入文件第一行包括N和M，表示長N寬M的矩陣。
第二行包括N-1個非負整數，分別表示沿着N-1條橫線切割的代價。
第二行包括M-1個非負整數，分別表示沿着M-1條豎線切割的代價。

=====================My code=======================

#include<fstream> using namespace std; ifstream fin("cut.in"); ofstream fout("cut.out"); const int N=2001,M=2001; int ans[M][N];//1維m2維n int hang[N],lie[M]; int quick_sort(int a[],int low,int high){ int i=low,j=high,temp=a[j]; if (i>=j)return 0; while(i<j){ while(i<j&&a[i]>=temp)i++; if (i<j)a[j--]=a[i]; while(i<j&&a[j]<=temp)j--; if (i<j)a[i++]=a[j]; } a[i]=temp; quick_sort(a,i+1,high); quick_sort(a,low,i-1); return 0; } int dp(int a[M][N],int n,int m){ int i,j; for (i=1;i<=m;i++) for (j=1;j<=n;j++) { a[i][j]=min(a[i-1][j]+lie[i]*(j+1),a[i][j-1]+hang[j]*(i+1)); } return a[m][n]; } int work1(int a[M][N],int n,int m){ int i=1; int s=0; quick_sort(hang,1,n); quick_sort(lie,1,m); for (i=1;i<=m;i++) { s+=lie[i]; a[i][0]=s; } s=0; for (i=1;i<=n;i++) { s+=hang[i]; a[0][i]=s; } return 0; } int main(){ int n,m; int i,j; fin>>n>>m; n--; m--; for (i=1;i<=n;i++) { fin>>hang[i]; } for (i=1;i<=m;i++) { fin>>lie[i]; } work1(ans,n,m); fout<<dp(ans,n,m); return 0; }