ZOJ-2365 Strong Defence 贪心,BFS

  题目链接:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=2365

  我没看懂题目。。。这样理解:一个有向图,要给一些边染色,使得所用的颜色最多,且S到T的任意路径的都包含所有颜色。  或者: 给定一个无向图,图中有一个起点S和一个终点T。要求选K个集合S1,S2,…,SK,每个集合都含有图中的一些边,任意两个不同的集合的交集为空。并且从图中任意去掉一个集合,S到T都没有通路。要求K尽量大。

  容易想到最短距离就是总共的数量tot,然后在图上求一个层次图,只要点的距离标号i<=tot的时候点标记为i就可以了,如果大于tot,任意1-tot都可以标记。证明比较简单,不多说。。

  1 //STATUS:C++_AC_10MS_19492KB
  2 #include <functional>
  3 #include <algorithm>
  4 #include <iostream>
  5 //#include <ext/rope>
  6 #include <fstream>
  7 #include <sstream>
  8 #include <iomanip>
  9 #include <numeric>
 10 #include <cstring>
 11 #include <cassert>
 12 #include <cstdio>
 13 #include <string>
 14 #include <vector>
 15 #include <bitset>
 16 #include <queue>
 17 #include <stack>
 18 #include <cmath>
 19 #include <ctime>
 20 #include <list>
 21 #include <set>
 22 #include <map>
 23 using namespace std;
 24 //using namespace __gnu_cxx;
 25 //define
 26 #define pii pair<int,int>
 27 #define mem(a,b) memset(a,b,sizeof(a))
 28 #define lson l,mid,rt<<1
 29 #define rson mid+1,r,rt<<1|1
 30 #define PI acos(-1.0)
 31 //typedef
 32 //typedef __int64 LL;
 33 //typedef unsigned __int64 ULL;
 34 //const
 35 const int N=410;
 36 const int INF=0x3f3f3f3f;
 37 const int MOD=100000,STA=8000010;
 38 //const LL LNF=1LL<<60;
 39 const double EPS=1e-8;
 40 const double OO=1e15;
 41 const int dx[4]={-1,0,1,0};
 42 const int dy[4]={0,1,0,-1};
 43 const int day[13]={0,31,28,31,30,31,30,31,31,30,31,30,31};
 44 //Daily Use ...
 45 inline int sign(double x){return (x>EPS)-(x<-EPS);}
 46 template<class T> T gcd(T a,T b){return b?gcd(b,a%b):a;}
 47 template<class T> T lcm(T a,T b){return a/gcd(a,b)*b;}
 48 template<class T> inline T lcm(T a,T b,T d){return a/d*b;}
 49 template<class T> inline T Min(T a,T b){return a<b?a:b;}
 50 template<class T> inline T Max(T a,T b){return a>b?a:b;}
 51 template<class T> inline T Min(T a,T b,T c){return min(min(a, b),c);}
 52 template<class T> inline T Max(T a,T b,T c){return max(max(a, b),c);}
 53 template<class T> inline T Min(T a,T b,T c,T d){return min(min(a, b),min(c,d));}
 54 template<class T> inline T Max(T a,T b,T c,T d){return max(max(a, b),max(c,d));}
 55 //End
 56 
 57 struct Edge{
 58     int u,v,id;
 59 }e[N*N];
 60 int first[N],next[N*N];
 61 int level[N],ans[N][10000],cnt[N],vise[N*N];
 62 int Case,n,m,S,T,mt,tot;
 63 
 64 void adde(int a,int b,int c)
 65 {
 66     e[mt].u=a,e[mt].v=b,e[mt].id=c;
 67     next[mt]=first[a],first[a]=mt++;
 68     e[mt].u=b,e[mt].v=a,e[mt].id=c;
 69     next[mt]=first[b],first[b]=mt++;
 70 }
 71 
 72 int bfs(int flag)
 73 {
 74     int i,j,u,v,d;
 75     queue<int> q;
 76     q.push(S);
 77     level[S]=1;
 78     while(!q.empty()){
 79         u=q.front();q.pop();
 80         for(i=first[u];i!=-1;i=next[i]){
 81             if(flag){
 82                 if(!vise[e[i].id]){
 83                     vise[e[i].id]=1;
 84                     d=level[e[i].u]>tot?1:level[u];
 85                     ans[d][cnt[d]++]=e[i].id;
 86                 }
 87             }
 88         //    else if(e[i].v==T)return level[e[i].u];
 89             if(!level[e[i].v]){
 90                 level[e[i].v]=level[u]+1;
 91                 q.push(e[i].v);
 92             }
 93         }
 94     }
 95     return level[T]-1;
 96 }
 97 
 98 int main()
 99 {
100  //   freopen("in.txt","r",stdin);
101     int i,j,a,b;
102     scanf("%d",&Case);
103     while(Case--)
104     {
105         scanf("%d%d%d%d",&n,&m,&S,&T);
106         mem(first,-1);mt=0;
107         for(i=1;i<=m;i++){
108             scanf("%d%d",&a,&b);
109             adde(a,b,i);
110         }
111 
112         mem(level,0);
113         tot=bfs(0);
114         mem(level,0);mem(cnt,0);mem(vise,0);
115         bfs(1);
116 
117         printf("%d\n",tot);
118         for(i=1;i<=tot;i++){
119             printf("%d",cnt[i]);
120             sort(ans[i],ans[i]+cnt[i]);
121             for(j=0;j<cnt[i];j++){
122                 printf(" %d",ans[i][j]);
123             }
124             putchar('\n');
125         }
126     }
127     return 0;
128 }

 

posted @ 2013-07-23 00:58  zhsl  阅读(238)  评论(0编辑  收藏  举报