[CSP-S模拟测试]:旅行计划(分块+DP)

题目传送门(内部题83)


输入格式

  第一行两个整数$n,m$
  接下来$m$行,每行三个整数,$u,v,w$,表示从$u$到$v$有一条权值为$w$的边
  接下来一行有一个整数$q$,表示$q$天
  接下来$q$行,每行三个整数,$s_i,t_i,k_i$,表示从$s_i$到$t_i$至少经过$k_i$条边


输出格式

  一共$q$行,每行一个整数,表示在第$i$天中的最短距离,如果没有符合要求的答案,输出$-1$。


样例

样例输入:

3 3
1 2 1
2 3 10
3 1 100
3
1 1 1
1 2 1
1 3 1

样例输出:

111
1
11


数据范围与提示

  对于$30\%$的数据,$k_ileqslant 50$
  对于另外$30\%$的数据,$qleqslant 50$
  对于$100\%$的数据,$nleqslant 50,mleqslant 10,000,wileqslant 10,000,qleqslant 100,000,kileqslant 10,000$


题解

$DP$和分块可能很难接合在一起,但是他们做到了。

不妨设$f[i][j][k]$表示从$i$到$j$恰好走$k$步的方案数,稍做修改就能做到至少走$k$步的方案数。

再设$p[i][j]$表示从$i$到$j$恰好走$100$步的方案数,这个数组也就是这道题的关键,将其分块了。

最后设$g[i][j][k]$表示从$i$到$j$恰好走$100 imes k$的方案数即可。

最后答案就是$minlimits_{i=1}^n(g[s][i][k/100]+f[i][t][k\%100])$。

时间复杂度:$Theta(100 imes n+q imes n)$。

期望得分:$100$分。

实际得分:$100$分。


代码时刻

#include<bits/stdc++.h>
using namespace std;
int n,m,q;
int d[51][51],p[51][51],g[51][51][200],f[51][51][200];
int main()
{
	memset(d,0x3f,sizeof(d));
	memset(g,0x3f,sizeof(g));
	memset(f,0x3f,sizeof(f));
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++)
	{
		int u,v,w;
		scanf("%d%d%d",&u,&v,&w);
		d[u][v]=min(d[u][v],w);
	}
	for(int i=1;i<=n;i++)
	{
		f[i][i][0]=0;
		for(int j=0;j<n+100;j++)
			for(int k=1;k<=n;k++)
				for(int l=1;l<=n;l++)
					f[i][l][j+1]=min(f[i][l][j+1],f[i][k][j]+d[k][l]);
	}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
		{
			p[i][j]=f[i][j][100];
			for(int k=n+100;~k;k--)
				f[i][j][k]=min(f[i][j][k],f[i][j][k+1]);
		}
	for(int i=1;i<=n;i++)
	{
		g[i][i][0]=0;
		for(int j=0;j<100;j++)
			for(int k=1;k<=n;k++)
				for(int l=1;l<=n;l++)
					g[i][l][j+1]=min(g[i][l][j+1],g[i][k][j]+p[k][l]);
	}
	scanf("%d",&q);
	while(q--)
	{
		int s,t,k,ans=0x3f3f3f3f;
		scanf("%d%d%d",&s,&t,&k);
		for(int i=1;i<=n;i++)ans=min(ans,f[i][t][k%100]+g[s][i][k/100]);
		if(ans==0x3f3f3f3f)puts("-1");
		else printf("%d
",ans);
	}
	return 0;
}

rp++

原文地址:https://www.cnblogs.com/wzc521/p/11730104.html