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的路径条数,从 ...
随机推荐
- springcloud <zuul2.0静态配置>
server: port: 9006 spring: application: name: cloud-zuul-wangbiao # zipkin: # base-url: http://local ...
- 回调与Promise
Promise 对象就是用于表示一个异步操作的最终状态(成功或失败).它的流程就是在什么状态下需要执行什么样的操作. resolve简单理解就是一步操作执行成功后的回调函数 then是Promise对 ...
- Qt5之反射机制（内省）
Qt的元对象系统除了提供信号/槽机制的特性之外,它还提供了以下特性: QObject::metaObject() 返回关联的元对象 QMetaObject::className() 在运行时状态下返回 ...
- vue 手机号码隐藏中间
template {{ phone | mobileFilter }} data export default { filters:{ mobileFilter(val){ ...
- Java的参数传递是值传递还是引用传递？
一.前言 首先先说结论,Java中方法参数传递方式是按值传递.如果参数是基本类型,传递的是基本类型的字面量值的拷贝.如果参数是引用类型,传递的是该参量所引用的对象在堆中地址值的拷贝. 接下来深入了解一 ...
- 依赖注入Bean属性——手动装配Bean
一.构造方法注入 其中,可以根据不同的参数列表调用不同的重载的构造方法: 其中,基本数据类型没有包,引用类型都有包路径,基本类型对应封装类: 二.通过property标签调用类的set方法注入 三.通 ...
- python库--pandas--Series
方法 返回数据类型 参数 说明 Series(一维) .Series() Series 实例s 创建一维数据类型Series data=None 要转化为Series的数据(也可用dict ...
- 这些解决 Bug 的套路，你都会了不？
最近整理了我原创的 140 篇编程经验和技术文章,欢迎大家阅读,一起成长!指路:https://t.1yb.co/ARnD 大家好,我是鱼皮. 学编程的过程中,我们会遇到各式各样的 Bug,也常常因为 ...
- 学习Tomcat（二）之容器概览
Tomcat容器的Server模块有管理容器的启动和关闭.管理了容器内的服务组件Service.管理了全局JNDI资源的功能,对Tomcat容器的生命周期管理有重要意义.Tomcat的服务组件则是To ...
- symfony生成路由
控制器里生成地址 $this->generateUrl('course_manage_show_test', array('id' => 1)) twig前端文件生成地址: {{ path ...