BZOJ 1746: [Usaco2005 open]Lazy Cows

1746

一道思维复杂度不高...代码细节很多的dp题

这题突破口显然就在于只有2行

那么我们考虑第i列的状态

共有五种：

1. 没有纸片覆盖

2.上边被覆盖

3.下边被覆盖

4.被同一纸片覆盖

5.被不同的纸片覆盖

离散后，显然不存在第一种状态

那么我们用dp[i][t]表示第i个有宝藏的列  状态为t的最小面积

于是就可以转移了

实现细节见代码……

（因为改来改去导致代码比较难看

#include<iostream>
#include<cstdio>
#include<queue>
#include<cmath>
#include<algorithm>
#include<map>
#include<cstring>
#define eps 1e-5
#define inf 1e9
#define For(i,j,k)  for(int i=j;i<=k;++i)
#define Dow(i,j,k)  for(int i=k;i>=j;--i)
#define ll long long
#define maxn 1000000
using namespace std;
inline ll read()
{
    ll t=0,f=1;char c=getchar();
    while(c<'0'||c>'9')   {if(c=='-') f=-1;c=getchar();}
    while(c<='9'&&c>='0') t=t*10LL+c-48LL,c=getchar();
    return t*f;
}
struct node{int x,l;}   q[2001],a[2001];
inline bool cmp(node x,node y){return x.x<y.x;}
//1 ->  *X
//2 ->  X*
//3 ->  **
//4 ->  *#
int dp[1100][1100][5],n,k,m,tot;
int main()
{
    n=read();k=read();m=read();
    For(i,1,n)  a[i].l=read(),a[i].x=read();
    sort(a+1,a+n+1,cmp);
    For(i,1,n)  if(a[i].x!=a[i-1].x)    q[++tot].x=a[i].x,q[tot].l=a[i].l;else  if(a[i].l+a[i-1].l==3)  q[tot].l=3;
    n=tot;
    For(i,1,n)  For(j,1,k)  For(t,1,4)  dp[i][j][t]=inf;
    if(q[1].l==1)   dp[1][1][1]=1;else if(q[1].l==2)    dp[1][1][2]=1;
    dp[1][2][4]=2;
    dp[1][1][3]=2;
    For(i,2,n)
    {
        For(j,1,k)
        {
            if(q[i].l==1)   
                dp[i][j][1]=min(dp[i-1][j][4],dp[i-1][j][1])+q[i].x-q[i-1].x;
            if(q[i].l==2)   
                dp[i][j][2]=min(dp[i-1][j][4],dp[i-1][j][2])+q[i].x-q[i-1].x;
            if(j!=1)    dp[i][j][4]=dp[i-1][j][4]+2*(q[i].x-q[i-1].x);
            dp[i][j][3]=dp[i-1][j][3]+2*(q[i].x-q[i-1].x);
            if(j==1)    continue;
            int mi=min(min(dp[i-1][j-1][1],dp[i-1][j-1][2]),min(dp[i-1][j-1][3],dp[i-1][j-1][4]));
            if(q[i].l==1)dp[i][j][1]=min(dp[i][j][1],mi+1);
            if(q[i].l==2)dp[i][j][2]=min(dp[i][j][2],mi+1);
            dp[i][j][3]=min(dp[i][j][3],mi+2);
            mi=min(dp[i-1][j-1][1],min(dp[i-1][j-1][2],dp[i-1][j-1][4]))+q[i].x-q[i-1].x+1;
            dp[i][j][4]=min(dp[i][j][4],mi);
            if(j==2)    continue;
            mi=min(min(dp[i-1][j-2][1],dp[i-1][j-2][2]),min(dp[i-1][j-2][3],dp[i-1][j-2][4]));
            dp[i][j][4]=min(dp[i][j][4],mi+2);
        }
    }
    int ans=1e9;
    For(i,1,k)  For(t,1,4)  ans=min(ans,dp[n][i][t]);
    printf("%d\n",ans);
}