P5960 差分约束算法模板
差分约束
差分约束,一般用来解决有\(n\)个未知数,\(m\)个不等式方程的问题,形如:
\ x_{a_1}-x_{b_1}\leq y_1\\
\ x_{a_2}-x_{b_2}\leq y_2\\
\ \cdots\\
\ x_{a_m}-x_{b_m}\leq y_m\\
\end{cases}
\]
可以判断有没有解,以及给出一组解
简单观察可以知道,每个未知数的系数都为\(1\),且不等式一边是两个未知数相减,另一边是一个常数
为了达到这种形式,一般都需要将题目中的关系进行变形
模板题:P5960 【模板】差分约束算法
如何求解
考虑对于每个不等式\(x_i-x_j\leq y\),可以变形成\(x_i\leq x_j+y\)
这个式子,可以容易联想到最短路算法中的不等式\(dis(v)\leq dis(u)+w\),(容易吗?
那么,可以用最短路的方法来求解,从\(j\)向\(i\)联一条边权为\(y\)的边(\(i\)对应终点\(v\))
然后,就要新建一个虚拟节点\(n+1\),从它向每个节点连一个权值为\(0\)的边,\(dis(i)\leq dis(n+1)+0\)
求出的各项\(dis(i)\)其实就是对应的\(x_i\)
又由于\(dis(n+1)=0\),所以此时的\(dis(i)\)均小于\(0\),那么我们求出来的这组解,就是所有小于\(0\)的解里,最大的一组
那么如果产生了负环,就会出现\(x<x\)的情况,然后肯定无解,证明应该不难,但不想证了
因为在一个没有负环的图上,最短路肯定不会经过多于\(n\)个点,也不会经过多于\(n-1\)条边
所以,我们用 spfa 来求解这个最短路问题,如果某个节点入队\(n\)次,就说明存在负环
还有一种,就是记录\(cnt_i\)为\(1\rightarrow i\)的最短路上的点的个数,如果\(cnt_i>n\)就说明存在负环
这种方法倒是更直观,据说跑起来也更快
当然,\(cnt_1=1,cnt_v=cnt_u+1\)
但是代码里放的是第一种
没错,spfa 肯定没死,他是唯一能处理负环的算法(起码我没听说过其它的
最短路
所以可以写出用最短路实现的代码:
#include<cstdio>
#include<queue>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<iomanip>
#include<cstring>
#define reg register
#define EN std::puts("")
#define LL long long
inline int read(){
int x=0,y=1;
char c=std::getchar();
while(c<'0'||c>'9'){if(c=='-') y=0;c=std::getchar();}
while(c>='0'&&c<='9'){x=x*10+(c^48);c=std::getchar();}
return y?x:-x;
}
int n,m;
int fir[5006],nex[10006],to[10006],w[10006],tot;
int cnt[5005],in[5005];
int dis[5005];
std::queue<int>q;
inline void add(int a,int b,int c){
to[++tot]=b;w[tot]=c;
nex[tot]=fir[a];fir[a]=tot;
}
inline int spfa(){
q.push(n+1);
std::memset(dis,0x3f,sizeof dis);
dis[n+1]=0;in[n+1]=1;
while(!q.empty()){
reg int u=q.front();q.pop();in[u]=0;
for(reg int v,i=fir[u];i;i=nex[i]){
v=to[i];
if(dis[v]>dis[u]+w[i]){
dis[v]=dis[u]+w[i];
if(!in[v]){
if(++cnt[v]==n) return 0;
in[v]=1;q.push(v);
}
}
}
}
return std::rand();
}
int main(){
n=read();m=read();
for(reg int i=1;i<=n;i++) add(n+1,i,0);
for(reg int a,b,c,i=1;i<=m;i++){
a=read();b=read();c=read();
add(b,a,c);
}
if(!spfa()) std::puts("NO");
else for(reg int i=1;i<=n;i++) std::printf("%d\n",dis[i]);
return 0;
}
最长路
其实不等式\(x_i-x_j\leq y\)还可以变形成\(x_j\geq x_i-y\)
所以就应该从\(i\)向\(j\)连一个权值为\(-y\)的边,然后跑最长路,判断正环
当然,虚拟节点也要建,从\(n+1\)连向\(i\),代表\(dis(i)\geq dis(n+1),dis(n+1)=0\)
因此,这种方法求出的是,当\(x_i\geq 0\),的最小解
所以代码是这样的:
#include<cstdio>
#include<queue>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<iomanip>
#include<cstring>
#define reg register
#define EN std::puts("")
#define LL long long
inline int read(){
int x=0,y=1;
char c=std::getchar();
while(c<'0'||c>'9'){if(c=='-') y=0;c=std::getchar();}
while(c>='0'&&c<='9'){x=x*10+(c^48);c=std::getchar();}
return y?x:-x;
}
int n,m;
int fir[5006],nex[10006],to[10006],w[10006],tot;
int cnt[5005],in[5005];
int dis[5005];
std::queue<int>q;
inline void add(int a,int b,int c){
to[++tot]=b;w[tot]=c;
nex[tot]=fir[a];fir[a]=tot;
}
inline int spfa(){
q.push(n+1);
for(reg int i=1;i<=n;i++) dis[i]=-5e8;
in[n+1]=1;
while(!q.empty()){
reg int u=q.front();q.pop();in[u]=0;
for(reg int v,i=fir[u];i;i=nex[i]){
v=to[i];
if(dis[v]<dis[u]+w[i]){
dis[v]=dis[u]+w[i];
if(!in[v]){
if(++cnt[v]==n) return 0;
in[v]=1;q.push(v);
}
}
}
}
return std::rand();
}
int main(){
n=read();m=read();
for(reg int i=1;i<=n;i++) add(n+1,i,0);
for(reg int a,b,c,i=1;i<=m;i++){
a=read();b=read();c=read();
add(a,b,-c);
}
if(!spfa()) std::puts("NO");
else for(reg int i=1;i<=n;i++) std::printf("%d\n",dis[i]);
return 0;
}
一些题目
在处理差分约束题目时,主要就是要把关系进行变形
看清楚边是\(i\rightarrow j\)还是\(j\rightarrow i\),边权是\(y\)还是\(-y\)
如果出现未知数系数不为\(1\)的情况,一般使用乘积最短路的形式,但这时候就要求没有常数,也就是\(x_i\leq kx_j\)
P1993 小K的农场,题解,比较经典的一题
P4926 [1007]倍杀测量者,这个就是刚才说的乘积最短路,但是具体的实现是要取\(\log\),然后变成加法的,洛谷上有一个名为“最短路”的题就是这个
SCOI2011 糖果
虫洞,uva上的一题
P5960 差分约束算法模板的更多相关文章
- 差分约束算法————洛谷P4878 [USACO05DEC] 布局
题目: 不难看出题意主要是给出ml+md个格式为xi-xj<=ak的不等式,xi-xj为i,j俩头牛的距离,要我们求x1-xn的最大值. 经过上下加减我们可以将这几个不等式化成x1-xn< ...
- 鉴于spfa基础上的差分约束算法
怎么搞? 1. 如果要求最大值 想办法把每个不等式变为标准x-y<=k的形式,然后建立一条从y到x权值为k的边,变得时候注意x-y<k =>x-y<=k ...
- P5960 【模板】差分约束算法
题目描述 给出一组包含 $m$ 个不等式,有 $n$ 个未知数的形如: 的不等式组,求任意一组满足这个不等式组的解. 输入格式 第一行为两个正整数 $n,m$,代表未知数的数量和不等式的数量. 接下来 ...
- P1993 小K的农场(差分约束)
小K的农场 题目描述 小K在MC里面建立很多很多的农场,总共n个,以至于他自己都忘记了每个农场中种植作物的具体数量了,他只记得一些含糊的信息(共m个),以下列三种形式描述: 农场a比农场b至少多种植了 ...
- 差分约束+spfa【模板】
相比dij,spfa优点是可处理含负边不含负圈的最短路问题,缺点是算法复杂度不太好[貌似可以使用两种优化.LLL和SLF] 差分约束就是将一些不等式转化为图中的带权边,然后求解最短路或最长路的方法 洛 ...
- bzoj 2330 [SCOI2011]糖果 差分约束模板
题目大意 幼儿园里有N个小朋友,lxhgww老师现在想要给这些小朋友们分配糖果,要求每个小朋友都要分到糖果.但是小朋友们也有嫉妒心,总是会提出一些要求,比如小明不希望小红分到的糖果比他的多,于是在分配 ...
- 算法复习——差分约束(ssoi种树)
题目: 题目描述 为了绿化乡村,H 村积极响应号召,开始种树了. H 村里有 n 幢房屋,这些屋子的排列顺序很有特点,在一条直线上.于是方便起见,我们给它们标上 1-n .树就种在房子前面的空地上. ...
- POJ-3159(差分约束+Dijikstra算法+Vector优化+向前星优化+java快速输入输出)
Candies POJ-3159 这里是图论的一个应用,也就是差分约束.通过差分约束变换出一个图,再使用Dijikstra算法的链表优化形式而不是vector形式(否则超时). #include< ...
- POJ 3159 Candies(差分约束,最短路)
Candies Time Limit: 1500MS Memory Limit: 131072K Total Submissions: 20067 Accepted: 5293 Descrip ...
随机推荐
- 「一闻秒懂」你了解goroutine和channel吗?
开源库「go home」聚焦Go语言技术栈与面试题,以协助Gopher登上更大的舞台,欢迎go home~ 背景介绍 大家都知道进程是操作系统资源分配的基本单位,有独立的内存空间,线程可以共享同一个进 ...
- Array(数组)对象-->数组的删除
1.数组的删除: 用delete操作符删除特定的元素 删除元素的位置只是被留空了,为undefined值 举例:删除下面数组中的第二个元素 var arr = [1,2,3,4,5]; /*删除第二个 ...
- 文档根元素 "beans" 必须匹配 DOCTYPE 根 "null"
文档根元素 "beans" 必须匹配 DOCTYPE 根 "null" (2011-11-20 21:26:41) 转载▼ 标签: 杂谈 分类: spring- ...
- AJ学IOS(26)UI之iOS抽屉效果小Demo
AJ分享,必须精品 先看效果 实现过程 第一步,把三个view设置好,还有颜色 #warning 第一步 - (void)addChildView { // left UIView *leftView ...
- linux下的信号量PV操作进阶之路
一.同步和互斥机制 信号量 互斥锁 同步:指多个任务按照约定的先后次序相互配合来完成一件事情. 比如读线程等待写线程写完之后再去读. 二.信号量-P/V操作 P(s)含义: if(信号量>0) ...
- Python程序设计实验报告四:循环结构程序设计(设计型实验)
安徽工程大学 Python程序设计 实验报告 班级 物流191 姓名 姚彩琴 学号3190505129 成绩 日期 2020.4.8 指导老师 修宇 [实验名称 ...
- 数据结构之栈—强大的四则复杂运算计算器(超过windows自带的科学计算器)【中缀转后缀表达式】
比windows自带计算器还强的四则复杂运算计算器! 实测随机打出两组复杂算式:-7.5 * 6 / ( -2 + ( -6.5 - -5.22 ) )与7.5+-3*8/(7+2) windows ...
- D. Minimax Problem Codeforces 1288D binary_search+二进制
题目大意:n*m的矩阵中,找到两行数,可以形成两个一维数组,数组1的位置i和数组2的位置i去最大构成新数组b的元素b[i],最终目的要使数组b中最小的数尽可能的大 题解: m的范围是(1,8),比较小 ...
- Rescue BFS+优先队列 杭电1242
思路 : 优先队列 每次都取最小的时间,遇到了终点直接就输出 #include<iostream> #include<queue> #include<cstring> ...
- 详解 DatagramSocket类
(请观看本人博文 -- <详解 网络编程>) DatagramSocket 概述: 这类代表一个发送和接收数据包的插座. 该类是遵循 UDP协议 实现的一个Socket类. 数据报套接字发 ...