poj2836 状态压缩~矩形覆盖

最新推荐文章于 2019-04-29 21:08:40 发布

原创最新推荐文章于 2019-04-29 21:08:40 发布 · 1.4k 阅读

1 ·

CC 4.0 BY-SA版权

动态规划专栏收录该内容

76 篇文章

订阅专栏

这篇博客介绍了如何利用状态压缩的方法解决POJ2836问题，通过预处理点集关系，进行动态规划求解，最终实现代码AC。

首先预处理每两个点为顶点包含的点集合，然后往后直接DP。

ACcode：

#include<cstdio>
#include<cstring>

const int NS=15;
const int INF=1<<30;

int n,lx,rx,uy,dy,lim;
int g[NS][NS],ar[NS][NS];
int x[NS],y[NS],dp[1<<NS];

int Max(int xa,int xb){
    return xa>xb?xa:xb;
}

int Min(int xa,int xb){
    return xa<xb?xa:xb;
}

int main()
{
    while (~scanf("%d",&n)&&n)
    {
        for (int i=0;i<n;i++)
        scanf("%d%d",&x[i],&y[i]);
        for (int i=0;i<n;i++)
        for (int j=0;j<i;j++)
        {
            g[i][j]=0;
            lx=Min(x[i],x[j]),rx=Max(x[i],x[j]);
            dy=Min(y[i],y[j]),uy=Max(y[i],y[j]);
            for (int k=0;k<n;k++)
            if (lx<=x[k]&&x[k]<=rx&&dy<=y[k]&&y[k]<=uy)
            g[i][j]|=1<<k;
            rx=rx-lx?rx-lx:1; uy=uy-dy?uy-dy:1;
            ar[i][j]=ar[j][i]=rx*uy;
        }

        lim=1<<n,dp[0]=0;
        for (int i=1;i<lim;i++)
        dp[i]=INF;
        for (int i=0;i<lim;i++)
        for (int j=0;j<n;j++)
        for (int k=0;k<j;k++)
        {
            lx=i|g[j][k];
            dp[lx]=Min(dp[lx],dp[i]+ar[j][k]);
        }
        printf("%d\n",dp[lim-1]);
    }
    return 0;
}