首页 > 代码库 > bzoj1093: [ZJOI2007]最大半连通子图

bzoj1093: [ZJOI2007]最大半连通子图

惨烈啊。。。int son[x]=>bool son[x]一直调不出来我也是醉了。!!!最新错法。。。 
缩点后有重边!!! 就是缩点之后找最长路然后找有多少条最长路树形dp一下。
#include<cstdio> #include<cstring> #include<cctype> #include<algorithm> #include<stack> using namespace std; #define rep(i,s,t) for(int i=s;i<=t;i++) #define dwn(i,s,t) for(int i=s;i>=t;i--) #define clr(x,c) memset(x,c,sizeof(x)) #define qwq(x) for(edge *o=head[x];o;o=o->next) #define qaq(x) for(edge *o=h[x];o;o=o->next) int read(){     int x=0,f=1;char c=getchar();     while(!isdigit(c)) {        if(c==‘-‘) f=-1;c=getchar();    }    while(isdigit(c)) x=x*10+c-‘0‘,c=getchar();     return x*f; } const int nmax=1e5+5; const int maxn=1e6+5; const int inf=0x7f7f7f7f; struct edge{     int to;edge *next; }; edge es[maxn<<1],*pt=es,*head[nmax]; void add(int u,int v){     pt->to=v;pt->next=head[u];head[u]=pt++; } int pre[nmax],scc[nmax],size[nmax],dfs_clock,scc_cnt,mod; stack<int>s; void mins(int &a,int b){     if(a>b) a=b; } void maxs(int &a,int b){     if(a<b) a=b; } int dfs(int x){     int lowu=pre[x]=++dfs_clock;s.push(x);     qwq(x) {         if(!pre[o->to]) mins(lowu,dfs(o->to));         else if(!scc[o->to]) mins(lowu,pre[o->to]);     }     if(lowu==pre[x]){         scc_cnt++;int o,cnt=0;         while(1){             o=s.top();s.pop();             scc[o]=scc_cnt;cnt++;             if(o==x) break;         }         size[scc_cnt]=cnt;     }     return lowu; }    edge *h[nmax]; bool vis[nmax];int son[nmax],sum[nmax],orz[nmax]; void adde(int u,int v){     pt->to=v;pt->next=h[u];h[u]=pt++; } void Dfs(int x){     vis[x]=1;int cnt=0;     qaq(x) {         if(!vis[o->to]) Dfs(o->to);maxs(cnt,son[o->to]);     }     son[x]=cnt+size[x]; } void DFS(int x){     vis[x]=1;    qaq(x){         if(!vis[o->to]) DFS(o->to);         if(son[x]==son[o->to]+size[x]&&orz[o->to]!=x) sum[x]=(sum[x]+sum[o->to])%mod,orz[o->to]=x;     }     if(!h[x]) sum[x]=1; } int main(){     int n=read(),m=read(),u,v;mod=read();     rep(i,1,m) u=read(),v=read(),add(u,v);     rep(i,1,n) if(!pre[i]) dfs(i);     //rep(i,1,n) printf("%d ",scc[i]);printf("\n");            rep(i,1,n) {         u=scc[i];         qwq(i) if(u!=scc[o->to]) adde(u,scc[o->to]);     }     int ans=0;     rep(i,1,scc_cnt) if(!vis[i]) Dfs(i),maxs(ans,son[i]);     //rep(i,1,scc_cnt) printf("%d ",son[i]);printf("\n");     int res=0;clr(vis,0);     rep(i,1,scc_cnt) if(!vis[i]) {         if(son[i]==ans) DFS(i),res=(res+sum[i])%mod;     }     printf("%d\n%d\n",ans,res);     return 0; } 

  

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

Time Limit: 30 Sec  Memory Limit: 162 MB
Submit: 2725  Solved: 1079
[Submit][Status][Discuss]

Description

  一个有向图G=(V,E)称为半连通的(Semi-Connected),如果满足:?u,v∈V,满足u→v或v→u,即对于图中任意
两点u,v,存在一条u到v的有向路径或者从v到u的有向路径。若G‘=(V‘,E‘)满足V‘?V,E‘是E中所有跟V‘有关的边,
则称G‘是G的一个导出子图。若G‘是G的导出子图,且G‘半连通,则称G‘为G的半连通子图。若G‘是G所有半连通子图
中包含节点数最多的,则称G‘是G的最大半连通子图。给定一个有向图G,请求出G的最大半连通子图拥有的节点数K
,以及不同的最大半连通子图的数目C。由于C可能比较大,仅要求输出C对X的余数。

Input

  第一行包含两个整数N,M,X。N,M分别表示图G的点数与边数,X的意义如上文所述接下来M行,每行两个正整
数a, b,表示一条有向边(a, b)。图中的每个点将编号为1,2,3…N,保证输入中同一个(a,b)不会出现两次。N ≤1
00000, M ≤1000000;对于100%的数据, X ≤10^8

Output

  应包含两行,第一行包含一个整数K。第二行包含整数C Mod X.

Sample Input

6 6 20070603
1 2
2 1
1 3
2 4
5 6
6 4

Sample Output

3
3

HINT

 

Source

 
[Submit][Status][Discuss]
 

bzoj1093: [ZJOI2007]最大半连通子图