POJ 2430 Lazy Cows 状压DP

时空隧道

---

题意：
两行格子，每个格子是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