HDU 4739 Zhuge Liang's Mines（DP）

该题一开始没想出怎么暴力， 其实我们可以先预处理出来所有的可能正方形，存到一个vector里，然后用dp的思想就行了。

因为最多只有20个点，所以我们可以状态压缩一下，然后状态转移就是d[S] = max(ans,dp(ss)+4);   如果当前这个正方形所组成的点完全在S中，那么ss就是S去掉这四个点后的状态，如此转移就可以了。

细节参见代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<string>
#include<vector>
using namespace std;
typedef long long ll;
const int INF = 1000000000;
const int maxn = 23;
int T,n,kase=0,len,m,d[1<<maxn], vis[1<<maxn];
struct node {
int x, y;
node(int xx=0, int yy=0):x(xx), y(yy) {}
bool operator < (const node& rhs) const {
return x < rhs.x || (x == rhs.x && y < rhs.y);
}
}a[maxn];
vector<int> g;
int dp(int S) {
if(vis[S] == kase) return d[S];
vis[S] = kase;
int& ans = d[S];
ans = 0;
for(int i=0;i<len;i++) {
int v = g[i];
if((S|v) == S) {
int ss = S ^ v;
ans = max(ans,dp(ss)+4);
}
}
return ans;
}
int main() {
while(~scanf("%d",&n)) {
if(n < 0) return 0;
g.clear();
for(int i=0;i<n;i++) scanf("%d%d",&a[i].x,&a[i].y);
sort(a,a+n);
for(int i=0;i<n;i++) {
for(int j=i+1;j<n;j++) {
if(a[i].x != a[j].x) break;
for(int ii=j+1;ii<n;ii++) {
for(int jj=ii+1;jj<n;jj++) {
if(a[ii].x != a[jj].x) break;
if(a[i].y == a[ii].y && a[j].y == a[jj].y && a[j].y-a[i].y == a[ii].x-a[i].x) {
int v = 0;
v |= (1<<i); v |= (1<<j);
v |= (1<<ii); v |= (1<<jj);
g.push_back(v);
}
}
}
}
}
len = g.size();
++kase;
printf("%d\n",dp((1<<n)-1));
}
return 0;
}