[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