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的路径条数,从 ...
随机推荐
- go语言 切片表达式
切片表达式 切片的底层就是一个数组,所以我们可以基于数组通过切片表达式得到切片. 切片表达式中的low和high表示一个索引范围(左包含,右不包含),得到的切片长度=high-low,容量等于得到的切 ...
- for循环操作(for...in、forEach)
1.for...in语句用于对数组或者对象的属性进行循环操作,是for循环的一种. 注意:该方法可用于数组或对象. 语法: for(变量 in 对象/数组){} 如: var obj = { nam ...
- 「萌新指南」SOA vs. 微服务:What’s the Difference?
实话实说,在我还没有实习之前,我是连 SOA 是啥都不知道的,只听说过微服务,毕竟微服务实在太火了,想不知道都难,我觉得实习的时候肯定也是微服务,进组之后发现是 SOA 架构,当时都懵了,看了很多文档 ...
- Shell中的运算
1.运算方式及运算符号 2.SHELL 中常用的运算命令 3.相关操作演示 1.用脚本写一个10秒倒计时 脚本的执行: 2.编写脚本,1分10秒的倒计时 执行脚本: 3.编写脚本,制作一个计算器 脚本 ...
- golang error错误处理
error定义 数据结构 go语言error是一普通的值,实现方式为简单一个接口. // The error built-in interface type is the conventional i ...
- python 文件批量改名重命名 rename
path = '/Volumes/Seagate/dev/imgs/' os.chdir(path) print('cwd: ', os.getcwd()) for f in os.listdir(' ...
- vue 输入框禁止输入空格 ,只能输入数字,禁止输入数字
正则表达式: @input="form.userName = form.userName.replace(/\s+/g,'')" ( 禁止输入空格) @input=&q ...
- Python - 通过PyYaml库操作YAML文件
PyYaml简单介绍 Python的PyYAML模块是Python的YAML解析器和生成器 它有个版本分水岭,就是5.1 读取YAML5.1之前的读取方法 def read_yaml(self, pa ...
- markdown的骚气操作(一)
markdown 系列其他内容 markdown的骚气操作(一)✓ latex的骚气操作(二) 本文目标 主要介绍markdown锚点.索引脚注.对勾及选择框.表格显示位置和符号显示位置.绘制 ...
- go的database/sql库中db.Exce()
db.Exec(query string, args ...interface{}) Db.Exec(`CREATE TABLE IF NOT EXISTS STU(ID int(8) PRIMARY ...