BZOJ 1093 [ZJOI2007]最大半连通子图

以前做过poj的一个判断图是否为弱连通的题，然后，这个题和poj那个差不多。

先强连通缩点，然后重新构图，然后找出包含点数最多的链，统计个数即可，可以用拓扑排序搞~

 

pS：重新构图时有重边，然后导致统计方案数的重复。。wa了好久。。还是wzc神犇告诉我这个蒟蒻的。。

 

View Code

  1 #include <iostream>

  2 #include <cstdio>

  3 #include <cstdlib>

  4 #include <cstring>

  5 #include <algorithm>

  6  

  7 #define N 200000

  8 #define M 5000000

  9 #define BUG system("pause")

 10  

 11 using namespace std;

 12  

 13 int head[N],to[M],next[M];

 14 int dfn[N],low[N];

 15 int st[M],ed[M];

 16 int n,m,cnt,ans,ansnum,mod;

 17 int divg,belong[N],t,p,stk[N],val[N];

 18 bool fg[N];

 19 int num[N],in[N],dp[N],q[M],vis[N];

 20  

 21 inline void add(int u,int v)

 22 {

 23     to[cnt]=v; next[cnt]=head[u]; head[u]=cnt++;

 24 }

 25  

 26 inline void read()

 27 {

 28     memset(head,-1,sizeof head); cnt=0;

 29     scanf("%d%d%d",&n,&m,&mod);

 30     for(int i=1;i<=m;i++)

 31     {

 32         scanf("%d%d",&st[i],&ed[i]);

 33         add(st[i],ed[i]);

 34     }

 35 }

 36  

 37 inline void dfs(int u)

 38 {

 39     low[u]=dfn[u]=++t;

 40     stk[++p]=u; fg[u]=true;

 41     for(int i=head[u];~i;i=next[i])

 42     {

 43         if(!dfn[to[i]])

 44         {

 45             dfs(to[i]);

 46             low[u]=min(low[u],low[to[i]]);

 47         }

 48         else if(fg[to[i]]) low[u]=min(low[u],dfn[to[i]]);

 49     }

 50     if(dfn[u]==low[u])

 51     {

 52         divg++;

 53         int tmp=-1;

 54         while(tmp!=u)

 55         {

 56             tmp=stk[p--];

 57             belong[tmp]=divg;

 58             val[divg]++;

 59             fg[tmp]=false;

 60         }

 61     }

 62 }

 63  

 64 inline void topsort()

 65 {

 66     int h=1,t=1,u;

 67     for(int i=1;i<=divg;i++)

 68         if(in[i]==0)

 69         {

 70             q[t++]=i;

 71             dp[i]=val[i];

 72             num[i]=1;

 73         }

 74     while(h<t)

 75     {

 76         u=q[h++];

 77         for(int i=head[u];~i;i=next[i])

 78         {

 79             in[to[i]]--;

 80             if(in[to[i]]==0) q[t++]=to[i];

 81             if(vis[to[i]]==u) continue;//有重边!! 

 82             if(dp[to[i]]<dp[u]+val[to[i]])

 83             {

 84                 dp[to[i]]=dp[u]+val[to[i]];

 85                 num[to[i]]=num[u];

 86             }

 87             else if(dp[to[i]]==dp[u]+val[to[i]])

 88             {

 89                 num[to[i]]=(num[to[i]]+num[u])%mod;

 90             }

 91             vis[to[i]]=u;

 92         }

 93     }

 94 }

 95  

 96 inline void go()

 97 {

 98     for(int i=1;i<=n;i++)

 99         if(!dfn[i]) dfs(i);

100     memset(head,-1,sizeof head); cnt=0;

101     for(int i=1;i<=m;i++)

102         if(belong[st[i]]!=belong[ed[i]])

103         {

104             add(belong[st[i]],belong[ed[i]]);

105             in[belong[ed[i]]]++;

106         }

107     topsort();

108     for(int i=1;i<=divg;i++)

109     {

110         if(dp[i]>ans)

111         {

112             ans=dp[i];

113             ansnum=num[i];

114         }

115         else if(dp[i]==ans) ansnum=(ansnum+num[i])%mod;

116     }

117     printf("%d\n%d\n",ans,ansnum);

118 }

119  

120 int main()

121 {

122     read();

123     go();

124     return 0;

125 }