差分约束—学习笔记

差分约束系统（system of difference constraints）

如果一个系统由n个变量和m个约束条件组成
其中每个约束条件形如xj-xi<=bk(i,j∈[1,n],k∈[1,m])
则称其为差分约束系统。
亦即，差分约束系统是求解关于一组变量的特殊不等式组的方法。

通俗一点地说，差分约束系统就是给定一些不等式的组
而我们的目标是通过给定的约束条件
求出最大值或者最小值或者差分约束系统是否有解

---

差分约束求解

差分约束系统可以通过图论求解

对于型如$x_v-x_u<=k$的不等式
我们从u到v连一条权值为k的边

若约束条件型如$x_v-x_u=k$
那么可以转化为$x_v-x_u<=k$且$x_v-x_u>=k\Rightarrow x_u-x_v<=-k$
可以连接两条边权相反的边
或者可以连接两条边权为0的边

建图后若图中存在负环，则无解
若要求差分约束系统的最小解，则在图中跑最长路
反之则跑最短路

特别要注意的图有可能不连通
这时候我们只需要连接一个超级源点即可

---

应用

洛谷P1993 小K的农场

若k==1，则有$x_a-x_b>=c\Rightarrow x_b-x_a<=-c$，从a向b连接边权为-c的边

若k==2，则有$x_a-x_b<=c$，从b向a连接边权为c的边

若k==2，则有$x_a==x_b$，从a、b间连接边权为0的无向边
最后判负环即可

#include<iostream>
#include<vector>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cstdio>
using namespace std;

int read()
{
    int f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
    while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
    return f*x;
}

void print(int x)
{
    if(x<0){putchar('-');x=-x;}
    if(x>9)print(x/10);
    putchar(x%10+'0');
}

int n,m;
int tot;
int s=0;
struct node{int v,dis,nxt;}E[1000010];
int head[1000010];
int d[1000010];
bool vis[1000010];
bool judge;

void add(int u,int v,int dis)	
{
    E[++tot].nxt=head[u];
    E[tot].dis=dis;
    E[tot].v=v;
    head[u]=tot;
}

void spfa(int u)
{
    vis[u]=true;
    for(int i=head[u];i;i=E[i].nxt)
    {
        int v=E[i].v;
        if(d[v]>d[u]+E[i].dis)
        {
            if(vis[v]||judge){judge=true;break;}
            d[v]=d[u]+E[i].dis;
            spfa(v);
        }
    }
    vis[u]=false;
}

int main()
{
    n=read();m=read();
    for(int i=1;i<=m;i++)
    {
        int k=read();
        int a=read(),b=read();
        if(k==1){ int c=read(); add(a,b,-c); }
        else if(k==2){ int c=read(); add(b,a,c); }
        else if(k==3){ add(a,b,0); add(b,a,0); }
    }
    for(int i=1;i<=n;i++)
    {
        spfa(i);
        if(judge)break;
    }
    
    if(!judge)cout<<"Yes";
    else if(judge)cout<<"No";
    return 0;
    
}

---

洛谷P3275 [SCOI2011]糖果【差分约束】

为了求出差分约束系统每个x的最小值
我们需要在建好的图上跑最长路

同时因为图有可能不连通，所里需要建立一个超级源点
最后答案就是超级源点到每个点最短路长度的和

#include<iostream>
#include<cstdio>
#include<vector>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;
typedef long long ll;

ll read()
{
    ll f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
    while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
    return f*x;
}

void print(ll x)
{
    if(x<0){putchar('-');x=-x;}
    if(x>9)print(x/10);
    putchar(x%10+'0');
}

int n,k;
ll tot;
struct node{ll v,dis,nxt;}E[1000010];
ll head[1000010];
ll d[1000010];
ll vis[1000010];
ll cnt[1000010];
ll ans;

void add(ll u,ll v,ll dis)
{
    E[++tot].nxt=head[u];
    E[tot].dis=dis;
    E[tot].v=v;
    head[u]=tot; 
}

void spfa()
{
    queue<ll> q; q.push(0);
    for(int i=1;i<=n;i++)d[i]=-1e9; d[0]=0;
    cnt[0]=1;
    
    while(!q.empty())
    {
        ll u=q.front();
        q.pop(); vis[u]=false;
        for(ll i=head[u];i;i=E[i].nxt)
        {
            ll v=E[i].v; 
            if(d[v]<d[u]+E[i].dis)
            {
                d[v]=d[u]+E[i].dis;
                if(!vis[v])
                {
                    cnt[v]++;
                    if(cnt[v]>=n) {print(-1);exit(0);}
                    q.push(v); vis[v]=true;
                } 
            }
        } 
    }
}

int main()
{
    n=read();k=read();
    while(k--)
    {
        ll x=read();
        ll a=read(),b=read();
        if(x==1)
        add(a,b,0),add(b,a,0);
        else if(x==2)
        {
            if(a==b) {print(-1);return 0;}
            else add(a,b,1);
        } 
        else if(x==3) add(b,a,0);
        else if(x==4)
        {
            if(a==b) {print(-1);return 0;}
            else add(b,a,1);
        }
        else if(x==5) add(a,b,0);
    }
    for(int i=n;i>=1;i--) add(0,i,1);
    
    spfa();
    for(int i=1;i<=n;i++)ans+=d[i];
    print(ans);
    return 0;
}

---

洛谷P2294 [HNOI2005]狡猾的商人【差分约束】

没什么好讲的，加个前缀和就可以差分约束套路建图了