一.什么是差分约束?
差分约束就是n元一次不等式的集合,通常用来求解最大值/最小值问题,或者解n元一次不等式。
二.差分约束的必备条件
1.拥有一个可以跑到全图各个点的超级源点
2.不等式的形式必须固定。
(1) 如求最大值:要用最短路
x1 <= 2
x2 <= 3 等
(2) 求最小值:要用最长路
x1 >= 2
x2 >= 3等
3.无解情况,最短路有负环,最长路有正环
差分约束模板
//最大值,要用最短路
//因为x1 <= 2
// x2 <= 3
// x3 <= 5
//那么最后的解一定是2,满足上界的最小值的
//最小值,要用最长路
//因为
// x1 >= 2
// x2 >= 3
// x3 >= 5
//那么最后的解一定是5,也就是满足下界的最大值
#include "bits/stdc++.h"
using namespace std;
const int N = 100010;
struct node{
int to;
int nt;
int wi;
}e[N * 3];
int he[N],cnt[N],vs[N],idx,stk[N * 3];
int n,m,op;
long long dis[N];
void add(int u,int v,int w)
{
e[++idx].nt = he[u];
e[idx].wi = w;
e[idx].to = v;
he[u] = idx;
}
long long spfa()
{
long long tot = 0;
memset(dis,-0x3f,sizeof(dis));
dis[0] = 0;vs[0] = true;stk[0] = 0;
int hh = 0,tt = 1;
while(hh < tt) {
//用栈优化负环
int v = stk[--tt];
vs[v] = false;
bool f = false;
for(int i = he[v];~i;i = e[i].nt) {
int j = e[i].to;
if(dis[j] < dis[v] + e[i].wi) {
dis[j] = dis[v] + e[i].wi;
cnt[j] = cnt[v] + 1;
//要判到n + 1,因为还要包括源点
if(cnt[j] >= n + 1) return -1;
if(!vs[j]) {
vs[j] = true;
stk[tt++] = j;
}
}
}
}
for(int i = 1;i <= n;i++) tot += dis[i];
return tot;
}
int main()
{
scanf("%d%d",&n,&m);
memset(he,-1,sizeof(he));
while(m--) {
int u,v;
scanf("%d%d%d",&op,&u,&v);
if(op == 1) add(u,v,0),add(v,u,0);
else if(op == 2) add(u,v,1);
else if(op == 3) add(v,u,0);
else if(op == 4) add(v,u,1);
else add(u,v,0);
}
//建立超级源点,这个1是根据题目要求的
for(int i = 1;i <= n;i++) add(0,i,1);
printf("%lld",spfa());
return 0;
}
三.列题分析
1.差分约束算法
(1)题目描述
(2)解题思路
这个题目只说让我们求一组可行解,那么我们只需要构造一下这个n元一次不等式组,然后建边跑一遍最短路即可。
比如第一行输入,我们可以这样转换一下1号点<=2号点+3,依此类推。我们设超级源点为0,让0向所有点连一条为0的边,转换为不等式为i号点<=0号点+0,相当于我们求这个不等式小于等于0的最大解。
(3)代码实现
#include "bits/stdc++.h"
using namespace std;
const int N = 5010,M = 15030;
int dis[N],q[N],idx,cnt[N],vs[N],he[N];
struct node {
int to;
int nt;
int wi;
}e[M];
int n,m;
void add(int a,int b,int c)
{
e[++idx].to = b;
e[idx].nt = he[a];
e[idx].wi = c;
he[a] = idx;
}
bool spfa()
{
int hh = 0,tt = 1;
memset(dis,0x3f,sizeof(dis));
vs[0] = true,dis[0] = 0,q[0] = 0;
while(hh != tt) {
int v = q[hh++];
if(hh == N) hh = 0;
vs[v] = false;
for(int i = he[v];~i;i = e[i].nt) {
int j = e[i].to,w = e[i].wi;
if(dis[j] > dis[v] + w) {
dis[j] = dis[v] + w;
cnt[j] = cnt[v] + 1;
if(cnt[j] >= n + 1) return true;
if(!vs[j]) {
vs[j] = true;
q[tt++] = j;
if(tt == N) tt = 0;
}
}
}
}
return false;
}
int main()
{
memset(he,-1,sizeof(he));
scanf("%d%d",&n,&m);
for(int i = 1;i <= m;i++) {
int u,v,w;
scanf("%d%d%d",&u,&v,&w);
add(v,u,w);
add(0,i,0);
}
if(!spfa()) {
for(int i = 1;i <= n;i++) printf("%d ",dis[i]);
puts("");
}
else puts("NO");
return 0;
}
2.糖果
(1)题目描述
(2)解题思路
已知题目要求最少的,那么我们果断判定跑最长路,然后不等式全部转换为>=类型的。
对于操作1:A>=B,B >= A
对于操作2:B>=A+1
对于操作3:A>=B
对于操作4:A>=B+1
对于操作5:B>=A
因此就没有什么难度了
(3)代码实现
//最大值,要用最短路
//因为x1 <= 2
// x2 <= 3
// x3 <= 5
//那么最后的解一定是2,满足上界的最小值的
//最小值,要用最长路
//因为
// x1 >= 2
// x2 >= 3
// x3 >= 5
//那么最后的解一定是5,也就是满足下界的最大值
#include "bits/stdc++.h"
using namespace std;
const int N = 100010;
struct node{
int to;
int nt;
int wi;
}e[N * 3];
int he[N],cnt[N],vs[N],idx,stk[N * 3];
int n,m,op;
long long dis[N];
void add(int u,int v,int w)
{
e[++idx].nt = he[u];
e[idx].wi = w;
e[idx].to = v;
he[u] = idx;
}
long long spfa()
{
long long tot = 0;
memset(dis,-0x3f,sizeof(dis));
dis[0] = 0;vs[0] = true;stk[0] = 0;
int hh = 0,tt = 1;
while(hh < tt) {
//用栈优化负环
int v = stk[--tt];
vs[v] = false;
for(int i = he[v];~i;i = e[i].nt) {
int j = e[i].to;
if(dis[j] < dis[v] + e[i].wi) {
dis[j] = dis[v] + e[i].wi;
cnt[j] = cnt[v] + 1;
//要判到n + 1,因为还要包括源点
if(cnt[j] >= n + 1) return -1;
if(!vs[j]) {
vs[j] = true;
stk[tt++] = j;
}
}
}
}
for(int i = 1;i <= n;i++) tot += dis[i];
return tot;
}
int main()
{
scanf("%d%d",&n,&m);
memset(he,-1,sizeof(he));
while(m--) {
int u,v;
scanf("%d%d%d",&op,&u,&v);
if(op == 1) add(u,v,0),add(v,u,0);
else if(op == 2) add(u,v,1);
else if(op == 3) add(v,u,0);
else if(op == 4) add(v,u,1);
else add(u,v,0);
}
//建立超级源点,这个1是根据题目要求的
for(int i = 1;i <= n;i++) add(0,i,1);
printf("%lld",spfa());
return 0;
}
3.区间
(1)题目描述
(2)解题思路
我们利用前缀和思想,用s[i]表示前i个被选中的个数,
那么肯定有
1.s[i] >= s[i - 1]
2.s[i - 1] <= s[i] -1
3.s[b]-s[a] >= c
又因为数据范围为50000,因此我们预处理到50001,最后答案就是s[50001]
(3)代码实现
#include "bits/stdc++.h"
using namespace std;
const int N = 50010,M = 150010;
int q[N],dis[N],he[N],idx,vs[N];
struct node {
int to;
int nt;
int wi;
}e[M];
void add(int a,int b,int c)
{
e[++idx].nt = he[a];
e[idx].to = b;
e[idx].wi = c;
he[a] = idx;
}
void spfa()
{
memset(dis,-0x3f,sizeof(dis));
dis[0] = 0;vs[0] = true;q[0] = 0;
int hh = 0,tt = 1;
while(hh != tt) {
int v = q[hh++];
if(hh == N) hh = 0;
vs[v] = false;
for(int i = he[v];~i;i = e[i].nt) {
int j = e[i].to,w = e[i].wi;
if(dis[j] < dis[v] + w) {
dis[j] = dis[v] + w;
if(!vs[j]) {
vs[j] = true;
q[tt++] = j;
if(tt == N) tt = 0;
}
}
}
}
}
int main()
{
int n;
scanf("%d",&n);
memset(he,-1,sizeof(he));
for(int i = 1;i < N;i++) {
add(i - 1,i,0);
add(i,i - 1,-1);
}
while(n --) {
int a,b,c;
scanf("%d%d%d",&a,&b,&c);
a ++,b ++;
add(a - 1,b,c);
}
spfa();
printf("%d",dis[50001]);
return 0;
}
4.排队布局
(1)题目描述
(2)解题思路
这里求最大距离,那么我们用最短路,考虑如何建立超级源点和如何转换不等式关系。
对于操作1:A <= B + L
对于操作2:B <= A - D
对于超级源点:i <= (i + 1) + 0
(3)代码实现
#include "bits/stdc++.h"
using namespace std;
const int N = 1010,M = 21010,INF = 0x3f3f3f3f;
int vs[N],cnt[N],he[N],idx,q[N],dis[N];
int n,m1,m2;
struct node {
int to;
int nt;
int wi;
}e[M];
void add(int a,int b,int c)
{
e[++idx].nt = he[a];
e[idx].to = b;
e[idx].wi = c;
he[a] = idx;
}
bool spfa(int st)
{
memset(dis,0x3f,sizeof(dis));
memset(vs,0,sizeof(vs));
memset(cnt,0,sizeof(cnt));
int hh = 0,tt = 0;
for(int i = 1;i <= st;i++) {
q[tt++] = i;
dis[i] = 0;
vs[i] = true;
}
while(hh != tt) {
int v = q[hh++];
if(hh == N) hh = 0;
vs[v] = false;
for(int i = he[v];~i;i = e[i].nt) {
int j = e[i].to,w = e[i].wi;
if(dis[j] > dis[v] + w) {
dis[j] = dis[v] + w;
cnt[j] = cnt[v] + 1;
if(cnt[j] >= n) return true;
if(!vs[j]) {
vs[j] = true;
q[tt++] = j;
if(tt == N) tt = 0;
}
}
}
}
return false;
}
int main()
{
int a,b,c;
scanf("%d%d%d",&n,&m1,&m2);
memset(he,-1,sizeof(he));
for(int i = 1;i < n;i++) add(i + 1,i,0);
while(m1 --) {
scanf("%d%d%d",&a,&b,&c);
if(a > b) swap(a,b);
add(a,b,c);
}
while(m2 --) {
scanf("%d%d%d",&a,&b,&c);
if(a > b) swap(a,b);
add(b,a,-c);
}
//判断负环
if(spfa(n)) puts("-1");
else {
spfa(1);
if(dis[n] == INF) puts("-2");
else printf("%d\n",dis[n]);
}
return 0;
}
5.雇佣收银员
(1)题目描述
(2)解题思路
我们用s[i]表示第i时刻上岗人数
则:
1.// 上岗人数不能负数,即 s[i]−s[i−1]>=0
2.// 实际上岗人数不能超过申请人数,即 s[i]−s[i−1]<=worker[i]
由于存在环
3.// 当 i>=8 时,s[i]−s[i−8]>=neet[i]
4.// 当 i<=7 时,s[i]+s[24]−s[16+i]>=need[i]
我们需要枚举24这个特殊的点,不然我们会多一个变量
5.add(0,24,k),add(24,0,-k);
(3)代码实现
#include "bits/stdc++.h"
using namespace std;
const int N = 30;
int need[N],worker[N],n,dis[N],q[N];
int idx,vs[N],cnt[N],he[N];
struct node {
int nt;
int to;
int wi;
}e[110];
void add(int a,int b,int c)
{
e[++idx].nt = he[a];
e[idx].to = b;
e[idx].wi = c;
he[a] = idx;
}
void build(int k)
{
memset(he,-1,sizeof(he));
idx = 0;
add(0,24,k),add(24,0,-k);
// 当 i>=8 时,s[i]−s[i−8]>=neet[i]
// 当 i<=7 时,s[i]+s[24]−s[16+i]>=need[i]
for(int i = 1;i <= 7;i++) add(16 + i,i,need[i] - k);
for(int i = 8;i <= 24;i++) add(i - 8,i,need[i]);
// 上岗人数不能负数,即 s[i]−s[i−1]>=0
// 实际上岗人数不能超过申请人数,即 s[i]−s[i−1]<=worker[i]
for(int i = 1;i <= 24;i++) {
add(i - 1,i,0);
add(i,i - 1,-worker[i]);
}
}
bool spfa(int k)
{
build(k);
memset(dis,-0x3f,sizeof(dis));
memset(vs,0,sizeof(vs));
memset(cnt,0,sizeof(cnt));
dis[0] = 0,vs[0] = true,q[0] = 0;
int hh = 0,tt = 1;
while(hh != tt) {
int v = q[hh++];
if(hh == N) hh = 0;
vs[v] = false;
for(int i = he[v];~i;i = e[i].nt) {
int j = e[i].to,w = e[i].wi;
if(dis[j] < dis[v] + w) {
dis[j] = dis[v] + w;
cnt[j] = cnt[v] + 1;
if(cnt[j] >= 25) return true;
if(!vs[j]) {
vs[j] = true;
q[tt++] = j;
if(tt == N) tt = 0;
}
}
}
}
return false;
}
int main()
{
int T;
scanf("%d",&T);
while(T --) {
for(int i = 1;i <= 24;i++) scanf("%d",&need[i]);
scanf("%d",&n);
memset(worker,0,sizeof(worker));
for(int i = 1;i <= n;i++) {
int x;
scanf("%d",&x);
worker[x + 1] ++;
}
bool ok = false;
for(int i = 0;i <= n;i++)
if(!spfa(i)) {
printf("%d\n",i);
ok = true;
break;
}
if(!ok) puts("No Solution");
}
return 0;
}
6.倍杀测量者
(1)题目描述
(2)解题思路
对于第一个操作:A >= B * (T - K)
对于第二个操作:A >= B * 1.0 / (T + K)
如果出现负环表示一定有人跳舞,二分最大T值即可。
(3)代码实现
#include "bits/stdc++.h"
using namespace std;
const double eps = 1e-10;
const int N = 1010,M = 3030;
double dis[N],score[N];
int q[N << 1],cnt[N],vs[N],idx,he[N];
struct node {
int to;
double wi;
int nt;
}e[M];
struct oper {
int op;
int p1;
int p2;
double k;
}pt[N];
int n,s,t;
int hh,tt;
void add(int a,int b,double c)
{
e[++idx].nt = he[a];
e[idx].wi = c;
e[idx].to = b;
he[a] = idx;
}
void build(double k)
{
memset(he,-1,sizeof(he));
idx = 0;
for(int i = 1;i <= s;i++) {
if(pt[i].op == 1) add(pt[i].p2,pt[i].p1,(pt[i].k - k));
else add(pt[i].p2,pt[i].p1,1.0 / (pt[i].k + k));
}
for(int i = 1;i <= n;i++) if(score[i]) add(0,i,score[i]),add(i,0,1.0 / score[i]);
}
bool spfa(double k)
{
build(k);
hh = 0,tt = 0;
memset(vs,0,sizeof(vs));
memset(cnt,0,sizeof(cnt));
queue<int> que;
for(int i = 0;i <= n;i++) dis[i] = 1,que.push(i);
vs[0] = true;
//手写队列没过,stl队列过了???
while(que.size()) {
int v = que.front();
que.pop();
vs[v] = false;
for(int i = he[v];~i;i = e[i].nt) {
int j = e[i].to;
double w = e[i].wi;
if(dis[j] < dis[v] * w) {
dis[j] = dis[v] * w;
cnt[j] = cnt[v] + 1;
if(cnt[j] >= n + 1) return true;
if(!vs[j]) {
vs[j] = true;
que.push(j);
}
}
}
}
return false;
}
int main()
{
int x1;
double x2;
double l = 0,r = 1e6,ans = -1;
scanf("%d%d%d",&n,&s,&t);
for(int i = 1;i <= s;i++) {
scanf("%d%d%d%lf",&pt[i].op,&pt[i].p1,&pt[i].p2,&pt[i].k);
if(pt[i].op == 1) r = min(r,pt[i].k);
}
for(int i = 1;i <= t;i++) {
scanf("%d%lf",&x1,&x2);
score[x1] = x2;
}
while(r - l > eps) {
double mid = (l + r) / 2;
if(spfa(mid)) ans = l = mid;
else r = mid;
}
printf("%.10lf",ans);
return 0;
}
列题链接: