[SCOI2005][BZOJ 1084]最大子矩陣
阿新 • • 發佈:2018-05-15
def bits its 轉載 ace cloud pap div pre
Description
這裏有一個n*m的矩陣,請你選出其中k個子矩陣,使得這個k個子矩陣分值之和最大。註意:選出的k個子矩陣
不能相互重疊。
Input
第一行為n,m,k(1≤n≤100,1≤m≤2,1≤k≤10),接下來n行描述矩陣每行中的每個元素的分值(每個元素的
分值的絕對值不超過32767)。
Output
只有一行為k個子矩陣分值之和最大為多少。
題解: 看了半天,突然發現,m小於等於2啊。 然後就亂dp一波,除了轉移寫起來很麻之外,就沒什麽了。 令f[i][j][k]表示當前第i行,以選中j個矩陣,當前行的取法為k的得分數,(取法只有5種啦)。 代碼:#include<bits/stdc++.h> using namespace std; inline int read(){ int x=0,f=1;char ch=getchar(); while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();} while(ch>=‘0‘&&ch<=‘9‘){x=(x<<3)+(x<<1)+ch-‘0‘;ch=getchar();} return x*f; } #define MN 105 #define inf 0x7f intn,m,k,a[MN][3],ans; int f[MN][15][6]; void rw(int &x,int y){if(y>x)x=y;} int main(){ n=read(),m=read(),k=read(); register int i,j; for(i=1;i<=n;i++)for(j=1;j<=m;j++) a[i][j]=read(); if(m==1){ for(i=1;i<=n;i++) for(j=1;j<=k;j++){ rw(f[i][j][0],f[i-1][j][0]); rw(f[i][j][0],f[i-1][j][1]); rw(f[i][j][1],f[i-1][j][1]+a[i][1]); rw(f[i][j][1],f[i-1][j-1][0]+a[i][1]); } printf("%d\n",max(f[n][k][1],f[n][k][0])); } else{ memset(f,-inf,sizeof f); for(i=0;i<=n;i++)for(j=0;j<=k;j++) f[i][j][1]=0; for(i=1;i<=n;i++)for(j=1;j<=k;j++){ for(int h=1;h<=5;h++) rw(f[i][j][1],f[i-1][j][h]); f[i][j][2]=max(f[i-1][j][2],f[i-1][j][5])+a[i][1]; f[i][j][3]=max(f[i-1][j][3],f[i-1][j][5])+a[i][2]; rw(f[i][j][2],f[i-1][j-1][4]+a[i][1]); rw(f[i][j][3],f[i-1][j-1][4]+a[i][2]); rw(f[i][j][2],max(f[i-1][j-1][1],f[i-1][j-1][3])+a[i][1]); rw(f[i][j][3],max(f[i-1][j-1][1],f[i-1][j-1][2])+a[i][2]); f[i][j][4]=max(f[i-1][j-1][1],f[i-1][j][4])+a[i][1]+a[i][2]; rw(f[i][j][4],max(f[i-1][j-1][2],f[i-1][j-1][3])+a[i][1]+a[i][2]); rw(f[i][j][4],f[i-1][j-1][5]+a[i][1]+a[i][2]); f[i][j][5]=f[i-1][j][5]+a[i][1]+a[i][2]; rw(f[i][j][5],max(f[i-1][j-1][2],f[i-1][j-1][3])+a[i][1]+a[i][2]); if(j>=2) rw(f[i][j][5],max(f[i-1][j-2][1],f[i-1][j-2][4])+a[i][1]+a[i][2]); } ans=-1e12; for(i=1;i<=5;i++) rw(ans,f[n][k][i]); printf("%d\n",ans); } return 0; }
來自PaperCloud的博客,未經允許,請勿轉載,TKS!
[SCOI2005][BZOJ 1084]最大子矩陣