浅谈差分约束系统——图论不等式的变形

Posted YJSheep

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了浅谈差分约束系统——图论不等式的变形相关的知识,希望对你有一定的参考价值。

浅谈差分约束系统——图论不等式的变形

                                                                 ----yangyaojia

 

版权声明:本篇随笔版权归作者YJSheep(www.cnblogs.com/yangyaojia)所有,转载请保留原地址!

一.定义

如若一个系统由n个变量和m个不等式组成,并且这m个不等式对应的系数矩阵中每一行有且仅有一个1和-1,其它的都为0,这样的系统称为差分约束( difference constraints )系统。

二.分析

简单来说就是给你n个变量,给m个形如x[i]-x[j]≥k①或x[i]-x[j]≤k②。求两个变量差的最大与最小值。我们可以对不等式稍稍变形,将①变为k+aj≤ai,②变为ai+(-k)≤aj。不难发现,这和我们的最短路问题中的不等式dis[i]+w[i,j]≤dis[j]极为相似。

最大值

那么,如果题目要求我们求出某两点之间的最大距离,我们可以这样:

  1. 对于x[i] - x[j]≤a[k],我们可以建一条由j到i的有向边,权值为k
  2. 对于x[i] - x[j]≥a[k],我们可以建一条由i到j的有向边,权值为-k

并进行如下判断。

  1. 如果存在最短路径可以到达这两点,答案就是最短路的长度。
  2. 如果不存在路径可以到达,那么两点间的距离可以是无限大。
  3. 如果图中存在负环,那么就无解。

最小值

如果题目要求求出两点间的最小距离,我们就建立形如x[i]-x[j]≥k,建立i到j,权值为k的路,求的最长路径。判断也跟上面一样,不过是判有没有正环。

 

下面来附上一道典型例题 POJ3169

#include<queue>
#include <algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#define Pair pair<int,int>
#define MAXN 1000+1
#define MAX 99999999
#define MAXM 40000+1
using namespace std;
int n,m1,m2,num,nnum,head[MAXN],s,o,pre[MAXN],v[MAXN],K;
int dis[MAXN],t[MAXN],ans=MAX;
struct Edge{
    int next,to,exi,from;
    int dis,cost;
}edge[MAXM];
int read()
{
    int in=0;char c;
    c=getchar();
    for(;c>\'9\'||c<\'0\';c=getchar());
    for(;c<=\'9\'&&c>=\'0\';c=getchar()) in=in*10+c-\'0\';
    return in;
}
void add(int from,int to,int dis)
{
    edge[++num].next=head[from];
    edge[num].to=to;
    edge[num].dis=dis;
    edge[num].from=from;
    head[from]=num;
}
int spfa()
{
    for(int i=1;i<=n;i++)
     dis[i]=99999999;
    queue<int> q;
    dis[s]=0;t[s]++;
    q.push(s);
    while(q.size()>0)
    {
        int k=q.front();q.pop();
        v[k]=0;
        for(int i=head[k];i;i=edge[i].next)
        {
            if(dis[edge[i].to]>dis[k]+edge[i].dis)
            {
                dis[edge[i].to]=dis[k]+edge[i].dis;
                t[edge[i].to]++;
                if(t[edge[i].to]>n) return 1;
                if(!v[edge[i].to]) v[edge[i].to]=1,q.push(edge[i].to);
                
            }
        }
    }
    return 0;
} 

int main()
{
    scanf("%d%d%d",&n,&m1,&m2);
    int a,b,c;
    for(int i=1;i<=m1;i++)
    {
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c); 
    }
    for(int i=1;i<=m2;i++)
    {
        scanf("%d%d%d",&a,&b,&c);
        add(b,a,-c);
    }s=1;
    if(spfa()) 
    {
        printf("-1\\n");
        return 0;
    }
    if(dis[n]==99999999) printf("-2\\n");
    else printf("%d\\n",dis[n]);
    return 0;
}

题目大意就是个一些约束条件求最大值,只要按着方法去建模,判断,就可以出解。

三.其他

         差分约束还分线性约束,区间约束与未知条件的约束。但本质上就是一样的,都是建立最短路模型来求解。关于未知条件的约束的题目,也是枚举k的大小求解。在这里也不再多说,毕竟是浅谈差分约束。

 

以上是关于浅谈差分约束系统——图论不等式的变形的主要内容,如果未能解决你的问题,请参考以下文章

浅谈差分约束

浅谈差分约束系统

浅谈差分约束系统

关于图论差分约束

差分约束系统

初涉差分约束系统