CodeForces 1197E : Culture Code 最短路计数 + 思维

最新推荐文章于 2021-11-15 21:42:57 发布

匿枫

最新推荐文章于 2021-11-15 21:42:57 发布

阅读量105

点赞数 1

分类专栏：大学ACM

本文链接：https://blog.csdn.net/tlyzxc/article/details/119465917

版权

大学ACM 专栏收录该内容

162 篇文章 5 订阅

订阅专栏

本文详细解析了一道关于最短路径的算法题，通过构造网络图，利用拓扑排序和最短路算法寻找满足条件的最小路径。作者首先将输入数据进行排序，然后构建图并应用特殊边权处理，最后通过spfa算法求解最短路并计算方案数。文章适合对图论和算法有一定了解的读者深入研究。

摘要由CSDN通过智能技术生成

传送门

题意

在这里插入图片描述

分析

这个题，太妙了，太妙了
首先我们去转化一下给的式子
$in_{i1} + (in_{i2} - out_{i1}) + ··· + (in_{ik} - out_{ik - 1})$
等价于
$in_{ik} + \sum\nolimits_{p = 1}^{k - 1}{(in_{ip} - out_{ip})}$
题目要求的是求这个的最小值，所以，等价于
$\sum\nolimits_{p = 1}^{k - 1}{(in_{ip} - out_{ip})}$ 最小
我们把所有点按照 $in_{i}$ 从小到大排序，设置点0，然后相邻两个点之间连边，权重为 $in_{i} - in{i - 1}$ ,这样， $1 - i$ 之间的的路径和就是 $in_{i} - in_{i - 1} + in_{i -1} - in_{i - 2} ... + in_{1} - in_{0}$ ，最后等价于 $in_{i}$ ，也就是 $in_{ik}$
那么后面那一坨我们怎么处理呢，我们可以把每个点 $i$ 后面第一个满足 $in_{j} > out_{i}$ 的点找到，然后 $i - j$ 之间连一条边，权值为 $in_{j} - out_{i}$ ，这样， $\sum\nolimits_{p = 1}^{k - 1}{(in_{ip} - out_{ip})}$ 也处理出来了
剩下的就是跑最短路了，并且因为要求方案数，所以还要最短路计数
因为这张图内存在边权为0的边，所有 $d i j k s t r a$ 和 $s p f a$ 都无法处理这个问题，但因为是拓扑图，所以我们可以顺着推
最后需要注意的是，如果 $i$ 点的 $o u t$ 小于 $n$ 点的 $i n$ 的话，说明还可以套娃，说明当前答案不符合要求

代码

#pragma GCC optimize(3)
#include <bits/stdc++.h>
#define debug(x) cout<<#x<<":"<<x<<endl;
#define dl(x) printf("%lld\n",x);
#define di(x) printf("%d\n",x);
#define _CRT_SECURE_NO_WARNINGS
#define pb push_back
#define mp make_pair
#define all(x) (x).begin(),(x).end()
#define fi first
#define se second
#define SZ(x) ((int)(x).size())
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> PII;
typedef vector<int> VI;
const ll INF = 0x3f3f3f3f3f3f3f3f;
const int N = 2e5 + 10,M = N << 1;
const ll mod = 1000000007;
const double eps = 1e-9;
const double PI = acos(-1);
template<typename T>inline void read(T &a) {
	char c = getchar(); T x = 0, f = 1; while (!isdigit(c)) {if (c == '-')f = -1; c = getchar();}
	while (isdigit(c)) {x = (x << 1) + (x << 3) + c - '0'; c = getchar();} a = f * x;
}
int gcd(int a, int b) {return (b > 0) ? gcd(b, a % b) : a;}
struct Node{
	int out,in;
	bool operator <(Node const &A) const{
		return in < A.in;
	}
}a[N];
int n;
int h[N],ne[M],e[M],w[M],idx;
ll d[N];
bool st[N];
ll f[N];

void add(int x,int y,int z){
	ne[idx] = h[x],e[idx] = y,w[idx] = z,h[x] = idx++;
}

int check(int x){
	int l = 1,r = n;
	while(l < r){
		int mid = (l + r) >> 1;
		if(a[mid].in >= x) r = mid;
		else l = mid + 1;
	}
	if(a[l].in >= x) return l;
	return -1;
}

void spfa(){
	memset(d,0x3f,sizeof d);
	d[0] = 0;
	f[0] = 1;
	for(int i = 0;i <= n;i++)
		for(int j = h[i];~j;j = ne[j]){
			int v = e[j];
			if(d[v] > d[i] + w[j]){
				d[v] = d[i] + w[j];
				f[v] = f[i];
			}
			else if(d[v] == d[i] + w[j]) f[v] = (f[v] + f[i]) % mod;
		}
}

int main() {
	memset(h,-1,sizeof h);
	read(n);
	for(int i = 1;i <= n;i++) read(a[i].out),read(a[i].in);
	sort(a + 1,a + 1 + n);
	for(int i = 0;i < n;i++) add(i,i + 1,a[i + 1].in - a[i].in);
	for(int i = 1;i <= n;i++){
		int p = check(a[i].out);
		if(p == -1) continue;
		add(i,p,a[p].in - a[i].out);
	}
	spfa();
	ll ans = 0,res = INF;
	for(int i = 1;i <= n;i++) 
		if(a[i].out > a[n].in) res = min(res,d[i]);
	for(int i = 1;i <= n;i++)
		if(a[i].out > a[n].in && res == d[i]) 
			ans = (ans + f[i]) % mod;
	dl(ans);
	return 0;
}