bzoj 2788(差分约束+tarjan)

传送门
题解:
限制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;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值