bzoj2788 festival 差分约束

填坑中……链接：http://www.lydsy.com/JudgeOnline/problem.php?id=2788

题意：

有$n$个正整数$X1,X2,...,Xn$，再给出$m1+m2$个限制条件，限制分为两类：
1. 给出$a,b(1<=a,b<=n)$，要求满足$Xa + 1 = Xb$
2. 给出$c,d (1<=c,d<=n)$，要求满足$Xc <= Xd$
在满足所有限制的条件下，求集合${Xi}$大小的最大值。

首先看情况我们也知道是差分约束……

但是这个差分约束有些不一样在于这个东西会出环……所以我们要得出强连通分量……

得出来之后对于每个$SCC$分别考虑就行了……

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxn=605,maxm=200005;
const int inf=0x3f3f3f3f;
int n,m1,m2,cnt,tot,scnt;
int head[maxn],MAP[maxn][maxn];
int dfn[maxn],low[maxn],belong[maxn];
struct node
{
    int from,to,dis,next;
}edge[maxm];
void addedge(int u,int v,int w)
{
    edge[++tot]=(node){u,v,w,head[u]};head[u]=tot;
}
#include<stack>
stack<int>s;
void tarjan(int root)
{
    dfn[root]=low[root]=++cnt;s.push(root);
    for(int i=head[root];i;i=edge[i].next)
    {
        int v=edge[i].to;
        if(!dfn[v])
        {
            tarjan(v);
            low[root]=min(low[root],low[v]);
        }
        else if(!belong[v])low[root]=min(low[root],dfn[v]);
    }
    if(low[root]==dfn[root])
    {
        scnt++;int k;
        do
        {
            k=s.top();s.pop();
            belong[k]=scnt;
        }while(k!=root);
    }
}
int haha()
{
    scanf("%d%d%d",&n,&m1,&m2);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)MAP[i][j]=-inf;
    for(int i=1;i<=n;i++)MAP[i][i]=0;
    for(int i=1;i<=m1;i++)
    {
        int x,y;scanf("%d%d",&x,&y);
        addedge(x,y,1),addedge(y,x,-1);
        MAP[x][y]=max(MAP[x][y],1);MAP[y][x]=max(MAP[y][x],-1);
    }
    for(int i=1;i<=m2;i++)
    {
        int x,y;scanf("%d%d",&x,&y);
        addedge(x,y,0);MAP[x][y]=max(MAP[x][y],0);
    }
    for(int i=1;i<=n;i++)
        if(!dfn[i])tarjan(i);
    int ans=0;
    for(int b=1;b<=scnt;b++)
    {
        int ret=0;
        for(int k=1;k<=n;k++)
        {
            if(belong[k]!=b)continue;
            for(int i=1;i<=n;i++)
            {
                if(belong[i]!=b)continue;
                if(MAP[i][k]==-inf)continue;
                for(int j=1;j<=n;j++)
                {
                    if(belong[j]!=b)continue;
                    if(MAP[k][j]==-inf)continue;
                    MAP[i][j]=max(MAP[i][j],MAP[i][k]+MAP[k][j]);
                }
            }
        }
        for(int i=1;i<=n;i++)
        {
            if(belong[i]!=b)continue;
            for(int j=1;j<=n;j++)
            {
                if(belong[j]!=b)continue;
                ret=max(ret,abs(MAP[i][j]));
            }
        }
        ans+=ret+1;
    }
    for(int i=1;i<=n;i++)
        if(MAP[i][i])
        {
            puts("NIE");
            return 0;
        }
    printf("%d
",ans);
}
int sb=haha();
int main(){;}