[bzoj3175] [Tjoi2013]攻击装置（最大独立集）

3175: [Tjoi2013]攻击装置

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 1129  Solved: 551
[Submit][Status][Discuss]

Description

给定一个01矩阵，其中你可以在0的位置放置攻击装置。每一个攻击装置（x,y）都可以按照“日”字攻击其周围的 8个位置（x-1,y-2）,(x-2,y-1),(x+1,y-2),(x+2,y-1),(x-1,y+2),(x-2,y+1), (x+1,y+2),(x+2,y+1)
求在装置互不攻击的情况下，最多可以放置多少个装置。

Input

第一行一个整数N，表示矩阵大小为N*N。接下来N行每一行一个长度N的01串，表示矩阵。

Output

一个整数，表示在装置互不攻击的情况下最多可以放置多少个装置。

Sample Input

3
010
000
100

Sample Output

4

HINT

100%数据 N<=200

Source

题解

二分图的最大独立集 和上篇博客 bzoj2510连边方法类似

一个点向他能到的八个点连边

要注意的是范围N<=200 也就是说最多40000个点

显然用邻接矩阵存不了 这个时候可以尝试用vector 

当然这不是正解 用前向星存是最好的 但是数据水。。就过了 其实是卡过的。。

最大独立集=顶点数-最大匹配数

最后因为双向边 ans要除2

#include<iostream>#include<cstdio>#include<cstdlib>#include<cstring>#include<algorithm>#include<cmath>#include<queue>#include<vector>using namespace std;int n;int dx[10]={0,-1,-2,1,2,-1,-2,1,2};int dy[10]={0,-2,-1,-2,-1,2,1,2,1};//应该用前向星去存 但是由于数据弱还是过了 如果把上面两行换成下面的就是超时 可见顺序还是很重要啊 //int dx[10]={0,2,2,-2,-2,1,1,-1,-1};  //int dy[10]={0,1,-1,1,-1,2,-2,2,-2};int t[205][205],mat[40005];int flag[40005];vector<int> mp[40005];bool find(int x){for(int i=mp[x].size()-1;i>=0;i--){int v=mp[x][i];if(!flag[v]){flag[v]=1;if(!mat[v]||find(mat[v])){mat[v]=x;return 1;}}}return 0;}int mk(int x,int y){return (x-1)*n+y;}int main(){//freopen("rand.txt","r",stdin);//freopen("my.txt","w",stdout);scanf("%d",&n);for(int i=1;i<=n;i++){char ch[205];scanf("%s",ch+1);for(int j=1;j<=n;j++)if (ch[j]=='0') t[i][j]=0;else t[i][j]=1;}for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)if(!t[i][j])for(int k=1;k<=8;k++){int x=i+dx[k],y=j+dy[k];if(x<1||x>n||y<1||y>n||t[x][y]) continue;mp[mk(i,j)].push_back(mk(x,y));}int ans=0,tot=0;for(int i=1;i<=n;i++)for(int j=1;j<=n;j++){if(t[i][j]==1) continue; int k=mk(i,j);tot++;memset(flag,0,sizeof(flag));if(find(k)) ans++;//printf("when k=%d,ans=%d\n",k,ans);}//printf("tot=%d ans=%d\n",tot,ans);printf("%d",tot-ans/2);//双向边 一组算了两次 要/2 return 0;}

唉就不应该用vector 

调了好久发现顺序问题。。可能有组数据卡了之前的dx,dy顺序吧

tle了好多次

18:35:112128787 xlj3175Accepted3392 kb8724 msC++/Edit1587 B2017-06-23

 18:34:102128780 xlj3175Wrong_Answer2232 kb4 msC++/Edit1585 B2017-06-23

 18:29:582128770 xlj3175Time_Limit_Exceed3180 kb11828 msC++/Edit1297 B2017-06-23

 18:13:122128768 xlj3175Time_Limit_Exceed3124 kb11140 msC++/Edit1436 B2017-06-23

 18:10:412128755 xlj3175Time_Limit_Exceed3140 kb11632 msC++/Edit1412 B2017-06-23 

16:49:162128643 xlj3175Time_Limit_Exceed3136 kb11260 msC++/Edit1401 B2017-06-23

 16:45:592128641 xlj3175Runtime_Error2872 kb32 msC++/Edit1400 B2017-06-23 16:44:56