博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
[Poi2012]Festival 差分约束+tarjan
阅读量:7226 次
发布时间:2019-06-29

本文共 1768 字,大约阅读时间需要 5 分钟。

差分约束建图,发现要在每个联通块里求最长路,600,直接O(n3) floyed

#include
#include
#include
#include
#include
#define N 650#define M 100050using namespace std;int g[N][N],n,m1,m2,f[N],ans;int e=1,head[N];struct edge{ int u,v,w,next;}ed[2*M];void add(int u,int v,int w){ ed[e].u=u;ed[e].v=v;ed[e].w=w; ed[e].next=head[u];head[u]=e++; g[u][v]=min(g[u][v],w);}int dfn[N],low[N],top,q[N],tot,id[N];bool bo[N];void tarjan(int x){ dfn[x]=low[x]=++top; q[top]=x;bo[x]=1; for(int i=head[x];i;i=ed[i].next){ int v=ed[i].v; if(!dfn[v]){ tarjan(v); low[x]=min(low[x],low[v]); } else if(bo[v]) low[x]=min(low[x],dfn[v]); } if(low[x]==dfn[x]){ int y;tot++; do{ y=q[top--]; bo[y]=0; id[y]=tot; }while(y!=x); }}int main(){ memset(g,0x3f,sizeof g); scanf("%d%d%d",&n,&m1,&m2); for(int i=1;i<=n;i++)g[i][i]=0; int u,v; for(int i=1;i<=m1;i++){ scanf("%d%d",&u,&v); add(u,v,1);add(v,u,-1); } for(int i=1;i<=m2;i++){ scanf("%d%d",&u,&v); add(v,u,0); } for(int i=1;i<=n;i++) if(!dfn[i]) tarjan(i); for(int k=1;k<=n;k++) for(int i=1;i<=n;i++)if(id[i]==id[k]) for(int j=1;j<=n;j++)if(id[j]==id[i]) g[i][j]=min(g[i][j],g[i][k]+g[k][j]); for(int i=1;i<=n;i++) if(g[i][i]<0){printf("NIE\n");return 0;} for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) if(id[i]==id[j]) f[id[i]]=max(f[id[i]],g[i][j]+1); for(int i=1;i<=tot;i++) ans+=f[i]; printf("%d\n",ans); return 0;}

转载于:https://www.cnblogs.com/Ren-Ivan/p/7746676.html

你可能感兴趣的文章
Bytescout Spreadsheet SDK for.NET
查看>>
我的友情链接
查看>>
Haproxy的三种保持客户端会话保持方式
查看>>
iOS的数学函数
查看>>
python 模块 chardet下载及介绍(转)
查看>>
能力工场--关于在JavaScript中使用EL表达式的问题
查看>>
NFS服务器设置
查看>>
s:iterator 中的status 使用方法
查看>>
cocos2d-x 源码剖析系列
查看>>
IT系统架构设计
查看>>
Nginx虚拟主机配置实践(一)
查看>>
细谈Spring(一)spring简介
查看>>
网络工程师的面试题
查看>>
nginx启动脚本
查看>>
常用输入法框架简介
查看>>
记录新机房建设。20130629
查看>>
安装ntop
查看>>
ssh远程登录讲解
查看>>
mysql的备份脚本
查看>>
linux下mysql的root密码忘记解决方法
查看>>