[COGS746] [网络流24题] 骑士共存
Posted 鄉勇
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[COGS746] [网络流24题] 骑士共存相关的知识,希望对你有一定的参考价值。
★★☆ 输入文件:knight.in
输出文件:knight.out
简单对比
时间限制:1 s
内存限制:128 MB
骑士共存问题
«问题描述:
在一个n*n个方格的国际象棋棋盘上,马(骑士)可以攻击的棋盘方格如图所示。棋盘
上某些方格设置了障碍,骑士不得进入。
«编程任务:
对于给定的n*n个方格的国际象棋棋盘和障碍标志,计算棋盘上最多可以放置多少个骑
士,使得它们彼此互不攻击。
«数据输入:
由文件knight.in给出输入数据。第一行有2 个正整数n 和m (1<=n<=200, 0<=m<=n*n)
分别表示棋盘的大小和障碍数。接下来的m 行给出障碍的位置。每行2 个正整数,表示障
碍的方格坐标。
«结果输出:
将计算出的共存骑士数输出到文件knight.out。
输入文件示例 输出文件示例
knight.in
3 2
1 1
3 3
knight.out
思路
网络流
http://www.cnblogs.com/shenben/p/6259661.html
代码实现
1 #include<cstdio> 2 #include<cstring> 3 const int maxn=800; 4 const int maxs=maxn*maxn; 5 inline int min_(int x,int y){return x<y?x:y;} 6 int n,m,s,t,ans; 7 int a,b,nh,nl; 8 int map[maxn][maxn]; 9 int hb[]={-2,-1,1,2,2,1,-1,-2}; 10 int lb[]={1,2,2,1,-1,-2,-2,-1}; 11 int h[maxs],hs=1,ew[maxs],et[maxs],en[maxs]; 12 void add(int u,int v){ 13 ++hs,ew[hs]=1,et[hs]=v,en[hs]=h[u],h[u]=hs; 14 ++hs,ew[hs]=0,et[hs]=u,en[hs]=h[v],h[v]=hs; 15 } 16 int q[maxs],head,tail; 17 int d[maxs]; 18 void bfs(){ 19 memset(d,0,sizeof(d)); 20 head=tail=0; 21 q[head++]=s,d[s]=1; 22 while(head>tail){ 23 a=q[tail++]; 24 for(int i=h[a];i;i=en[i]) 25 if(!d[et[i]]&&ew[i]){ 26 d[et[i]]=d[a]+1; 27 if(et[i]==t) return; 28 q[head++]=et[i]; 29 } 30 } 31 } 32 int ap(int k,int w){ 33 if(k==t) return w; 34 int dw=w; 35 for(int i=h[k];i&&dw;i=en[i]) 36 if(ew[i]&&d[et[i]]==d[k]+1){ 37 int wt=ap(et[i],min_(dw,ew[i])); 38 if(wt) ew[i]-=wt,ew[i^1]+=wt,dw-=wt; 39 else d[et[i]]=0; 40 } 41 return w-dw; 42 } 43 void Dinic(){while(bfs(),d[t]) ans-=ap(s,ans);} 44 int main(){ 45 freopen("knight.in","r",stdin); 46 freopen("knight.out","w",stdout); 47 scanf("%d%d",&n,&m); 48 s=0,t=1,ans=n*n-m; 49 for(int i=1;i<=n;i++) 50 for(int j=1;j<=n;j++) 51 map[i][j]=i&1^j&1; 52 for(int i=1;i<=m;i++){ 53 scanf("%d%d",&a,&b); 54 map[a][b]=2; 55 } 56 for(int i=1;i<=n;i++) 57 for(int j=1;j<=n;j++) 58 if(!map[i][j]){ 59 add(s,i*n+j); 60 for(int k=0;k<8;k++){ 61 nh=i+hb[k],nl=j+lb[k]; 62 if(nh>0&&nh<=n&&nl>0&&nl<=n&&map[nh][nl]==1) 63 add(i*n+j,nh*n+nl); 64 } 65 } 66 else if(map[i][j]==1) add(i*n+j,t); 67 Dinic(); 68 printf("%d\\n",ans); 69 return 0; 70 }
以上是关于[COGS746] [网络流24题] 骑士共存的主要内容,如果未能解决你的问题,请参考以下文章