NOIP2017 D1T3 逛公园
发现 \(K\) 很小,不妨设置一个 \(O(NK)\) 的 \(DP\)。
发现可行的最短路必须满足是 \(d <= dis <= d + K\)。
由逆向思维,则是从某点出发,可以消耗 \(K\) 个单位的冗余长度,最终到达 \(n\)。
如何快速的计算出有走这条边冗余长度呢
首先建反向图跑 \(Dijkstra\),求出 \(dis[i]\) 表示从 \(i\) 到 \(n\) 的最短路距离。
假设有一条边 \((u, v)\),边长为 \(w\)
那么我现在从 \(u\) 到 \(n\) 的预算距离是 \(dis[v] + w\),最短距离是 \(dis[u]\)。
那么多走的就是 \(dis[v] + w - dis[u]\)。
设计 \(DP\) 状态表示:
\(f[i][j]\) 表示从 \(i\) 到 \(n\),以及消耗了 \(j\) 个单位的冗余长度的方案数。
状态转移:
设有边 \((u, v)\),边长为 \(w\)
则有 \(f[u][j] += f[v][j - (dis[v] + w - dis[u])]\)。
边界 \(f[n][0] = 1\),其余为 \(0\)。
答案 \(\sum_{i = 0}^{K} f[1][i]\)。
无穷解的判断:
发现有无穷解,当且仅当:
有一个总权为 \(0\) 的环。
我们可以在 \(DP\) 的时候搞一个 \(vis\) 数组判断,就不需要单独判无穷解了。
时间复杂度
整个过程用记忆化搜索实现,由于一共有 \(NK\) 个状态,每个点被枚举 \(K\) 次,即每条边总体被枚举 \(K\) 次。
所以复杂度 \(O(K(N + M))\)。
\(Tips:\)
- 可能走到 \(n\) 再折回去,所以 $u = n $ 时不能直接 \(return\)
- 可能存在经过 \(n\) 点的 \(0\) 环,所以到 \(n\) 点时顺便判一下环。
#include <cstdio>
#include <iostream>
#include <cstring>
#include <queue>
#include <vector>
#include <cstdlib>
using namespace std;
typedef pair<int, int> PII;
const int N = 100005, M = 400005, S = 51;
int n, m, K, P, f[N][S], dis[N];
int head[N], rhead[N], numE[2];
bool st[N], vis[N][S], ep = false;
struct E {
int next, v, w;
}e[M], r[M];
//建图
void inline add(E g[], int h[], int u, int v, int w, int p) {
g[++numE[p]] = (E) { h[u], v, w };
h[u] = numE[p];
}
//多组数据初始化
void inline init() {
ep = false;
memset(dis, 0x3f, sizeof dis);
memset(st, false, sizeof st);
numE[0] = numE[1] = 0;
memset(head, 0, sizeof head);
memset(rhead, 0, sizeof head);
memset(f, -1, sizeof f);
}
// Dijkstra 最短路
priority_queue<PII, vector<PII>, greater<PII> > q;
void inline dijkstra() {
dis[n] = 0; q.push(make_pair(0, n));
while(!q.empty()) {
PII u = q.top(); q.pop();
if(st[u.second]) continue;
st[u.second] = true;
for (int i = rhead[u.second]; i; i = r[i].next) {
int v = r[i].v;
if(dis[u.second] + r[i].w < dis[v]) {
dis[v] = dis[u.second] + r[i].w;
q.push(make_pair(dis[v], v));
}
}
}
}
//记忆化搜索
int dfs(int u, int j) {
if(vis[u][j]) { ep = true; return 0; }
if(f[u][j] != -1) return f[u][j];
vis[u][j] = true;
int &val = f[u][j] = 0;
for (int i = head[u]; i; i = e[i].next) {
int v = e[i].v, w = e[i].w;
//消耗的冗余长度折算
int k = j - (dis[v] + w - dis[u]);
if(0 <= k && k <= K) (val += dfs(v, k)) %= P;
if(ep) return 0;
}
vis[u][j] = false;
if(u == n && j == 0) val = 1;
return val;
}
int main() {
int T; scanf("%d", &T);
while(T--) {
init();
scanf("%d%d%d%d", &n, &m, &K, &P);
for (int i = 1, u, v, w; i <= m; i++) {
scanf("%d%d%d", &u, &v, &w);
add(e, head, u, v, w, 0); add(r, rhead, v, u, w, 1);
}
dijkstra();
int ans = 0;
for (int i = 0; i <= K; i++) {
memset(vis, false, sizeof vis);
(ans += dfs(1, i)) %= P;
if(ep) break;
}
if(ep) puts("-1");
else printf("%d\n", ans);
}
}
NOIP2017 D1T3 逛公园的更多相关文章
- [luogu P3953] [noip2017 d1t3] 逛公园
[luogu P3953] [noip2017 d1t3] 逛公园 题目描述 策策同学特别喜欢逛公园.公园可以看成一张$N$个点$M$条边构成的有向图,且没有 自环和重边.其中1号点是公园的入口,$N ...
- NOIP2017 D1T3逛公园
DP+最短路 两遍最短路判零环 DP转移f[i][j] 到点i的距离比最短路多j时的方案数 #include<bits/stdc++.h> using namespace std; ; s ...
- 【NOIP2017】逛公园 拆点最短路+拓扑(记忆化搜索
题目描述 策策同学特别喜欢逛公园.公园可以看成一张N个点M条边构成的有向图,且没有 自环和重边.其中1号点是公园的入口,N号点是公园的出口,每条边有一个非负权值, 代表策策经过这条边所要花的时间. 策 ...
- 【NOIP2017】逛公园(最短路图,拓扑排序,计数DP)
题意: 策策同学特别喜欢逛公园. 公园可以看成一张 N 个点 M 条边构成的有向图,且没有自环和重边.其中 1 号点是公园的入口, N 号点是公园的出口,每条边有一个非负权值,代表策策经过这条边所要花 ...
- NOIP2017:逛公园
Sol 发现\(NOIP2017\)还没\(AK\)??? 赶紧改 考场上明明打出了\(DP\),没时间了,没判环,重点是没初始化数组,爆\(0\) \(TAT\) 先最短路,然后\(f[i][j]\ ...
- 【NOIP2017】逛公园 最短路+DP
诶,去年场上不会处理$0$的环,只拿了$60$有点可惜. 我们先不管边边权为$0$的边. 我们先跑一次最短路,令$dis[u]$表示从$1$至$u$的最短路的长度. 那么根据题目的要求,从起点走到$u ...
- LOJ2316. 「NOIP2017」逛公园【DP】【最短路】【思维】
LINK 思路 因为我想到的根本不是网上的普遍做法 所以常数出奇的大,而且做法极其暴力 可以形容是带优化的大模拟 进入正题: 首先一个很显然的思路是如果在合法的路径网络里面存在零环是有无数组解的 然后 ...
- 【NOIP2017】逛公园 D1 T3
记忆化搜索 跑一次反向的最短路求出MinDis(u,n)MinDis(u,n)MinDis(u,n) f[u][k]f[u][k]f[u][k]表示dis(u,n)<=MinDis(u,n)+d ...
- 【LOJ2316】「NOIP2017」逛公园
[题目链接] [点击打开链接] [题目概括] 对给定\(K\),起点\(1\)到终点\(n\)中对长度为\([L,L+K]\)的路径计数. \(L\)为\(1\)到\(n\)的最短路长度. [思路要点 ...
随机推荐
- 1. 线性DP 300. 最长上升子序列 (LIS)
最经典单串: 300. 最长上升子序列 (LIS) https://leetcode-cn.com/problems/longest-increasing-subsequence/submission ...
- 微服务通信之feign的配置隔离
前言 由上文我们知道针对某一个Feign接口,我们可以给他设置特定的配置类.那如果现在有一个服务,我们只想对A服务配置一个拦截器拦截请求而不影响其他服务,那应该怎么做呢? 一.feign接口配置 由前 ...
- Innodb之线程独享内存
引用链接: https://blog.csdn.net/miyatang/article/details/54881547 https://blog.csdn.net/wyzxg/article/de ...
- ceph osd tree的可视化
前言 很久没有处理很大的集群,在接触一个新集群的时候,如果集群足够大,需要比较长的时间才能去理解这个集群的结构,而直接去看ceph osd tree的结果,当然是可以的,这里是把osd tree的结构 ...
- 关于backfill参数建议
前言 在做一个比较满的集群的扩容的时候,遇到了一些问题,在这里做下总结,一般来说很难遇到,扩容要趁早,不然出的问题都是稀奇古怪的一些问题 建议 环境一般来说在70%左右就需要考虑扩容了,这个时候的扩容 ...
- webug第八关:CSRF
第八关:CSRF 用tom用户登陆,进入更改密码界面 用burp构造csrf页面 将构造好的页面放在服务器上,然后点击, admin的密码被更改
- ABBYY FineReader 14创建PDF文档功能解析
使用ABBYY FineReader,您可以轻松查看和编辑任何类型的 PDF,真的是一款实至名归的PDF编辑转换器,您知道的,它能够保护.签署和编辑PDF文档,甚至还可以创建PDF文档,本文和小编一起 ...
- lca(lowestCommonAncestor)
- LNMP 一键安装脚本
这个脚本是使用shell编写,为了快速在生产环境上部署lnmp/lamp/lnmpa(Linux.Nginx/Tengine/OpenResty.MySQL/MariaDB/Percona.PHP), ...
- 虚拟机VM14.X安装Mac10.12启动出现问题的解决方法
虚拟机安装Mac系统,会出现的问题太多,于是乎变记录下来,方便以后使用或者方便大家解决问题. 一:VM14.X安装Mac10.12虚拟机,启动出现下面无限重启问题 解决方法: 亲测有效 在OS X 1 ...