BZOJ 1066 POJ 2711 [SCOI2007]蜥蜴

1066: [SCOI2007]蜥蜴

Time Limit: 1 Sec Memory Limit: 162 MB
Submit: 2947 Solved: 1471
[Submit][Status][Discuss]

Description

　　在一个r行c列的网格地图中有一些高度不同的石柱，一些石柱上站着一些蜥蜴，你的任务是让尽量多的蜥蜴逃
到边界外。每行每列中相邻石柱的距离为1，蜥蜴的跳跃距离是d，即蜥蜴可以跳到平面距离不超过d的任何一个石
柱上。石柱都不稳定，每次当蜥蜴跳跃时，所离开的石柱高度减1（如果仍然落在地图内部，则到达的石柱高度不
变），如果该石柱原来高度为1，则蜥蜴离开后消失。以后其他蜥蜴不能落脚。任何时刻不能有两只蜥蜴在同一个
石柱上。

Input

　　输入第一行为三个整数r，c，d，即地图的规模与最大跳跃距离。以下r行为石竹的初始状态，0表示没有石柱
，1~3表示石柱的初始高度。以下r行为蜥蜴位置，“L”表示蜥蜴，“.”表示没有蜥蜴。

Output

　　输出仅一行，包含一个整数，即无法逃离的蜥蜴总数的最小值。

Sample Input

5 8 2
00000000
02000000
00321100
02000000
00000000
........
........
..LLLL..
........
........

Sample Output

HINT

100%的数据满足：1<=r, c<=20, 1<=d<=4

从隔壁Orion_Rigel博客搬了张思路图过来~~（我是盗图王）~~

感觉代码规模差不多，为何为TLE呢？

先放TLE的代码，以后AC了再更新

————UPD2017.5

一拖就拖了一年的节奏？

其实只要换成邻接表就可以了

  1 #include<iostream>
  2 #include<cstdio>
  3 #include<cstring>
  4 #include<algorithm>
  5 #include<queue>
  6 using namespace std;
  7 const int INF=0x3f3f3f3f;
  8 const int mxn=1205;
  9 struct edge{
 10     int v,nxt,f;
 11 }e[mxn*mxn];
 12 int hd[mxn],mct=1;
 13 void add_edge(int u,int v,int f){
 14     e[++mct].v=v;e[mct].nxt=hd[u];e[mct].f=f;hd[u]=mct;return;
 15 }
 16 void insert(int u,int v,int f){
 17     if(!v)return;
 18 //    printf("%d to %d : %d
",u,v,f);
 19     add_edge(u,v,f); add_edge(v,u,0);
 20     return;
 21 }
 22 int n,m,D,S,T;
 23 int d[mxn];
 24 bool BFS(){
 25     queue<int>q;
 26     memset(d,0,sizeof d);
 27     d[S]=1;q.push(S);
 28     while(!q.empty()){
 29         int u=q.front();q.pop();
 30         for(int i=hd[u];i;i=e[i].nxt){
 31             int v=e[i].v;
 32             if(!d[v] && e[i].f){
 33 //                printf("%d to %d :%d
",u,v,e[i].f);
 34                 d[v]=d[u]+1;
 35                 q.push(v);
 36             }
 37         }
 38     }
 39     return d[T];
 40 }
 41 int DFS(int u,int lim){
 42     if(u==T)return lim;
 43     int f=0,tmp;
 44     for(int i=hd[u],v;i;i=e[i].nxt){
 45         v=e[i].v;
 46         if(e[i].f && d[v]==d[u]+1 && (tmp=DFS(e[i].v,min(lim,e[i].f)))){
 47             e[i].f -=tmp;
 48             e[i^1].f+=tmp;
 49             f+=tmp;
 50             lim -=tmp;
 51             if(!lim)return f;
 52         }
 53     }
 54     d[u]=0;
 55     return f;
 56 }
 57 int Dinic(){
 58     int res=0;
 59     while(BFS())res+=DFS(S,INF);//,printf("res:%d
",res);
 60     return res;
 61 }
 62 int id[25][25],sz=0,lct=0;
 63 int mp[25][25];
 64 void ADD(int x,int y,int out){
 65     bool flag=0;
 66     for(int i=x-D;i<=x+D;i++){
 67         for(int j=y-D;j<=y+D;j++){
 68             if(i==x && j==y)continue;
 69             if( (i-x)*(i-x)+(j-y)*(j-y)<=D*D ){
 70                 if(i<1 || i>n || j<1 || j>m){if(!flag)insert(out,T,INF),flag=1;}
 71                 else insert(out,id[i][j],INF);
 72             }
 73         }
 74     }
 75     return;
 76 }
 77 void Build(){
 78     T=sz+sz+lct+1;
 79     for(int i=1;i<=n;i++)
 80         for(int j=1;j<=m;j++){
 81             int out=id[i][j]+lct+sz;
 82             if(mp[i][j]){
 83 //                printf("[%d %d ]%d
",i,j,mp[i][j]);
 84                 insert(id[i][j],out,mp[i][j]);    
 85                 ADD(i,j,out);
 86             }
 87         }
 88     return;
 89 }
 90 char sc[25];
 91 int main(){
 92 //    freopen("in.txt","r",stdin);
 93     int i,j;
 94     scanf("%d%d%d",&n,&m,&D);
 95     for(i=1;i<=n;i++){
 96         scanf("%s",sc+1);
 97         for(j=1;j<=m;j++){
 98             mp[i][j]=sc[j]-'0';
 99             if(mp[i][j])id[i][j]=++sz;
100         }
101     }
102     S=0;
103     for(i=1;i<=n;i++){
104         scanf("%s",sc+1);
105         for(j=1;j<=m;j++){
106             if(sc[j]=='L'){
107                 ++lct;
108                 insert(S,lct+sz,1);
109                 insert(lct+sz,id[i][j],1);
110             }
111         }
112     }
113     Build();
114     int ans=Dinic();
115 //    printf("%d
",ans);
116 /*    for(i=hd[2];i;i=e[i].nxt){
117         int v=e[i].v;
118         printf("%d %d : %d
",2,v,e[i].f);
119     }*/
120     printf("%d
",lct-ans);
121     return 0;
122 }

  1 #include<iostream>
  2 #include<cstdio>
  3 #include<algorithm>
  4 #include<cmath>
  5 #include<cstring>
  6 using namespace std;
  7 const int INF=0x6ffffff;
  8 struct NODE{
  9     int w;//容量(在此题中为石柱高度)
 10     int f;//流量 
 11 }e[1000][1000];//邻接矩阵流量边 
 12 int r,c;
 13 int s,t;
 14 int dis;//跳跃距离 
 15 int mp[60][60],//坐标为[x][y]的点是否为结点，是的话记录节点号
 16     ht[60][60];//高度
 17 int cnode=0,clizard=0;//结点数 蜥蜴数 
 18 //
 19 int add_eg(int u,int v,int h){
 20     e[u][v].w=h;
 21 }
 22 int pd(int x1,int y1,int x2,int y2){
 23     return (x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)<=dis*dis;
 24 }
 25 //
 26 int q[6000],d[60];
 27 bool flag[60];
 28 int BFS(){
 29     int hd=1,tl=1;
 30     int u,v;
 31     memset(flag,0,sizeof(flag));
 32     q[hd]=s;d[s]=0;flag[s]=1;
 33     while(hd<=tl){
 34         u=q[hd];
 35         for(v=0;v<=t;v++){
 36             if(!flag[v] && e[u][v].w>e[u][v].f){
 37                 flag[v]=1;
 38                 d[v]=d[u]+1;
 39                 q[++tl]=v;
 40             }
 41             if(flag[t])return 1;
 42         }
 43         hd++;
 44     }
 45     return 0;
 46 
 47 }
 48 int DFS(int v,int low){
 49     int i;    int flow;
 50     if(v==t)return low;
 51     for(i=0;i<=t;i++){
 52         if(e[v][i].w>e[v][i].f && d[i]==d[v]+1){
 53             if(flow=DFS(i,min(low,e[v][i].w-e[v][i].f))){
 54                 e[v][i].f+=flow;
 55                 e[i][v].f=-e[v][i].f;
 56                 return flow;
 57             }
 58         }
 59     }
 60     return 0;
 61 }
 62 int dinic(){
 63     int ans=0;
 64     while(BFS()){
 65         int flow;
 66         while(flow=DFS(s,INF)){
 67             ans+=flow;
 68         }
 69     }
 70     return ans;
 71 }
 72 int main(){
 73     scanf("%d%d%d",&r,&c,&dis);
 74     int i,j;
 75     char ch[30];
 76     for(i=1;i<=r;i++){//读石柱图 
 77         scanf("%s",ch);
 78         for(j=0;j<c;j++){
 79             ht[i][j+1]=ch[j]-'0';
 80             if(ht[i][j+1]>0)mp[i][j+1]=++cnode;
 81         }
 82     }
 83     s=0;t=cnode*2+1;
 84     for(i=1;i<=r;i++){//读蜥蜴图 
 85         scanf("%s",ch);
 86         for(j=0;j<c;j++){
 87             if(ch[j]=='L'){
 88                 add_eg(s,mp[i][j+1],1);
 89                 clizard++;
 90             }
 91         }
 92     }
 93     for(i=1;i<=r;i++){
 94         for(j=1;j<=c;j++){
 95             
 96             if(mp[i][j]){
 97             
 98                 //建石柱自身到自身镜像的容量为石柱高度的边 
 99                 add_eg(mp[i][j],mp[i][j]+cnode,ht[i][j]);
100                 
101             if(i<=dis || i+dis>r || j<=dis || j+dis>c)
102                 add_eg(mp[i][j]+cnode,t,INF);//建出区域的边 
103             int a,b;
104             for(a=max(i-dis,1);a<=min(i+dis,r);a++)
105               for(b=max(j-dis,1);b<=min(j+dis,c);b++){
106                   if(i!=a || b!=j)
107                   if(pd(i,j,a,b)) add_eg(mp[i][j]+cnode,mp[a][b],INF);//建石柱之间的边 
108               }
109             }
110         }
111     }
112     int ans=dinic();
113 //    printf("%d %d",clizard,ans);
114     printf("%d
",clizard-ans);
115     return 0;
116 }

TLE