noip2017D1T3逛公园(拓扑图上dp,记忆化搜索)
QWQ前几天才刚刚把这个D1T3写完
看着题解理解了很久,果然我还是太菜了QAQ
题目大意就是
给你一个n个点,m条边的图,保证1能到达n,求从1到n的 (设1到n的最短路长度是d)路径长度在[d,d+k]之间的路径有多少条,答案要对p取膜
下面附上数据范围的大表哥!
首先对于30%的数据,我们可以直接跑最短路计数来实现QWQ
这里最短路计数就不作详细解释了!
一定注意的是 当更新dis[to[i]]时,要记得把ans[to[i]]赋值成ans[x] 千万不要手残写成1!!!
上代码
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
#define pa pair<int,int>
using namespace std;
const int maxn = 100010;
const int maxm = 400010;
inline int read()
{
int x=0,f=1;char ch=getchar();
while (!isdigit(ch)){ if (ch=='-') f=-1;ch=getchar();
}while (isdigit(ch)){ x=(x << 3)+(x << 1) + ch-'0';ch=getchar();
}return x*f;
}
struct Node{
int num;
int id;
};
priority_queue <pa,vector<pa>,greater<pa> > q;
int point[maxn],nxt[maxm],to[maxm],val[maxm];
int dis[maxn];
Node qq[1000100];
int tt[maxn],vis[maxn];
int aa[maxn];
int cnt,ans;
int n,m;
void addedge(int x,int y,int w)
{
nxt[++cnt]=point[x];
to[cnt]=y;
val[cnt]=w;
point[x]=cnt;
}
void dijkstra(int s,int pp)
{
for (int i=1;i<=n;i++)
dis[i]=2e9;
memset(aa,0,sizeof(aa));
aa[s]=1;
dis[s]=0;
memset(vis,0,sizeof(vis));
q.push(make_pair(0,s));
while (!q.empty())
{
int x = q.top().second;
q.pop();
if (vis[x]) continue;
vis[x]=1;
for (int i=point[x];i;i=nxt[i])
{
int p=to[i];
if (dis[p]>dis[x]+val[i])
{
dis[p]=dis[x]+val[i];
aa[p]=aa[x]%pp;
q.push(make_pair(dis[p],p));
}
else
if (dis[p]==dis[x]+val[i])
{
aa[p]=(aa[p]+aa[x])%pp;
}
}
}
}
int t,k,p;
void bfs(int d,int pp)
{
int head=0,tail=1;
qq[tail].id=1;
while (head<=tail)
{
head++;
int x=qq[head].id;
for (int i=point[x];i;i=nxt[i])
{
int p=to[i];
if (tt[p]>n) continue;
qq[++tail].id=p;
qq[tail].num=qq[head].num+val[i];
tt[p]++;
if (p==n&&qq[tail].num<=d+k)
{
ans=(ans+1)%pp;
}
}
}
}
int main()
{
t=read();
while (t--)
{
scanf("%d%d%d%d",&n,&m,&k,&p);
cnt=0;
ans=0;
memset(point,0,sizeof(point));
memset(qq,0,sizeof(qq));
memset(tt,0,sizeof(tt));
for (int i=1;i<=m;i++)
{
int u,v,w;
u=read();v=read();w=read();
addedge(u,v,w);
}
dijkstra(1,p);
printf("%d\n",aa[n]);
}
}
这是(修改后的)考场源代码QWQ可能有点丑陋
而对于其他数据QAQ emmmmmm
这个嘛~
我们就需要考虑dp
我这里用的dp状态是
f[i][j]表示从1到i这个点,比最短路长了j的方案数
对于一条边 u - > v QAQ我们不难发现
f[v][dis[u]+k+val[i]-dis[v]]+=f[u][k]; (0<=dis[u]+k+val[i]-dis[v]<=k)
好啦!这不就可以转移了嘛?
别急QWQ 貌似还有0环的问题。
这里就需要思考一下0环的性质
如果有0环的话.....这些边应该一定会出现在最短路图上吧,那么我们只需要在最短路图上跑拓扑排序~如果到最后发现无法构成DAG 那么应该就是有0环
同时拓扑排序也是为了在最短路上的点在后面的dp中,制定一个顺序
例如x->y->z 更新顺序一定是x y z
而对于0环上的点,如果dis[i]+disn[i](到n的最短路)<=dis[n]+k 那么它就可以无限制的更新下去(可以理解为一直在0环上,从而使方案数变为无限)
如果遇到这种情况 就直接输出-1了
下面dp的部分也没什么好说的了
枚举这个偏移量(就是比最短路长多少)
就是分成两部分,先更新最短路的点,然后再用当前的偏移量的u,去更新更大偏移量的v
上代码
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#define pa pair < int , int >
using namespace std;
inline int read()
{
int x=0,f=1;char ch=getchar();
while (!isdigit(ch)){if (ch=='-') f=-1;ch=getchar();}
while (isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return x*f;
}
//f[to[i]][dis[u]+k+w-dis[v]]+=f[u][k]
const int maxn = 1e5+1e2;
const int maxm = 5e5+1e2;
const int inf = 1e9;
int f[maxn][60];
int nxt[maxm],to[maxm],point[maxn],val[maxm];
queue<int> que;
int x[maxm],y[maxm],w[maxm];
int vis[maxn],dis[maxn],disn[maxn];
int in[maxn];
int n,m,k,p;
bool flag;
int cnt=0;
priority_queue< pa , vector<pa>,greater<pa> > q;
int mod;
void addedge(int x,int y,int w)
{
nxt[++cnt]=point[x];
to[cnt]=y;
val[cnt]=w;
point[x]=cnt;
}
void init()
{
cnt=0;
memset(point,0,sizeof(point));
memset(f,0,sizeof(f));
memset(in,0,sizeof(in));
flag=true;
}
int dijkstra(int s)
{
memset(vis,0,sizeof(vis));
for (register int i=1;i<=n;i++) dis[i]=inf;
q.push(make_pair(0,s));
dis[s]=0;
while (!q.empty())
{
int x = q.top().second;
q.pop();
if (vis[x]) continue;
vis[x]=1;
for (register int i=point[x];i;i=nxt[i])
{
int p = to[i];
if (dis[p]>dis[x]+val[i])
{
dis[p]=dis[x]+val[i];
q.push(make_pair(dis[p],p));
}
}
}
}
int dijkstran(int s)
{
memset(vis,0,sizeof(vis));
for (register int i=1;i<=n;i++) disn[i]=inf;
disn[s]=0;
q.push(make_pair(0,s));
while (!q.empty())
{
int x = q.top().second;
q.pop();
if (vis[x]) continue;
vis[x]=1;
for (register int i=point[x];i;i=nxt[i])
{
int p = to[i];
if (disn[p]>disn[x]+val[i])
{
disn[p]=disn[x]+val[i];
q.push(make_pair(disn[p],p));
}
}
}
}
int t;
int top[maxn];
int tmp;
void tpsort()
{
tmp=0;
for (register int x=1;x<=n;++x)
for (register int i=point[x];i;i=nxt[i])
{
int p = to[i];
if (dis[p]==dis[x]+val[i]) in[p]++;
}
for (register int i=1;i<=n;++i)
if (in[i]==0) que.push(i),top[++tmp]=i;
while (!que.empty())
{
int x = que.front();
que.pop();
for (register int i=point[x];i;i=nxt[i])
{
int p = to[i];
if (dis[p]==dis[x]+val[i])
{
in[p]--;
if (in[p]==0)
{
que.push(p);
top[++tmp]=p;
}
}
}
}
}
void dp() //之所以要拓扑排序,是因为在更新最短路上的点的时候,有一个先后顺序 就好比是u->v 必须先算u,再算v
{
f[1][0]=1;
for (register int i=0;i<=k;i++)
{
//更新最短路上的点
// cout<<1<<endl;
for (register int j=1;j<=tmp;++j)
{
int x=top[j];
for (register int ii=point[x];ii;ii=nxt[ii])
{
int p= to[ii];
if (dis[p]==dis[x]+val[ii]) f[p][i]=(f[x][i]+f[p][i])%mod;
}
}
for (register int x=1;x<=n;++x)
for (register int ii=point[x];ii;ii=nxt[ii])
{
int p= to[ii];
int now = dis[x]+val[ii]+i-dis[p];
if (dis[p]!=dis[x]+val[ii] && now<=k) f[p][now]=(f[x][i]+f[p][now])%mod;
}
}
// cout<<2<<endl;
}
int main()
{
cin>>t;
while (t--)
{
init();
n=read();m=read();k=read();mod=read();
for (register int i=1;i<=m;++i)
{
x[i]=read();
y[i]=read();
w[i]=read();
}
for (register int i=1;i<=m;++i)
addedge(y[i],x[i],w[i]);
dijkstran(n);
init();
for (register int i=1;i<=m;++i)
addedge(x[i],y[i],w[i]);
dijkstra(1);
//for (int i=1;i<=n;i++)
// cout<<dis[i]<<" ";
//cout<<endl;
//for (int i=1;i<=n;i++)
//// cout<<disn[i]<<" ";
//cout<<endl;
tpsort();
for (register int i=1;i<=n;++i)
if (in[i]>0 && dis[i]+disn[i]<=dis[n]+k)
{
printf("-1\n");
flag=false;
break;
}
if (!flag) continue;
dp();
int ans=0;
for (register int i=0;i<=k;++i)
ans=(ans+f[n][i])%mod;
printf("%d\n",ans);
}
return 0;
}
QAQ这种写法,代码常数特别大,需要卡常,才能A掉
QWQ
这个题的另一个做法,记忆化搜索
正着建图求好dis后,
然后反向建图,将f[1][0]=1
从n开始做记忆化搜索
f[x][k]=f[x][k]+f[to[i]][k-(dis[to[i]]+val[i]-dis[x])]
一定一定一定一定注意!!!!!!!
这种方法要将f数组初始化成-1
在dfs 的时候
用一个中间变量保存f的值
最后再赋值,详情看代码吧QWQ被这个点坑了很久
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
using namespace std;
inline int read()
{
int x=0,f=1;char ch=getchar();
while (!isdigit(ch)) {if (ch=='-') f=-1;ch=getchar();}
while (isdigit(ch)) {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
return x*f;
}
const int maxn = 2e5+1e2;
const int maxm = 1e6+1e2;
int point[maxn],nxt[maxm],to[maxm],val[maxm];
int dis[maxn],vis[maxn];
int f[maxn][61];
int n,m,k,mod;
int x[maxm],y[maxm],w[maxm];
queue<int> q;
int g[maxn][61];
bool flag;
int kk;
int cnt;
void addedge(int x,int y,int w)
{
nxt[++cnt]=point[x];
to[cnt]=y;
val[cnt]=w;
point[x]=cnt;
}
void init()
{
cnt=0;
flag=true;
memset(point,0,sizeof(point));
memset(f,-1,sizeof(f));
memset(g,0,sizeof(g));
}
int spfa(int s)
{
memset(vis,0,sizeof(vis));
memset(dis,127/3,sizeof(dis));
vis[s]=1;
dis[s]=0;
q.push(s);
while (!q.empty())
{
int x = q.front();
q.pop();
vis[x]=0;
for (int i=point[x];i;i=nxt[i])
{
int p = to[i];
if (dis[p]>dis[x]+val[i])
{
dis[p]=dis[x]+val[i];
if (!vis[p])
{
vis[p]=1;
q.push(p);
}
}
}
}
}
int dfs(int x,int k)
{
//cout<<1<<endl;
int ret=0;
if (g[x][k]) {
flag=false;
return 0;
}
if (f[x][k]!=-1) return f[x][k]; //这里如果写成if (f[x][k]) 会re 因为f[x][k]==0的状态有很多
g[x][k]=1;
if (!flag) return 0;
for (int i=point[x];i;i=nxt[i])
{
int p = to[i];
int cnt=k-(dis[p]+val[i]-dis[x]);
if (cnt<0 || cnt>kk) continue;
ret=(ret+dfs(p,cnt))%mod;
if (!flag) return 0;
}
if (!flag) return 0;
g[x][k]=0;
if (x==1 && k==0)
{
f[x][k]=1;
}
else
f[x][k]=ret;
return f[x][k];
}
int t;
int main()
{
cin>>t;
while (t--)
{
init();
n=read(),m=read(),k=read(),mod=read();
for (int i=1;i<=m;i++)
{
x[i]=read();
y[i]=read();
w[i]=read();
addedge(x[i],y[i],w[i]);
}
spfa(1);
//for (int i=1;i<=n;i++) cout<<dis[i]<<" ";
// cout<<endl;
init();
for (int i=1;i<=m;i++) addedge(y[i],x[i],w[i]);
int ans=0;
kk=k;
for (int i=0;i<=k;i++)
{
//memset(g,0,sizeof(g));
int tmp = dfs(n,i)%mod;
if (!flag) break;
ans=(ans+tmp)%mod;
}
if (!flag) {
cout<<-1<<endl;
continue;
}
else
{
printf("%d\n",ans);
}
}
return 0;
}
noip2017D1T3逛公园(拓扑图上dp,记忆化搜索)的更多相关文章
- Luogu 3953[NOIP2017] 逛公园 堆优化dijkstra + 记忆化搜索
题解 首先肯定是要求出单源最短路的,我用了堆优化dijikstra ,复杂度 mlogm,值得拥有!(只不过我在定义优先队列时把greater 打成了 less调了好久 然后我们就求出了$i$到源点的 ...
- [NOIP2017] 逛公园 (最短路,动态规划&记忆化搜索)
题目链接 Solution 我只会60分暴力... 正解是 DP. 状态定义: \(f[i][j]\) 代表 \(1\) 到 \(i\) 比最短路长 \(j\) 的方案数. 那么很显然最后答案也就是 ...
- 【BZOJ】1415 [Noi2005]聪聪和可可 期望DP+记忆化搜索
[题意]给定无向图,聪聪和可可各自位于一点,可可每单位时间随机向周围走一步或停留,聪聪每单位时间追两步(先走),问追到可可的期望时间.n<=1000. [算法]期望DP+记忆化搜索 [题解]首先 ...
- POJ 1088 DP=记忆化搜索
话说DP=记忆化搜索这句话真不是虚的. 面对这道题目,题意很简单,但是DP的时候,方向分为四个,这个时候用递推就好难写了,你很难得到当前状态的前一个真实状态,这个时候记忆化搜索就派上用场啦! 通过对四 ...
- 【bzoj5123】[Lydsy12月赛]线段树的匹配 树形dp+记忆化搜索
题目描述 求一棵 $[1,n]$ 的线段树的最大匹配数目与方案数. $n\le 10^{18}$ 题解 树形dp+记忆化搜索 设 $f[l][r]$ 表示根节点为 $[l,r]$ 的线段树,匹配选择根 ...
- [题解](树形dp/记忆化搜索)luogu_P1040_加分二叉树
树形dp/记忆化搜索 首先可以看出树形dp,因为第一个问题并不需要知道子树的样子, 然而第二个输出前序遍历,必须知道每个子树的根节点,需要在树形dp过程中记录,递归输出 那么如何求最大加分树——根据中 ...
- poj1664 dp记忆化搜索
http://poj.org/problem?id=1664 Description 把M个相同的苹果放在N个相同的盘子里,同意有的盘子空着不放,问共同拥有多少种不同的分法?(用K表示)5.1.1和1 ...
- 状压DP+记忆化搜索 UVA 1252 Twenty Questions
题目传送门 /* 题意:给出一系列的01字符串,问最少要问几个问题(列)能把它们区分出来 状态DP+记忆化搜索:dp[s1][s2]表示问题集合为s1.答案对错集合为s2时,还要问几次才能区分出来 若 ...
- ACM International Collegiate Programming Contest, Tishreen Collegiate Programming Contest (2017)- K. Poor Ramzi -dp+记忆化搜索
ACM International Collegiate Programming Contest, Tishreen Collegiate Programming Contest (2017)- K. ...
- zoj 3644(dp + 记忆化搜索)
题目链接:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=4834 思路:dp[i][j]表示当前节点在i,分数为j的路径条数,从 ...
随机推荐
- Learning ROS: Getting started with roswtf (检查ROS系统,找出问题)
本文主要部分来源于ROS官网的Tutorials. roswtf是ROS的检查工具,用于检查ROS安装和运行系统. Checking your installation&Offline mak ...
- 安全强化机制——SELinux
1.基本 SELINUX 安全性概念 SELINUX(Security Enhanced Linux),意思是安全增强型Linux, 是可保护你系统安全性的额外机制 在某种程度上 , 它可以被看作是与 ...
- 微信小程序开发正常,真机预览模式错误
原因是开发时设置了不检查域名是否合法,去后台设置上就可以了
- 大数据最后一公里——2021年五大开源数据可视化BI方案对比
个人非常喜欢这种说法,最后一公里不是说目标全部达成,而是把整个路程从头到尾走了一遍. 大数据在经过前几年的野蛮生长以后,开始与数据中台的概念一同向着更实际的方向落地.有人问,数据可视化是不是等同于数据 ...
- Sonarqube C# 配置资料整理
c#配置方式: http://www.cnblogs.com/CoderAyu/p/9416376.html http://www.cnblogs.com/jingridong/p/6513884.h ...
- CodeForce-792C Divide by Three(数学)
Divide by Three CodeForces - 792C 有一个正整数 n 写在黑板上.它有不超过 105 位. 你需要通过删除一些位使得他变成一个美丽的数,并且需要删除尽量少的位数.删除的 ...
- 最长公共前缀 js 实现代码
编写一个函数来查找字符串数组中的最长公共前缀: 输入 : ["abca","abc","abca","abc",&quo ...
- Jmeter扩展组件开发(2) - 扩展开发第一个demo的实现
maven工程src目录介绍 main:写代码 main/java:写Java代码 main/resources:写配置文件 test:写测试代码 test/java demo实现 创建Package ...
- python 建站教程
主端:安装nginx uwsgi django pymysql mysql安装mysql到http://repo.mysql.com/里面找 mysql57-community-release-el7 ...
- CF19E-Fairy【树形结构,差分】
正题 题目链接:https://www.luogu.com.cn/problem/CF19E 题目大意 给出\(n\)个点\(m\)条边的一张无向图,求有多少条边去掉后可以使得图变成一张二分图. \( ...