【连通分量】BZOJ 1093: [ZJOI2007]最大半连通子图
来源:互联网 发布:北京美工培训班学费 编辑:程序博客网 时间:2024/03/29 13:29
BZOJ 1093: [ZJOI2007]最大半连通子图
Description
Input
第一行包含两个整数N,M,X。N,M分别表示图G的点数与边数,X的意义如上文所述。接下来M行,每行两个正整数a, b,表示一条有向边(a, b)。图中的每个点将编号为1,2,3…N,保证输入中同一个(a,b)不会出现两次。
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
对于100%的数据, N ≤100000, M ≤1000000;对于100%的数据, X ≤10^8。
Solution
tarjan缩点
然后topo dp
维护一个dis 一个count
Code
#include<cstdio>#include<cstring>#include<algorithm>#include<cmath>#include<stack>#include<queue>using namespace std;#define maxn 1000001queue<int>q;stack<int>s;struct edge{ int to,lst;}e[maxn];int last[maxn],tot,n,m,low[maxn],dfn[maxn],Tot,num[maxn],size,f[maxn],mod,in[maxn];bool ins[maxn];void add(int u,int v){ e[++tot]=(edge){v,last[u]};last[u]=tot;}void tarjan(int poi){ dfn[poi]=low[poi]=++Tot; s.push(poi);ins[poi]=1; for(int i=last[poi];i;i=e[i].lst) { int u=e[i].to; if(!dfn[u]) { tarjan(u); low[poi]=min(low[poi],low[u]); } else if(ins[u])low[poi]=min(low[poi],low[u]); } if(dfn[poi]==low[poi]) { size++; int u; do{ u=s.top(); s.pop(); ins[u]=0; f[u]=size; num[size]++; }while(poi!=u); }}struct ED{ int u,v;}E[maxn];bool cmp(const ED A,const ED B){ return f[A.u]==f[B.u]?f[A.v]<f[B.v]:f[A.u]<f[B.u];}int top=0,F[maxn],cot[maxn];void topo_dp(){ for(int i=1;i<=size;i++)if(!in[i])q.push(i),F[i]=num[i],cot[i]=1; while(q.size()) { int u=q.front();q.pop(); for(int i=last[u];i;i=e[i].lst) { int v=e[i].to; in[v]--; if(!in[v])q.push(v); if(F[v]<F[u]+num[v])F[v]=F[u]+num[v],cot[v]=cot[u]; else if(F[v]==F[u]+num[v])cot[v]=(cot[v]+cot[u])%mod; } } int mx=0,sum=0; for(int i=1;i<=size;i++) { if(F[i]>F[mx])sum=cot[i]%mod,mx=i; else if(F[i]==F[mx])sum=(cot[i]+sum)%mod; } printf("%d\n%d",F[mx],sum);}int main(){ freopen("semi.in","r",stdin); freopen("semi.out","w",stdout); scanf("%d%d%d",&n,&m,&mod); for(int i=1;i<=m;i++) { int u,v; scanf("%d%d",&u,&v); add(u,v); E[i]=(ED){u,v}; } for(int i=1;i<=n;i++)if(!dfn[i])tarjan(i); sort(1+E,1+E+m,cmp); memset(last,0,sizeof(last));tot=0; for(int i=1;i<=m;i++) { int A=f[E[i].u],B=f[E[i].v]; if(A==f[E[i-1].u]&&B==f[E[i-1].v])continue; if(A==B)continue; add(A,B); in[B]++; } topo_dp(); return 0;}
0 0
- 【连通分量】BZOJ 1093: [ZJOI2007]最大半连通子图
- BZOJ 1093: [ZJOI2007]最大半连通子图
- BZOJ 1093: [ZJOI2007]最大半连通子图
- 【bzoj 1093】: [ZJOI2007]最大半连通子图
- 【BZOJ 1093】 [ZJOI2007]最大半连通子图
- BZOJ 1093: [ZJOI2007]最大半连通子图
- BZOJ 1093: [ZJOI2007]最大半连通子图
- BZOJ 1093: [ZJOI2007]最大半连通子图 强连通分量缩点,最长链,拓扑排序,DP
- 1093: [ZJOI2007]最大半连通子图
- 1093: [ZJOI2007]最大半连通子图
- 1093: [ZJOI2007]最大半连通子图
- 1093: [ZJOI2007]最大半连通子图
- 1093: [ZJOI2007]最大半连通子图
- [bzoj1093][ZJOI2007]最大半连通子图
- bzoj1093【ZJOI2007】最大半连通子图
- [BZOJ1093][ZJOI2007]最大半连通子图
- 【bzoj1093】 [ZJOI2007]最大半连通子图
- BZOJ1093: [ZJOI2007]最大半连通子图
- 检测远程端口是否打开
- 3351: [ioi2009]Regions
- 数据结构基础之栈与递归
- 数组方法归纳
- Win10关闭休眠方法
- 【连通分量】BZOJ 1093: [ZJOI2007]最大半连通子图
- 【Prufer数列/组合数学】[HNOI2008][HYSBZ/BZOJ1005]明明的烦恼
- Python开发简单爬虫(笔记)
- 为什么虚拟机里面的ubuntu不能共享windows的文件. /mnt/hgfs无效
- 记录下集中谷歌搜索方法
- 计算机网络3--计算机网络结构
- ubuntu 安装 paramon 集群监控前台
- MATLAB图像处理工具箱
- [Andoird自定义View]打分控件