巧妙的算法——差分约束

题目传送门，bzoj3436

这道几乎可以算是差分约束的板子题,那么什么是差分约束呢？

如果一个系统由n个变量和m个约束条件组成，形成m个形如ai-aj≤k的不等式，则称其为差分约束系统(system of difference constraints)。我们可以利用差分约束系统，构建成为一张有向图，从而用最短路（或最长路）解决问题。

---

例如我选的这道例题

输入为(1,a,b,c)时，显然此时满足的关系为a-b>=c,而a-b>=c等价于b-a<=-c，所以我们从a到b连一条为-c的边。

输入为(2,a,b,c)时，显然此时满足的关系为a-b<=c，所以我们从a到b连一条c的边。

输入为(3,a,b)时，显然此时满足的关系为a=b，则我们将a->b,b->a都连一条为0的边。

那么很显然，如果这张图中存在负环，就是矛盾的，输出No,否则不矛盾，输出Yes。

#include<bits/stdc++.h>
#include<vector>
using namespace std;
int i,j,k,n,m,tot,ans,r1,r2,r3,flag,dist[10005],vis[10005];
vector<pair<int,int> > f[10005];
int read(){
    char c;int x;while(c=getchar(),c<'0'||c>'9');x=c-'0';
    while(c=getchar(),c>='0'&&c<='9') x=x*10+c-'0';return x;
}
void tarjan(int x){
    int to,power;
    if(flag) return;
    for(int i=0;i<f[x].size();i++){
        to=f[x][i].first;power=f[x][i].second;
        if(dist[x]+power<dist[to]){
            if(vis[to]){
                flag=1;return;
            }
            else{
                vis[to]=1;
                dist[to]=dist[x]+power;
                tarjan(to);
            }
        }
    }
    vis[x]=0;
}
int main()
{
    n=read();m=read();flag=0;
    for(int i=1;i<=m;i++){
        int x=read();
        if(x==1) {r1=read(),r2=read(),r3=read();f[r2].push_back(make_pair(r1,-r3));}
        if(x==2) {r1=read(),r2=read(),r3=read(),f[r1].push_back(make_pair(r2,r3));}
        if(x==3) {r1=read(),r2=read();f[r1].push_back(make_pair(r2,0));f[r2].push_back(make_pair(r1,0));}
    }
    for(int i=1;i<=n;i++){memset(dist,0,sizeof(dist));tarjan(i);if(flag)break;}
    if(flag) puts("No");
    else puts("Yes");
    return 0;
}

注释：为什么存在负环即为矛盾，我们可以这样理解：i最多比j多f[i][j],如果图中没有矛盾，那一定满足u+f[u][u]>=u,但如果存在负环，则对于环上每一点u，都存在f[u][u]<0，u加上小于0的数肯定小于u，所以与u+f[u][u]>=u矛盾，所以这张图是不存在的。