PAT天梯赛练习题 L3-011. 直捣黄龙(多关键字SPFA+DFS)

L3-011. 直捣黄龙

时间限制
150 ms
内存限制
65536 kB
代码长度限制
8000 B
判题程序
Standard
作者
陈越

本题是一部战争大片 —— 你需要从己方大本营出发,一路攻城略地杀到敌方大本营。首先时间就是生命,所以你必须选择合适的路径,以最快的速度占领敌方大本营。当这样的路径不唯一时,要求选择可以沿途解放最多城镇的路径。若这样的路径也不唯一,则选择可以有效杀伤最多敌军的路径。

输入格式:

输入第一行给出2个正整数N(2 <= N <= 200,城镇总数)和K(城镇间道路条数),以及己方大本营和敌方大本营的代号。随后N-1行,每行给出除了己方大本营外的一个城镇的代号和驻守的敌军数量,其间以空格分隔。再后面有K行,每行按格式“城镇1 城镇2 距离”给出两个城镇之间道路的长度。这里设每个城镇(包括双方大本营)的代号是由3个大写英文字母组成的字符串。

输出格式:

按照题目要求找到最合适的进攻路径(题目保证速度最快、解放最多、杀伤最强的路径是唯一的),并在第一行按照格式“己方大本营->城镇1->...->敌方大本营”输出。第二行顺序输出最快进攻路径的条数、最短进攻距离、歼敌总数,其间以1个空格分隔,行首尾不得有多余空格。

输入样例:
10 12 PAT DBY
DBY 100
PTA 20
PDS 90
PMS 40
TAP 50
ATP 200
LNN 80
LAO 30
LON 70
PAT PTA 10
PAT PMS 10
PAT ATP 20
PAT LNN 10
LNN LAO 10
LAO LON 10
LON DBY 10
PMS TAP 10
TAP DBY 10
DBY PDS 10
PDS PTA 10
DBY ATP 10
输出样例:
PAT->PTA->PDS->DBY
3 30 210

 

 

注意一下当判断某个关键字条件时要保证优先级更高的关键字均已判断。其他就是比较麻烦的记录前驱,map记录标号和反向标号

代码:

#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<sstream>
#include<cstring>
#include<bitset>
#include<cstdio>
#include<string>
#include<deque>
#include<stack>
#include<cmath>
#include<queue>
#include<set>
#include<map>
using namespace std;
#define INF 0x3f3f3f3f
#define CLR(x,y) memset(x,y,sizeof(x))
#define LC(x) (x<<1)
#define RC(x) ((x<<1)+1)
#define MID(x,y) ((x+y)>>1)
typedef pair<int,int> pii;
typedef pair<string,int>psi;
typedef long long LL;
const double PI=acos(-1.0);

const int N=210;

int d[N],pre[N],lib[N];
int val[N],sumval[N],vis[N];
int cnt;
int n,k,s,t,z,va,roads;
vector<pii>E[N];
map<string,int>ref;
map<int,string>num;
string st,en,city,x,y;

void init()
{
	CLR(d,INF);
	CLR(val,0);
	cnt=0;
	for (int i=0; i<N; ++i)
		E[i].clear();
	ref.clear();
	num.clear();
	CLR(lib,0);
	CLR(sumval,0);
	CLR(pre,-1);	
	roads=0;
	CLR(vis,0);
}
void spfa(const int &s)
{
	int i;
	priority_queue<pii>Q;
	d[s]=0;
	lib[s]=0;
	sumval[s]=0;
	Q.push(pii(-d[s],s));
	while (!Q.empty())
	{
		int now=Q.top().second;
		Q.pop();
		int SZ=E[now].size();
		for (i=0; i<SZ; ++i)
		{
			int v=E[now][i].first;
			int w=E[now][i].second;
			if(d[v]>d[now]+w)
			{
				d[v]=d[now]+w;
				sumval[v]=sumval[now]+val[v];
				pre[v]=now;
				lib[v]=lib[now]+1;
				Q.push(pii(-d[v],v));
			}
			else if(d[v]==d[now]+w&&lib[v]<lib[now]+1)
			{
				sumval[v]=sumval[now]+val[v];		
				pre[v]=now;
				lib[v]=lib[now]+1;
			}
			else if(d[v]==d[now]+w&&lib[v]==lib[now]+1&&sumval[v]<sumval[now]+val[v])
			{
				sumval[v]=sumval[now]+val[v];		
				pre[v]=now;
				lib[v]=lib[now]+1;
			}
		}
	}
}
void dfs(int now)
{
	if(now==t)
	{
		++roads;
		return ;
	}
	for (int i=0;i<E[now].size(); ++i)
	{
		int v=E[now][i].first;
		int w=E[now][i].second;
		if(!vis[v]&&d[now]+w==d[v])
		{
			vis[v]=1;
			dfs(v);
			vis[v]=0;
		}
	}
}
int main(void)
{
	int i,j;
	while (cin>>n>>k>>st>>en)
	{
		init();
		if(ref.find(st)==ref.end())
		{
			ref.insert(psi(st,cnt));
			num[cnt]=st;
			++cnt;
		}
		for (i=0; i<n-1; ++i)
		{
			cin>>city>>va;
			ref[city]=cnt;
			val[cnt]=va;
			num[cnt]=city;
			++cnt;
		}
		for (i=0; i<k; ++i)
		{
			cin>>x>>y>>z;
			int u=ref[x],v=ref[y];
			E[u].push_back(pii(v,z));
			E[v].push_back(pii(u,z));
		}
		s=ref[st],t=ref[en];
		spfa(s);
		int flag=0;
		stack<string>pos;
		for (int tpre=t; ~tpre; tpre=pre[tpre],flag++)
			pos.push(num[tpre]);
		while (!pos.empty())
		{
			--flag;
			printf("%s%s",pos.top().c_str(),flag?"->":"\n");
			pos.pop();
		}
		vis[s]=1;
		dfs(s);
		vis[s]=0;
		printf("%d %d %d\n",roads,d[t],sumval[t]);
	}
	return 0;
}

转载于:https://www.cnblogs.com/Blackops/p/5783201.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是exp8-16.cpp的代码: ```c++ #include <iostream> #include <cstring> #include <algorithm> using namespace std; const int MAXN = 1005; const int MAXM = 20005; const int INF = 0x3f3f3f3f; int n, m, s, t; int head[MAXN], nxt[MAXM], to[MAXM], w[MAXM], c[MAXM], cnt; int dis[MAXN], cost[MAXN], pre[MAXN], vis[MAXN]; void add_edge(int u, int v, int d, int p) { nxt[++cnt] = head[u]; to[cnt] = v; w[cnt] = d; c[cnt] = p; head[u] = cnt; } bool spfa() { memset(dis, INF, sizeof(dis)); memset(cost, INF, sizeof(cost)); memset(vis, false, sizeof(vis)); memset(pre, -1, sizeof(pre)); dis[s] = 0; cost[s] = 0; vis[s] = true; while (true) { bool flag = false; for (int i = head[s]; i; i = nxt[i]) { int v = to[i]; if (dis[v] > dis[s] + w[i] && c[i] > 0) { dis[v] = dis[s] + w[i]; cost[v] = cost[s] + c[i]; pre[v] = i; if (!vis[v]) { vis[v] = true; flag = true; } } } for (int u = 1; u <= n; u++) { if (u != s && vis[u]) { vis[u] = false; for (int i = head[u]; i; i = nxt[i]) { int v = to[i]; if (dis[v] > dis[u] + w[i] && c[i] > 0) { dis[v] = dis[u] + w[i]; cost[v] = cost[u] + c[i]; pre[v] = i; if (!vis[v]) { vis[v] = true; flag = true; } } } } } if (!flag) { break; } } return pre[t] != -1; } int main() { cin >> n >> m >> s >> t; for (int i = 1; i <= m; i++) { int u, v, d, p; cin >> u >> v >> d >> p; add_edge(u, v, d, p); add_edge(v, u, d, p); } int ans1 = 0, ans2 = 0; while (spfa()) { int flow = INF; for (int i = t; i != s; i = to[pre[i] ^ 1]) { flow = min(flow, c[pre[i]]); } for (int i = t; i != s; i = to[pre[i] ^ 1]) { c[pre[i]] -= flow; c[pre[i] ^ 1] += flow; } ans1 += flow; ans2 += flow * cost[t]; } cout << ans1 << " " << ans2 << endl; return 0; } ``` 代码解释: - 第1行到第3行定义了常量MAXN、MAXM和INF,其中MAXN表示最大的点数,MAXM表示最大的边数,INF表示无穷大。 - 第5行到第7行定义了变量n、m、s和t,其中n表示点数,m表示边数,s表示起点,t表示终点。 - 第9行到第19行定义了邻接表存图所需要的变量和函数。变量head数组表示每个点的第一条边,nxt数组表示每条边的下一条边,to数组表示每条边的另一个端点,w数组表示每条边的长度,c数组表示每条边的花费,cnt表示边的数量。函数add_edge用来添加一条从u到v长度为d、花费为p的边。 - 第21行到第50行是SPFA算法的实现。变量dis数组表示从起点s到每个点的最短路径长度,cost数组表示从起点s到每个点的最小花费,pre数组表示每个点在最短路径上的前驱边,vis数组表示每个点是否在队列中。在SPFA算法中,我们先将dis、cost、pre和vis数组初始化,然后将起点s加入队列中。接下来,我们不断从队列中取出元素,更新其可达的点的dis、cost、pre和vis数组。具体来说,我们对于每个点u,遍历其所有相邻的点v,如果dis[v]>dis[u]+w[i]且c[i]>0,则更新dis[v]、cost[v]、pre[v]和vis[v],并将v加入队列中。然后,我们再次遍历所有点u,如果u在队列中且vis[u]为true,则遍历它所有相邻的点v,进行同样的更新操作。如果队列为空,则SPFA算法结束。最后,我们返回pre[t]!=-1判断是否存在从起点s到终点t的路径。 - 第52行到第56行是主函数。我们先读入n、m、s和t,并使用add_edge函数添加所有边。然后,我们使用最小费用最大流算法求解最短路径问题。具体来说,我们不断使用spfa函数找到从起点s到终点t的最短路径,然后在这条路径上找到最小剩余容量flow,将这个flow加入最大流答案ans1中,将flow乘以路径总花费加入最小费用ans2中,然后更新每条边的剩余容量和反向边的容量。如果不存在从起点s到终点t的路径,则最小费用最大流算法结束。最后,我们输出ans1和ans2即可。 注意:本程序中使用了SPFA算法和最小费用最大流算法,需要读者具备相关知识。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值