题意:
两行格子,每个格子是1或者0,共有n个1,现在要求使用k个矩形覆盖所有的1,求出最小矩形面积(总和)…
分析:
每一列就2个格子,妥妥的状压dp,因为n很小b(列数)很大…所以需要离散化…
f[i][j][k]代表前i列使用j个矩形覆盖最后一列状态为k的最小矩形面积…
k=0代表只覆盖上面的格子
k=1代表只福噶下面的格子
k=2代表分别覆盖上下两个格子
k=3代表同时覆盖上下两个格子
转移比较麻烦,看代码吧…
代码如下:
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
//by NeighThorn
#define inf 0x3f3f3f3f
#define Min(a,b,c) min(a,min(b,c))
using namespace std;
const int maxn=1000+5;
int n,k,b,mp[maxn],cnt,vis[3][maxn],pos[maxn],f[maxn][maxn][4];
struct M{
int x,y;
}s[maxn];
signed main(void){
scanf("%d%d%d",&n,&k,&b);cnt=0;
for(int i=1,x,y;i<=n;i++)
scanf("%d%d",&s[i].x,&s[i].y),mp[++cnt]=s[i].y;
// cout<<"lala"<<endl;
sort(mp+1,mp+cnt+1);memset(vis,0,sizeof(vis));
int len=unique(mp+1,mp+cnt+1)-mp-1;
for(int i=1;i<=n;i++)
pos[lower_bound(mp+1,mp+len+1,s[i].y)-mp]=s[i].y,vis[s[i].x][lower_bound(mp+1,mp+len+1,s[i].y)-mp]=1;
memset(f,inf,sizeof(f));
f[1][1][0]=1,f[1][1][1]=1;
f[1][2][2]=2,f[1][1][3]=2;
if(vis[1][1])
f[1][1][1]=inf;
if(vis[2][1])
f[1][1][0]=inf;
// for(int i=1;i<=len;i++)
// cout<<pos[i]<<" "<<vis[1][i]<<" "<<vis[2][i]<<endl;
// cout<<endl;
for(int i=2;i<=len;i++)
for(int j=1;j<=k;j++){
// cout<<i<<" "<<j<<endl;
if(!vis[1][i])
f[i][j][1]=Min(f[i][j][1],f[i-1][j][1]+pos[i]-pos[i-1],f[i-1][j-1][1]+1),
f[i][j][1]=Min(f[i][j][1],f[i-1][j-1][0]+1,f[i-1][j-1][2]+1),
f[i][j][1]=Min(f[i][j][1],f[i-1][j-1][3]+1,f[i-1][j][2]+pos[i]-pos[i-1]);
if(!vis[2][i])
f[i][j][0]=Min(f[i][j][0],f[i-1][j][0]+pos[i]-pos[i-1],f[i-1][j-1][0]+1),
f[i][j][0]=Min(f[i][j][0],f[i-1][j-1][1]+1,f[i-1][j-1][2]+1),
f[i][j][0]=Min(f[i][j][0],f[i-1][j-1][3]+1,f[i-1][j][2]+pos[i]-pos[i-1]);
f[i][j][2]=Min(f[i][j][2],j-2<0?inf:f[i-1][j-2][2]+2,f[i-1][j-1][2]+1+pos[i]-pos[i-1]),
f[i][j][2]=Min(f[i][j][2],j-2<0?inf:f[i-1][j-2][3]+2,f[i-1][j][2]+2*(pos[i]-pos[i-1])),
f[i][j][2]=Min(f[i][j][2],j-2<0?inf:f[i-1][j-2][0]+2,f[i-1][j-1][0]+pos[i]-pos[i-1]+1),
f[i][j][2]=Min(f[i][j][2],j-2<0?inf:f[i-1][j-2][1]+2,f[i-1][j-1][1]+pos[i]-pos[i-1]+1);
f[i][j][3]=Min(f[i][j][3],f[i-1][j-1][3]+2,f[i-1][j][3]+2*(pos[i]-pos[i-1])),//cout<<f[i][j][3]<<" ",
f[i][j][3]=min(f[i][j][3],f[i-1][j-1][2]+2),//cout<<f[i][j][3]<<" ",
f[i][j][3]=Min(f[i][j][3],f[i-1][j-1][0]+2,f[i-1][j-1][1]+2);//cout<<f[i][j][3]<<endl;
}
// for(int i=1;i<=len;i++,cout<<endl)
// for(int j=1;j<=k;j++,cout<<endl)
// for(int k=0;k<4;k++)
// cout<<f[i][j][k]<<" ";
int ans=0x3f3f3f3f;
for(int i=1;i<=k;i++)
for(int j=0;j<4;j++)
ans=min(ans,f[len][i][j]);
cout<<ans<<endl;
return 0;
}
by >_< NeighThorn