BZOJ2330: [SCOI2011]糖果

发现自己学了假的差分约束qaq
有几种关系
$x=y/x>y/x>=y$
分别对应$x>=y,x<=y/x-1>=y/x>=y$

所有关系都是形如$x+c>=y$的形式，要求总数最小，可以把关系变成$x>=y-c$的形式，把x，y看成到x，y的最长路，发现就是一条y连向x，权为$-c$的边，此时因为对于每个x，取得值是满足所有不等式$x>=di$的下界，和是最小的

此题直接用spfa判正环会T，注意到原图只有0/1两类边，可以先建0边，tarjan缩点，在缩点后的图建1边，此时如果有环就说明有正环

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

inline void up(int &x,const int &y){if(x<y)x=y;}
inline void down(int &x,const int &y){if(x>y)x=y;}
const int maxn = 110000;
const int maxm = 210000;

int n,m;
int e[maxm][3],tot;
struct edge{int y,nex;}a[maxm]; int len,fir[maxn];
inline void ins(const int x,const int y){a[++len]=(edge){y,fir[x]};fir[x]=len;}
vector<int>V[maxn],vc[maxn];

int dfn[maxn],low[maxn],id;
int t[maxn],tp; bool insta[maxn];
int bel[maxn],bn[maxn],cnt;
void tarjan(const int x)
{
    low[x]=dfn[x]=++id; insta[t[++tp]=x]=true;
    for(int k=fir[x],y=a[k].y;k;k=a[k].nex,y=a[k].y)
    {
        if(!dfn[y]) tarjan(y),down(low[x],low[y]);
        else if(insta[y]) down(low[x],dfn[y]);
    }
    if(low[x]==dfn[x])
    {
        ++cnt; int la=0;
        while(la!=x)
        {
            insta[la=t[tp--]]=false;
            bn[bel[la]=cnt]++;
        }
    }
}
int ind[maxn];
int f[maxn];
queue<int>q;

int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        int op,x,y; scanf("%d%d%d",&op,&x,&y);
        if(op==1) 
        {
            e[++tot][0]=x,e[tot][1]=y; ins(x,y);
            e[++tot][0]=y,e[tot][1]=x; ins(y,x);
        }
        else if(op==2) e[++tot][0]=x,e[tot][1]=y,e[tot][2]=1;
        else if(op==3) e[++tot][0]=y,e[tot][1]=x,ins(y,x);
        else if(op==4) e[++tot][0]=y,e[tot][1]=x,e[tot][2]=1;
        else e[++tot][0]=x,e[tot][1]=y,ins(x,y);
    }
    for(int i=1;i<=n;i++) if(!dfn[i]) tarjan(i);

    for(int i=1;i<=tot;i++) 
    {
        if(e[i][2])
        {
            int x=bel[e[i][0]],y=bel[e[i][1]];
            if(x==y) return puts("-1"),0;
            V[x].push_back(y); vc[x].push_back(1); ind[y]++;
        }
        else
        {
            int x=bel[e[i][0]],y=bel[e[i][1]];
            if(x!=y) V[x].push_back(y),vc[x].push_back(0),ind[y]++;
        }
    }
    for(int i=1;i<=cnt;i++) if(!ind[i]) f[i]=1,q.push(i);
    while(!q.empty())
    {
        const int x=q.front(); q.pop();
        for(int i=0;i<V[x].size();i++)
        {
            int y=V[x][i];
            up(f[y],f[x]+vc[x][i]);
            ind[y]--; if(!ind[y]) q.push(y);
        }
    }
    ll re=0;
    for(int i=1;i<=cnt;i++)
    {
        if(ind[i]) return puts("-1"),0;
        re+=(ll)f[i]*bn[i];
    }
    printf("%lld\n",re);

    return 0;
}