bzoj 2788(差分约束+tarjan)

来源:互联网 发布:如何卸载mac第三方软件 编辑:程序博客网 时间:2024/06/11 23:25

传送门
题解:
限制1连边(a,b,1),(b,a,-1),限制2连边(d,c,0)。
tarjan缩点后每个(强连通分量内的最短路最大值+1)之和就是答案(集合大小是指不同的元素个数)。
P.S.status里ac的分为两派,10000+ms和1000+ms,经过和同班神犇的对比,并没有在算法和实现方法上发现任何区别。迷之自带大常数orz。。。

#include<bits/stdc++.h>using namespace std;const int MAXN=605,MAXM=1e5+5;int n,m1,m2;int head[MAXN],edge=0;struct EDGE {    int v,nxt;}e[MAXM<<1];int low[MAXN],dfn[MAXN],tim=0,snt=0;bool ins[MAXN];stack<int> S;int blo[MAXN][MAXN];int f[MAXN][MAXN],ans=0;inline int cal(int p) {    int ret=0;    for (int k=1;k<=blo[p][0];++k)        for (int i=1;i<=blo[p][0];++i)            for (int j=1;j<=blo[p][0];++j)                f[blo[p][i]][blo[p][j]]=min(f[blo[p][i]][blo[p][j]],f[blo[p][i]][blo[p][k]]+f[blo[p][k]][blo[p][j]]);    for (int i=1;i<=blo[p][0];++i)        for (int j=1;j<=blo[p][0];++j)            ret=max(ret,f[blo[p][i]][blo[p][j]]);    return ret;}void adde(int u,int v) {    e[++edge].nxt=head[u],e[edge].v=v,head[u]=edge;}void dfs(int p) {    low[p]=dfn[p]=++tim;    ins[p]=true;    S.push(p);    for (int i=head[p];~i;i=e[i].nxt) {        int v=e[i].v;        if (!dfn[v]) {            dfs(v);            low[p]=min(low[p],low[v]);        }        else if (ins[v]) low[p]=min(low[p],dfn[v]);    }    if (low[p]==dfn[p]) {        ++snt;        while (1) {            int t=S.top();            S.pop();            ins[t]=false;            blo[snt][++blo[snt][0]]=t;            if (t==p) break;        }        ans+=cal(snt)+1;    }}int main() {//  freopen("bzoj 2788.in","r",stdin);    memset(f,0x3f,sizeof(f));    memset(head,-1,sizeof(head));    scanf("%d%d%d",&n,&m1,&m2);    for (int i=1;i<=n;++i) f[i][i]=0;    for (register int i=0;i<m1;++i) {        int a,b;        scanf("%d%d",&a,&b);        adde(a,b),adde(b,a);        f[a][b]=min(f[a][b],1),f[b][a]=min(f[b][a],-1);    }    for (register int i=0;i<m2;++i) {        int c,d;        scanf("%d%d",&c,&d);        adde(d,c);        f[d][c]=min(f[d][c],0);    }    for (int i=1;i<=n;++i)        if (!dfn[i]) dfs(i);    for (int i=1;i<=n;++i)        if (f[i][i]<0) return puts("NIE"),0;    printf("%d\n",ans);    return 0;}
原创粉丝点击