HDU4739Zhuge Liang's Mines(状压）

链接

预处理出只有四个1的情况存入数组中 然后状压下

#include <iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<stdlib.h>
using namespace std;
#define N 1050000
int dp[2][N];
int q[2][N];
int pt[N];
bool f[N];
struct node
{
    int x,y;
}pp[25];
bool cmp(node a,node b)
{
    if(a.x==b.x)
    return a.y<b.y;
    return a.x<b.x;
}
int main()
{
    int n,i,j,k,o,p[25][2];
    while(scanf("%d",&n)!=EOF)
    {
        if(n==-1)
        break;
        memset(dp,0,sizeof(dp));
        for(i = 0; i < n ; i++)
        scanf("%d%d",&pp[i].x,&pp[i].y);
        sort(pp,pp+n,cmp);
        k = 0;
        for(i = 0 ; i < 1<<n ; i++)
        {
            o = 0;
            for(j = 0 ; j < n ; j++)
            {
                if(i&(1<<j))
                {
                    o++;
                    p[o][0] = pp[j].x;
                    p[o][1] = pp[j].y;
                }
            }
            if(o==4)
            {
                int x1,x2;
                if(p[1][0]==p[2][0]&&p[3][0]==p[4][0]&&p[2][1]==p[4][1]&&p[1][1]==p[3][1])
                {
                    x1 = p[2][1]-p[1][1];
                    x2 = p[3][0]-p[1][0];
                    if(x1==x2)
                    {
                        k++;
                        q[1][k] = i;
                        dp[1][i] = 4;
                        pt[k] = i;
                    }
                }
            }
        }
        int tt = k;
        o=k;
        for(i = 2; i <= n/4+1 ; i++)
        {
            k = tt;tt=0;
            for(j = 1 ; j <= k ; j++)
            {
                for(int g = 1 ; g <= o ; g++)
                {
                    if(q[(i-1)%2][j]&pt[g])
                    continue;
                    int x = q[(i-1)%2][j]+pt[g];
                    dp[i%2][x] = max(dp[i%2][x],dp[(i-1)%2][q[(i-1)%2][j]]+4);
                    if(!f[x])
                    {
                        f[x] = 1;
                        tt++;
                        q[i%2][tt] = x;
                    }
                }
            }
        }
        int ans=0;
        for(i = 0 ; i < 1<<n ; i++)
        {
            ans = max(ans,max(dp[0][i],dp[1][i]));
        }
        printf("%d\n",ans);
    }
    return 0;
}