P5960【模板】差分约束

【模板】差分约束

题目描述

给出一组包含 m m m 个不等式,有 n n n 个未知数的形如:

{ x c 1 − x c 1 ′ ≤ y 1 x c 2 − x c 2 ′ ≤ y 2 ⋯ x c m − x c m ′ ≤ y m \begin{cases} x_{c_1}-x_{c'1}\leq y_1 \\x{c_2}-x_{c'2} \leq y_2 \\ \cdots\\ x{c_m} - x_{c'_m}\leq y_m\end{cases} ⎩ ⎨ ⎧xc1−xc1′≤y1xc2−xc2′≤y2⋯xcm−xcm′≤ym

的不等式组,求任意一组满足这个不等式组的解。

输入格式

第一行为两个正整数 n , m n,m n,m,代表未知数的数量和不等式的数量。

接下来 m m m 行,每行包含三个整数 c , c ′ , y c,c',y c,c′,y,代表一个不等式 x c − x c ′ ≤ y x_c-x_{c'}\leq y xc−xc′≤y。

输出格式

一行, n n n 个数,表示 x 1 , x 2 ⋯ x n x_1 , x_2 \cdots x_n x1,x2⋯xn 的一组可行解,如果有多组解,请输出任意一组,无解请输出 NO

样例 #1

样例输入 #1

复制代码
3 3
1 2 3
2 3 -2
1 3 1

样例输出 #1

复制代码
5 3 5

提示

样例解释

{ x 1 − x 2 ≤ 3 x 2 − x 3 ≤ − 2 x 1 − x 3 ≤ 1 \begin{cases}x_1-x_2\leq 3 \\ x_2 - x_3 \leq -2 \\ x_1 - x_3 \leq 1 \end{cases} ⎩ ⎨ ⎧x1−x2≤3x2−x3≤−2x1−x3≤1

一种可行的方法是 x 1 = 5 , x 2 = 3 , x 3 = 5 x_1 = 5, x_2 = 3, x_3 = 5 x1=5,x2=3,x3=5。

{ 5 − 3 = 2 ≤ 3 3 − 5 = − 2 ≤ − 2 5 − 5 = 0 ≤ 1 \begin{cases}5-3 = 2\leq 3 \\ 3 - 5 = -2 \leq -2 \\ 5 - 5 = 0\leq 1 \end{cases} ⎩ ⎨ ⎧5−3=2≤33−5=−2≤−25−5=0≤1

数据范围

对于 100 % 100\% 100% 的数据, 1 ≤ n , m ≤ 5 × 1 0 3 1\leq n,m \leq 5\times 10^3 1≤n,m≤5×103, − 1 0 4 ≤ y ≤ 1 0 4 -10^4\leq y\leq 10^4 −104≤y≤104, 1 ≤ c , c ′ ≤ n 1\leq c,c'\leq n 1≤c,c′≤n, c ≠ c ′ c \neq c' c=c′。

评分策略

你的答案符合该不等式组即可得分,请确保你的答案中的数据在 int 范围内。

如果并没有答案,而你的程序给出了答案,SPJ 会给出 There is no answer, but you gave it,结果为 WA;

如果并没有答案,而你的程序输出了 NO,SPJ 会给出 No answer,结果为 AC;

如果存在答案,而你的答案错误,SPJ 会给出 Wrong answer,结果为 WA;

如果存在答案,且你的答案正确,SPJ 会给出 The answer is correct,结果为 AC。

解法

SPFA

cpp 复制代码
#include <bits/stdc++.h> 
using namespace std;
struct road {
	int v,w;
};
vector<road> v[10086];
int n,m,ta,tb,tc,dis[10086],fang[10086];
bool vis[10086];
queue <int> q;
int main() {
	scanf("%d%d",&n,&m);
	for (int i=1;i<=m;i++) {
		scanf("%d%d%d",&ta,&tb,&tc);
		v[tb].push_back({ta,tc});
	}
	for (int i=1;i<=n;i++) {
		v[0].push_back({i,0});
		dis[i]=0x7fffffff;
	}
	q.push(0);
	dis[0]=0;
	while (!q.empty()) {
		int u=q.front();
		q.pop();
		vis[u]=0;
		int l=v[u].size();
		for (int i=0;i<l;i++) {
			if (dis[u]+v[u][i].w<dis[v[u][i].v]) {
				dis[v[u][i].v]=dis[u]+v[u][i].w;
				if (!vis[v[u][i].v]) {
					fang[v[u][i].v]++;
					if (fang[v[u][i].v]>n+1) {
						cout <<"NO";
						return 0;
					}
					q.push(v[u][i].v);
					vis[v[u][i].v]=1;
				}
			}
		}
	}
	for (int i=1;i<=n;i++) {
		cout <<dis[i] <<" ";
	}
	return 0;
}
相关推荐
tanyongxi6617 分钟前
C++ 特殊类设计与单例模式解析
java·开发语言·数据结构·c++·算法·单例模式
qq_5139704418 分钟前
力扣 hot100 Day76
算法·leetcode·职场和发展
Moshow郑锴1 小时前
机器学习相关算法:回溯算法 贪心算法 回归算法(线性回归) 算法超参数 多项式时间 朴素贝叶斯分类算法
算法·机器学习·回归
liulilittle1 小时前
BFS寻路算法解析与实现
开发语言·c++·算法·宽度优先·寻路算法·寻路
剪一朵云爱着1 小时前
PAT 1065 A+B and C (64bit)
算法·pat考试
喜欢吃燃面2 小时前
C++算法竞赛:位运算
开发语言·c++·学习·算法
项目申报小狂人2 小时前
算法应用上新!自适应更新策略差分进化算法求解球形多飞行器路径规划问题,附完整MATLAB代码
开发语言·算法·matlab
PAK向日葵8 小时前
【算法导论】PDD 0817笔试题题解
算法·面试
地平线开发者11 小时前
ReID/OSNet 算法模型量化转换实践
算法·自动驾驶