【NOIP2017】逛公园【最短路DAG】【dp】【拓扑排序】

题意:给一张帯权有向图,求 111nnn 长度不超过最短路长度 +k+k+k 的路径条数 模 PPP。有无数条输出 −1-11

n≤105,m≤2×105,k≤50n\leq 10^5,m\leq 2\times 10^5,k\leq 50n105,m2×105,k50,边权非负

先往最短路图想

发现 kkk 很小,考虑硬上 dp。

f(u,k)f(u,k)f(u,k) 表示从 uuunnn 多走 kkk 的路径数。

转移枚举下一步往哪里走。

f(u,k)=∑(u,v)∈Ef(v,k−(w(u,v)+disu−disv))f(u,k)=\sum_{(u,v)\in E}f(v,k-(w(u,v)+dis_u-dis_v))f(u,k)=(u,v)Ef(v,k(w(u,v)+disudisv))

其中 disudis_udisu111uuu 最短路长度。

在没有零边的基础上,右边这块 w(u,v)+disu−disvw(u,v)+dis_u-dis_vw(u,v)+disudisv000 当且仅当 (u,v)(u,v)(u,v) 在最短路图上。因为这是个 DAG,所以在最短路图上跑拓扑排序,然后枚举 kkk 更新就可以了。这样可以拿到 707070 分的好成绩。

有零边时可能会出现 disn+kdis_n+kdisn+k 内可达的零环,这样答案为无穷大。如果拓扑排序后点 iii 入度不为 000 说明在零环上,此时如果 disi+disi′≤disn+kdis_i+dis'_i \leq dis_n+kdisi+disidisn+k 就输出 −1-11。(dis′dis'dis 为到 nnn 的最短距离)

否则的话反正都走不到这里,直接把伪的拓扑序继续往后做。

复杂度 O(nlog⁡n+nk)O(n\log n+nk)O(nlogn+nk)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#include <algorithm>
#include <utility>
#include <queue>
#define MAXN 100005
#define MAXM 200005
using namespace std;
int n,m,k,P;
inline int add(const int& x,const int& y){return x+y>=P? x+y-P:x+y;}
inline int read()
{int ans=0;char c=getchar();while (!isdigit(c)) c=getchar();while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();return ans;
}
struct edge{int u,v,w;}e[MAXM];
int head[MAXN],nxt[MAXM],cnt;
inline void addnode(int u,int v,int w)
{e[++cnt]=(edge){u,v,w};nxt[cnt]=head[u];head[u]=cnt;
}
template<typename T,typename cmp=less<T> >
struct heap
{T val[MAXM];T* end;heap(){end=val;}inline bool empty(){return end==val;}inline void push(const T& v){*(end++)=v;push_heap(val,end,cmp());}inline T pop(){pop_heap(val,end,cmp());return *(--end);}	
};
typedef pair<int,int> pi;
heap<pi,greater<pi> > q;
int dis[MAXN],disn[MAXN],deg[MAXN];
inline void dij(int* dis,int s=1)
{memset(dis,0x3f,sizeof(disn));dis[s]=0;q.push(make_pair(0,s));while (!q.empty()){int u=q.pop().second;for (int i=head[u];i;i=nxt[i])if (dis[u]+e[i].w<dis[e[i].v])dis[e[i].v]=dis[u]+e[i].w,q.push(make_pair(dis[e[i].v],e[i].v));}
}
int lis[MAXN],vis[MAXN],tim;
void dfs(int u)
{vis[u]=1;for (int i=head[u];i;i=nxt[i])if (dis[u]+e[i].w==dis[e[i].v]){++deg[e[i].v];if (!vis[e[i].v]) dfs(e[i].v);}
}
inline void topsort()
{queue<int> q;tim=0;q.push(1);while (!q.empty()){int u=q.front();q.pop();lis[++tim]=u;for (int i=head[u];i;i=nxt[i])if (dis[u]+e[i].w==dis[e[i].v]&&(--deg[e[i].v]==0))q.push(e[i].v);}
}
int f[55][MAXN],u[MAXM],v[MAXM],w[MAXM];
int main()
{for (int T=read();T;T--){cnt=tim=0;memset(head,0,sizeof(head));memset(nxt,0,sizeof(nxt));memset(deg,0,sizeof(deg));memset(vis,0,sizeof(vis));n=read(),m=read(),k=read(),P=read();for (int i=1;i<=m;i++) u[i]=read(),v[i]=read(),w[i]=read(),addnode(v[i],u[i],w[i]);dij(disn,n);cnt=0;memset(head,0,sizeof(head));memset(nxt,0,sizeof(nxt));		for (int i=1;i<=m;i++) addnode(u[i],v[i],w[i]);dij(dis);dfs(1);topsort();int ans=0;for (int i=1;i<=n;i++) if (deg[i]&&dis[i]+disn[i]<=dis[n]+k) {puts("-1");goto end;}memset(f,0,sizeof(f));f[0][n]=1;for (int t=0;t<=k;t++)for (int x=tim;x>=1;x--){int u=lis[x];for (int i=head[u];i;i=nxt[i]){int del=e[i].w-dis[e[i].v]+dis[u];if (t>=del) f[t][u]=add(f[t][u],f[t-del][e[i].v]);}}for (int i=0;i<=k;i++) ans=add(ans,f[i][1]);printf("%d\n",ans);end:;}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/314935.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Codeforces Round #599 (Div. 2) E. Sum Balance 图转换 + 子集dp + 环

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先我们知道如果所有数的和summodk!0sum\bmod k!0summodk!0那么此时无解&#xff0c;否则我们设needsum/kneedsum/kneedsum/k。 看到kkk这么小&#xff0c;自然的想到是否能状压&#xff…

资深开发者们是如何读书的?---线下读书会记录

读了多少书&#xff1f;怎么读书&#xff1f;——请问大家最近读书了吗&#xff1f;一年大概读了多少本书。有超过20本&#xff0c;超过30本&#xff0c;超过50本吗&#xff1f;在7月14日由长沙互联网活动基地和长沙.NET技术社区组织的小型社区活动&#xff0c;《开发者爱读书•…

【APIO2018】Duathlon 铁人两项 【圆方树】

题意&#xff1a;给一张 nnn 点 mmm 边的简单无向图&#xff0c;求有多少个三元组 (s,c,f)(s,c,f)(s,c,f) &#xff0c;满足存在一条从 sss 到 fff 经过 ccc 的简单路径。 n≤105,m≤2105n\leq 10^5,m\leq 2\times 10^5n≤105,m≤2105 首先这个 “经过 ccc 的简单路径” &…

Codeforces Round #598 (Div. 3) E. Yet Another Division Into Teams dp + 输出方案

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个长度为nnn的序列aaa&#xff0c;你需要将其分成若干组&#xff0c;每组的价值为max⁡(ai)−min(ai)\max(a_i)-min(a_i)max(ai​)−min(ai​)&#xff0c;求如何分组才能使得代价最小&#xff0c;输出…

架构杂谈《三》

一致性问题前面的《架构杂谈一》和《架构杂谈二》 杂谈了从服务化到微服务架构的演进&#xff0c;并肯定了服务化和微服务架构是一脉相承的。微服务在服务化架构的基础上&#xff0c;对服务化的细节和方案进行了优化和细化&#xff0c;重点突出了无中心化管理的微服务架构&…

【SDOI2018】战略游戏【圆方树】【虚树】

题意&#xff1a;给一张 nnn 点 mmm 边的连通无向图&#xff0c;qqq 次询问&#xff0c;每次给出一个点集 SSS &#xff0c;求有多少个不在 SSS 中的点满足删除后 SSS 中存在两个点不连通。 n≤105,m≤2105,∑∣S∣≤2105n\leq 10^5,m\leq 2\times 10^5,\sum |S|\leq 2\times 1…

Codeforces Round #598 (Div. 3) F. Equalizing Two Strings 思维 + 逆序对

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你两个长度为nnn的串a,ba,ba,b&#xff0c;每次可以同时翻转a,ba,ba,b中任意一段长度为L(1≤L≤n)L(1\le L\le n)L(1≤L≤n)的子串&#xff0c;问能否通过若干次操作使两个串相同。 思路&#xff1a; 首…

腾讯发布 Tencent SCF Toolkit VS Code 插件,轻松开发无服务器云函数

在之前的文章中&#xff0c;我们提到了亚马逊、谷歌、IBM 等大厂都上了 Visual Studio Code 的船。阿里巴巴也在近日发布了 Aliyun Serverless VS Code 插件。近期&#xff0c;腾讯也上了 VS Code 的船&#xff0c;发布了基于 VS Code 的 Tencent SCF Toolkit 插件&#xff0c;…

Codeforces Round #726 (Div. 2) D. Deleting Divisors 博弈

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个数nnn&#xff0c;有两个人博弈&#xff0c;每次可以将nnn减去一个nnn的因子&#xff0c;这个因子不能为111或nnn。当不能操作的人输掉游戏。问你先手赢还是后手赢。 思路&#xff1a; 这个题多写…

【POI2007】OSI-Axes of Symmetry【计算几何】【manacher】

题意&#xff1a;给一个 nnn 个点的多边形&#xff0c;求对称轴个数。 n≤105n\leq 10^5n≤105 显然对称轴一定在顶点或边的中点上。 但你 n2n^2n2 枚举完全没有一点能过的样子。 冷静分析&#xff0c;发现有 “中点”&#xff0c;“对称轴”&#xff0c;很自然个鬼地想到了…

端到端测试实践:Jenkins集成TestCafe

上一篇《对产品质量的一点思考》中说到自动化测试的重要性&#xff0c;本文简单介绍下怎样在实际项目中实现端到端测试的自动化&#xff0c;在这里我们使用的端到端测试工具是TestCafe。环境Jenkisn&#xff1a;2.183TestCafe&#xff1a;1.3.0为什么采用TestCafe做自动化测试要…

Codeforces Round #727 (Div. 2) F. Strange Array 线段树 + 区间合并 + 排序优化

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个长度为nnn的数组&#xff0c;对每个位置iii求一个最大价值&#xff0c;价值计算方式如下&#xff1a;选择一个包含iii的[l,r][l,r][l,r]&#xff0c;让后将其拿出来排序&#xff0c;之后价值就是当前…

【CF487E】Tourists【圆方树】【树链剖分】【multiset】

题意&#xff1a;给一张 nnn 点 mmm 边的连通无向图&#xff0c;点帯权&#xff0c;qqq 次操作&#xff1a; 修改一个点的权值。询问两点间所有简单路的最小权值的最小值。 n,m,q≤105n,m,q\leq 10^5n,m,q≤105 显然建出圆方树然后询问路径最小值。多半要树链剖分了。 对于方…

再见Jenkins,从Gitlab代码提交到k8s服务持续交付只需七毛三

Gitlab runner 快速搭建CICD pipeline背景日常开发中&#xff0c;相信大家已经做了很多的自动化运维环境&#xff0c;用的最多的想必就是利用Jenkins实现代码提交到自动化测试再到自动化打包&#xff0c;部署全流水线Jenkins在devops担任了很重要的角色&#xff0c;但是另一方面…

Codeforces Round #726 (Div. 2) F. Figure Fixing 二分图 + 思维

传送门 文章目录题意&#xff1a;思路题意&#xff1a; 给你一张nnn个点mmm条边的图&#xff0c;每个点都有一个当前值aia_iai​&#xff0c;目标值bib_ibi​&#xff0c;每次可以选择一条边(i,j)(i,j)(i,j)&#xff0c;将ai,aja_i,a_jai​,aj​都加上任意一个数&#xff0c;问…

【NOI2020】命运【树形dp】【线段树合并】

题意&#xff1a;给一棵 nnn 个点的树&#xff0c;每条边需要染成黑白两种颜色中的一种。给出 mmm 个条件&#xff0c;每个条件给出 u,vu,vu,v&#xff0c;其中 uuu 是 vvv 的祖先&#xff0c;要求 uuu 到 vvv 的链上至少一条黑边。求方案数 模 998244353998244353998244353。 …

.NET Core 3.0之深入源码理解Kestrel的集成与应用(二)

前言前一篇文章主要介绍了.NET Core继承Kestrel的目的、运行方式以及相关的使用&#xff0c;接下来将进一步从源码角度探讨.NET Core 3.0中关于Kestrel的其他内容&#xff0c;该部分内容&#xff0c;我们无需掌握&#xff0c;依然可以用好Kestrel&#xff0c;本文只是将一些内部…

Codeforces Round #726 (Div. 2) E2. Erase and Extend (Hard Version) 贪心

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个长度为nnn的串sss&#xff0c;你有两个操作可以使用&#xff1a; (1)(1)(1)从sss的结尾删除一个字母。 (2)sss(2)sss(2)sss。 让你通过若干次操作使其变成一个长度为kkk的串&#xff0c;且其字典序最…

【BZOJ3451】Normal【期望线性性】【点分治】【NTT卷积】

题意&#xff1a;随机分治中心点分治的期望操作次数 n≤3104n\leq 3\times 10^4n≤3104 即求点分树的 siz 之和的期望 即祖孙关系对数期望 考虑一有序点对 (u,v)(u,v)(u,v)&#xff0c;uuu 在点分树上是 vvv 祖先当且仅当 uuu 是 u∼vu\sim vu∼v 路径上第一个被选为分治中心…

译 | 使用Roslyn分析器高效编写更好的代码

原文&#xff1a;Mika Dumont翻译&#xff1a;汪宇杰Roslyn 是 .NET 编译器平台&#xff0c;即使在运行代码之前&#xff0c;它也能帮助您捕获 Bug。例如内置于 Visual Studio 中的Roslyn 拼写检查分析器。假设您正在创建一个静态方法&#xff0c;并将 static 单词拼错为 statc…