差分约束系统

最新推荐文章于 2023-09-16 17:16:43 发布
最新推荐文章于 2023-09-16 17:16:43 发布
阅读量506 收藏
点赞数
文章标签： 差分约束系统 poj dijkstra
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u013475704/article/details/44024001
版权

         首先说说定义，差分约束系统是指一组不等式，且每个不等式都由两个变量和一个常数组成，形如v<=u+w或v>=u+w，u和v是变量，w是常数。每个不等式都是一个约束条件，约束目的是使目标函数xs-xt最大或最小。

        差分约束系统属于线性规划问题，不过神奇的是可以转换为图论问题，先做些简单的说明，就像有一组不等式：

         b-a<=k1,

         c-b<=k2,

         c-a<=k3,

         上述不等式组就是一个差分约束系统，每个不等式都是两个变量之差小于（当然也可以大于）等于某个常数。现在要求c-a的最大值，明显，这里可以用线性规划解决，但我们发现，上述每个不等式都与求解最短路径问题时的三角不等式相像，即边u->v 的权值为w的话，有d[v]<=d[u]+w，变换下就有d[v]-d[u]<=w，因此，我们可以将a,b,c三个变量看做三个节点，k1,k2,k3看做三个边的权值，于是可以构建出下图：

                                                                

         求c-a的最大值也就是求从a到c的最短路，事实上，由b-a<=k1和c-b<=k2可得c-a<=k1+k2，则k1+k2和k3的最小值就是c-a的最大值，从这里我们也可以从直观上理解为什么从a到c的最短路就是c-a的最大值。

        下面再以另一个例子做进一步介绍。

      X1 - X2 <= 0
      X1 - X5 <= -1
      X2 - X5 <= 1
      X3 - X1 <= 5
      X4 - X1 <= 4
      X4 - X3 <= -1
      X5 - X3 <= -3
      X5 - X4 <= -3

         这也是一个差分约束系统，总共有5个变量，8个约束条件。要求满足所有不等式的一组解。

        由前一个例子知道，差分约束系统的解法就是利用了解单源最短路径时的三角不等式d[v]-d[u]<=w，对这类问题，在找出所有不等式后，就是要建图。5个变量就是5个节点，8个约束条件就是8条边，每个约束条件xi-xj<=k代表一条有向边<xj,xi>的权值为k(如果是xi-xj>=k，变为xj-xi<=-k)，求出新建图的单源最短路径就是一组解。

       另外，差分约束系统要么无解，要么有无数组解，因为当在所有不等式两边同时加上某个数时，就可以得到另外一组解。比如在本例中，{-5，-3，0，-1，-4}是一组解，加5后，{0,2,5,4,1}是另一组解。

      当然，既然是求单源最短路径，自然要有一个源点，而这个源点需要我们自己加上去，否则图可能不连通。通常这个源点与其它节点的距离设为0（当然也可以设为任何比INF小的值，加入源点只是为了要使图保持连通性，同时也更好理解一些）。假设加入的源点为x0，从不等式的角度出发，也可以理解为xi-x0<=0，i=1,2,3,4,5。于是本例构建的图如下：

                                                

       接下来只需要求解从源点v0出发，到达其它各点的最短路就行了。

       另外需要注意的几点：

      1：要求目标函数vt-vs的最大值，需要把所有的不等式化为vi-vj<=k的形式:，建图后求最短路。求最小值的话就化为vi-vj>=k的形式，建图后求最长路。

      2：若图中存在负权回路，这差分约束系统无解。

      3：若只需要判断差分约束系统是否有解，那么可以不用加入源点，用bellman-ford算法判断是否有负权回路，但要注意的是这时不能用spfa算法，因为前者不管图是否连通，都会对所有边进行松弛操作，而后者就不行了。

       下面通过poj上的一些题来作为练习。

       poj 3159

       很典型的差分约束系统的题。输入表示有n个点，m个约束条件，每个约束条件给出为 A B c，表示d[B]-d[A]<=c，要求d[n]-d[1]的最大值。直接建图求最短路就行了。

      代码如下：

#include <cstdio>
#include <stack>
#include <set>
#include <iostream>
#include <string>
#include <vector>
#include <queue>
#include <functional>
#include <cstring>
#include <algorithm>
#include <cctype>
#include <string>
#include <map>
#include <iomanip>
#include <cmath>
#define LL long long
#define ULL unsigned long long
#define SZ(x) (int)x.size()
#define Lowbit(x) ((x) & (-x))
#define MP(a, b) make_pair(a, b)
#define MS(arr, num) memset(arr, num, sizeof(arr))
#define PB push_back
#define F first
#define S second
#define ROP freopen("input.txt", "r", stdin);
#define MID(a, b) (a + ((b - a) >> 1))
#define LC rt << 1, l, mid
#define RC rt << 1|1, mid + 1, r
#define LRT rt << 1
#define RRT rt << 1|1
#define BitCount(x) __builtin_popcount(x)
#define BitCountll(x) __builtin_popcountll(x)
#define LeftPos(x) 32 - __builtin_clz(x) - 1
#define LeftPosll(x) 64 - __builtin_clzll(x) - 1
const double PI = acos(-1.0);
const int INF = 0x3f3f3f3f;
using namespace std;
const double eps = 1e-8;
const int MAXN = 300 + 10;
const int MOD = 1000007;
const int M=20010;
const int N=300010;
const int dir[][2] = { {-1, 0}, {1, 0}, {0, -1}, {0, 1} };
typedef pair<int, int> pii;
int n,m,r,d[N];
int w[N],u[N],v[N],first[N],next[N];
bool vis[N];
void build()
{
    int i,j;
    MS(first,-1);
    for (i=0;i<m;i++) {
        scanf("%d%d%d",u+i,v+i,w+i);
        next[i]=first[u[i]];
        first[u[i]]=i;
    }
}
int dijkstra()
{
    int i,j;
    priority_queue<pii,vector<pii>,greater<pii> > q;
    for (i=1;i<=n;d[i++]=INF);
    d[1]=0;
    q.push(MP(d[1],1));
    while(!q.empty()) {
        int t=INF,pos;
        //for (j=1;j<=n;j++) if (!vis[j] && t>d[j]) t=d[pos=j];
        pii temp=q.top(); q.pop();
        pos=temp.S;
        if (d[pos]<temp.F) continue;
        vis[pos]=true;
        for (j=first[pos];j!=-1;j=next[j]) if (d[v[j]]>d[pos]+w[j]) {
            d[v[j]]=d[pos]+w[j];
            q.push(MP(d[v[j]],v[j]));
        }
    }
    return d[n];
}
int main()
{
    int i,j;
    while(~scanf("%d%d",&n,&m))
    {
        build();
        MS(vis,false);
        cout<<dijkstra()<<endl;
    }
}

poj 1364

      这道题表示读题读到吐血。。。弄了老半天都不知道它要干嘛，简直气的牙花子都撮出来了。

      题意简单的说就是给出n和m，n表示一个数字序列的长度，m表示有多少个约束条件。每个约束条件形如si ni oi ki，当oi为gt时，表示a[si]+a[si+1]+....+a[si+ni]>ki，当oi为lt时，表示a[si]+a[si+1]+....+a[si+ni]<ki。求是否存在一个长度为n的数字序列能满足所有约束条件。

     首先把子序列的和用前n项和之差来表示，并且将>转换为>=，将<转换为<=，即si ni gt ki转换为sum[si+ni] - sum[si-1] >= ki-1，si ni lt ki转换为sum[si+ni] - sum[si-1] <= ki-1，然后建图，判断是否有负权回路。若有，则不存在，否则存在。

    代码如下：

   

#include <cstdio>
#include <stack>
#include <set>
#include <iostream>
#include <string>
#include <vector>
#include <queue>
#include <functional>
#include <cstring>
#include <algorithm>
#include <cctype>
#include <string>
#include <map>
#include <iomanip>
#include <cmath>
#define LL long long
#define ULL unsigned long long
#define SZ(x) (int)x.size()
#define Lowbit(x) ((x) & (-x))
#define MP(a, b) make_pair(a, b)
#define MS(arr, num) memset(arr, num, sizeof(arr))
#define PB push_back
#define F first
#define S second
#define ROP freopen("input.txt", "r", stdin);
#define MID(a, b) (a + ((b - a) >> 1))
#define LC rt << 1, l, mid
#define RC rt << 1|1, mid + 1, r
#define LRT rt << 1
#define RRT rt << 1|1
#define BitCount(x) __builtin_popcount(x)
#define BitCountll(x) __builtin_popcountll(x)
#define LeftPos(x) 32 - __builtin_clz(x) - 1
#define LeftPosll(x) 64 - __builtin_clzll(x) - 1
const double PI = acos(-1.0);
const int INF = 0x3f3f3f3f;
using namespace std;
const double eps = 1e-8;
const int MAXN = 300 + 10;
const int MOD = 1000007;
const int M=20010;
const int N=300010;
const int dir[][2] = { {-1, 0}, {1, 0}, {0, -1}, {0, 1} };
typedef pair<int, int> pii;
int n,m,si,ni,ki,d[N];
char s[N];
int w[N],u[N],v[N],first[N],next[N],cnt[N];
bool vis[N];
void add(int a,int b,int c,int i)
{
    u[i]=a; v[i]=b; w[i]=c;
    next[i]=first[a];
    first[a]=i;
}
void build()
{
    int i,j;
    MS(first,-1);
    for (i=0;i<m;i++){
        scanf("%d %d %s %d",&si,&ni,s,&ki);
        int a=si+ni,b=si-1;
        if (s[0]=='g') add(a,b,-ki-1,i);
        else add(b,a,ki-1,i);
    }
    for (j=0;j<=n;j++,i++) add(n+1,j,0,i); // 将n+1作为源点加入图中，这里0不能作为源点，因为图中已经包括了0这个节点
}
bool spfa()
{
    int i,j;
    queue<int> q;
    q.push(n+1);
    MS(vis,false);
    MS(cnt,0);
    for (i=0;i<=n;i++) d[i]=INF;
    d[n+1]=0; cnt[n+1]++;
    while(!q.empty())
    {
        int pos=q.front(); q.pop();
        vis[pos]=false;
        for (i=first[pos];i!=-1;i=next[i]){
            if (d[v[i]]>d[u[i]]+w[i]){
                d[v[i]]=d[u[i]]+w[i];
                if (!vis[v[i]]){
                    vis[v[i]]=true;
                    if (++cnt[v[i]]>n+1) return false;
                    q.push(v[i]);
                }
            }
        }
    }
    return true;
}
int main()
{
    int i,j;
    while(~scanf("%d",&n),n)
    {
        scanf("%d",&m);
        build();
        if (spfa()) puts("lamentable kingdom"); // 存在
        else puts("successful conspiracy"); // 不存在
    }
}

     poj 3169

    题意：有n个奶牛，编号从1到n，按顺序站为一排等待进食。在排队过程中，有好感的两个奶牛之间距离不能超过给定值，没有好感的两个奶牛之间距离不能小于某个值，求1号奶牛与n号奶牛之间的最大距离。

    输入给出n，ml，md，表示有n个奶牛，其中ml对奶牛相互有好感，md对奶牛相互没有好感。接着ml行，每行表示哪两个奶牛相互有好感以及给定距离值，再是md行，表示哪两个奶牛相互没有好感以及给定距离值。

   设d[i]表示i号奶牛距离1号奶牛的距离，于是对有好感的奶牛有d[j]-d[i]<=k（i<j，k是给定距离值），没有好感的奶牛有d[j]-d[i]>=k，由于是求最大值，将d[j]-d[i]>=k转换为d[i]-d[j]<=-k(表示有向弧<j,i>的权值为-k)，这样就有了ml+md个形式一样的约束条件。

   不过需要注意的是题目中还有一个隐含的约束条件：d[i]<=d[j]，即d[i]-d[j]<=0。于是约束条件总共有ml+md+n-1个。

   找出所有约束条件后，剩下的就是建图求最短路了。

代码如下：

#include <cstdio>
#include <stack>
#include <set>
#include <iostream>
#include <string>
#include <vector>
#include <queue>
#include <functional>
#include <cstring>
#include <algorithm>
#include <cctype>
#include <string>
#include <map>
#include <iomanip>
#include <cmath>
#define LL long long
#define ULL unsigned long long
#define SZ(x) (int)x.size()
#define Lowbit(x) ((x) & (-x))
#define MP(a, b) make_pair(a, b)
#define MS(arr, num) memset(arr, num, sizeof(arr))
#define PB push_back
#define F first
#define S second
#define ROP freopen("input.txt", "r", stdin);
#define MID(a, b) (a + ((b - a) >> 1))
#define LC rt << 1, l, mid
#define RC rt << 1|1, mid + 1, r
#define LRT rt << 1
#define RRT rt << 1|1
#define BitCount(x) __builtin_popcount(x)
#define BitCountll(x) __builtin_popcountll(x)
#define LeftPos(x) 32 - __builtin_clz(x) - 1
#define LeftPosll(x) 64 - __builtin_clzll(x) - 1
const double PI = acos(-1.0);
const int INF = 0x3f3f3f3f;
using namespace std;
const double eps = 1e-8;
const int MAXN = 300 + 10;
const int M=20010;
const int N=50010;
const int dir[][2] = { {-1, 0}, {1, 0}, {0, -1}, {0, 1} };
typedef pair<int, int> pii;
int n,ml,md,d[N];
struct edge
{
    int u,v,w;
}e[N];
int bellman()
{
    int i,j;
    for (i=1;i<=n;i++) d[i]=INF;
    d[1]=0;
    for (i=0;i<n;i++){
        for (j=0;j<ml+md+n-1;j++){
            if (d[e[j].u]!=INF && d[e[j].v]>d[e[j].u]+e[j].w) {
                d[e[j].v]=d[e[j].u]+e[j].w;
                //if (i==n-1) return -1;
            }
        }
    }
    if (d[1]<0) return -1;
    else if (d[n]==INF) return -2;
    else return d[n];
}
int main()
{
    int i,j;
    while(~scanf("%d%d%d",&n,&ml,&md))
    {
        for (i=0;i<ml;i++) scanf("%d%d%d",&e[i].u,&e[i].v,&e[i].w);
        for (j=i;j<i+md;j++) { scanf("%d%d%d",&e[j].v,&e[j].u,&e[j].w); e[j].w=-e[j].w; }
        for (i=1;i<=n-1;i++) {
            e[j].u=i+1;
            e[j].v=i;
            e[j].w=0;
            j++;
        }
        int t=bellman();
        cout<<t<<endl;
    }
}






流沙-岁月
关注
codevs 1138
dianning8393的博客
09-14 60
题意：略 解法：看代码（这是c++） #include<stdio.h> #include<algorithm> #define L long long #define N 200010 using namespace std; int w[N],v[N],l[N],r[N],n,m; L s1[N],s2[N],ans=1ll...
C/C++中__builtin_popcount()的使用及原理
热门推荐
dlvector的博客
11-10 2万+
这个函数功能：返回输入数据中，二进制中‘1’的个数。 对于不同的使用类型，可以采用采用以下函数： __builtin_popcount = int __builtin_popcountl = long int __builtin_popcountll = long long 1. 二分法，源码采用的方法 主要思路是：将相邻两位相加，可以实现用二进制来表示输入数据中‘1’的个数。然后依次将上半部分和下半部分相和并实现计数。 unsigned popcount (unsigned u) { u = (
__builtin_popcount
冷月寒笙的博客
06-19 720
__builtin_popcount __builtin_popcountll 该函数的主要作用是计算一个数字的二进制中有多少个1，返回值就是其中1的个数。 它使用一张基于表的方法来进行位搜索，因此这个操作的执行效率很高 ...
__builtin_popcountll指令
dongyu_1989的博客
01-03 3664
#include #include #include typedef unsigned long long u64; typedef unsigned int u32; inline u32 popcount64_1(u64 x) { return __builtin_popcountll(x); } inline u32 popcount64_2(u64 x) {
差分约束系统例题详解[总结].pdf
10-12
差分约束系统(Differential Constraint System, 简称DCS)是一种在算法和数学规划中常见的模型，用于处理一系列变量间的关系。在这个例题中，我们面对的是一个序列的子序列约束问题，需要找到满足所有约束条件的序列...
Bellman-Ford算法与差分约束系统.ppt
07-12
Bellman-Ford算法与差分约束系统
差分约束系统题解1
08-03
【差分约束系统】是一种用于解决一系列线性不等式约束问题的模型，常用于最优化问题，如求解最短路径、最长路径或者判断是否存在满足条件的解。在这个模型中，我们通常需要将问题转化为寻找图中的最短路径或最长路径...
C++ __builtin_popcount函数作用
m0_54615144的博客
08-02 872
_builtin_popcount函数是系统自带的一个返回值是int/long/long long二进制'1'的个数的函数。拿__builtin_popcount（unsigned int n）举个例子。
GCC内建函数C语言实现__builtin_popcountll，__builtin_ctzll，__buitlin_clzll
最新发布
SongYuLong
09-16 470
C语言编程有时需要用到__builtin_popcountll，__builtin_ctzll，__buitlin_clzll等GCC内建函数，这里给出了简单的实现代码，算法实现肯定不是最优的但可以使用，大端模式CPU没有验证！
__builtin_xxx指令学习【3】__builtin_popcount & __builtin_popcountll
hanhan的博客
04-08 873
__builtin_popcount是GCC和Clang编译器提供的一个内置函数，用于计算一个整数中二进制位为1的个数。该函数的使用背景是在一些位运算和计算机视觉等领域中，需要对二进制数据进行处理和分析
C语言popcount函数,关于c ++：我可以使用带符号的整数作为__builtin_popcount()的参数吗？...
weixin_33212263的博客
05-21 355
这个网站：https://gcc.gnu.org/onlinedocs/gcc/Other-Builtins.html说它是为无符号整数定义的。 在某些情况下，将其用于带符号的int会导致错误的结果吗？也许。 如果说它仅是为未签名定义的，则签名可能有效。 不要依靠它。为什么在调用它之前不只是将其转换为unsigned int(该操作完全由语言定义)？ 您是否打算传递负数，如果是，您想要什么答案？_...
glibc的几个有用的处理二进制位的内置函数
weixin_30664539的博客
04-19 189
— Built-in Function: int __builtin_ffs (unsigned int x)Returns one plus the index of the least significant 1-bit of x, or if x is zero, returns zero. 返回右起第一个‘1’的位置。— Built-in Function: int __built...
__builtin_popcount 分析
无尽星空--帝王铠
12-28 801
glibc— Built-in Function: int __builtin_popcount (unsigned int x) Returns the number of 1-bits in x. 返回‘1’的个数。 int f(int s,int n) //s源数字，n是数字长度 { int cnt = 0; for(int i = 0;i < n;i++)
CF lesson 并查集
micaudience的博客
09-18 106
求最右 支持两种操作： 删除节点操作 询问某一节点右侧没有被删除的第一个节点。 //直线 const int N = 1e6 + 10; int f[N]; int findf(int x) { if(f[x] == x) return x; f[x] = findf(f[x]); return f[x]; } int main() { int n, m; scanf("%d%d", &n, &m); for(int i = 1; i <= n; i++) f[i] =
error C3861
qq_39620453的博客
06-01 3941
 问题1：error C3861: “__builtin_popcount”:  找不到标识符解决方案1（亲测可用）：Add these lines to math_functions.cpp to fix __builtin_popcount and __builtin_popcountl error:#define __builtin_popcount __popcnt #define __b...
C语言popcount函数,C/C++中__builtin_popcount()的使用及原理详解
weixin_29292535的博客
05-21 5996
__builtin_popcount()用于计算一个 32 位无符号整数有多少个位为1Counting out the bits可以很容易的判断一个数是不是2的幂次：清除最低的1位(见上面)并且检查结果是不是0.尽管如此，有的时候需要直到有多少个被设置了，这就相对有点难度 了。GCC有一个叫做__builtin_popcount的内建函数，它可以精确的计算1的个数。尽管如此，不同于__builti...
差分约束系统详解（转化为最短路）
10-02 2992
转自：http://www.cnblogs.com/void/archive/2011/08/26/2153928.html
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值