PAT1003 三种方法附注释

1.朴素的点权做法（邻接表）

#include<bits/stdc++.h>
using std::cin;
using std::cout;
using std::endl;
using std::vector;
using std::fill;
struct node
{
	int pointer;
	int distance;
	node(){}
	node(int _pointer,int _distance):pointer(_pointer),distance(_distance) {}//构造函数
	node(int _pointer):pointer(_pointer) {}//构造函数
};
vector<node> adj[505];//存储图
bool visit[505];//是否被访问过
int dis[505];//到第i个点最短路径
int rescue_num[505];//存储到第i个点的营救人数的最大值
int city_rescue[505];//存储每个点（城市）的点权（即题目已经给出的每个城市的营救人员数）
int route_num[505];//存储到第i个点最短路径数目
int n,m,sta,fin;//n个城市，m条边，sta起点，fin终点
void dj(int root)
{
	fill(dis,dis+n+1,0x3fffffff);
	dis[root]=0;
	route_num[root]=1;
	rescue_num[root]=city_rescue[root];//初始化三个数组，1.dis[root]=0，其余设为极大值，2.route_num[root]=1,其余设为0，3.rescue_num[root]=city_rescue[root]，其余设置为0
	for(int i=0;i<n;i++)
	{
		int u=-1;
		int minn=0x3fffffff;
		for(int j=0;j<n;j++)
		{
			if(visit[j]==false&&dis[j]<minn)
			{
				u=j;
				minn=dis[j];
			}
		}
		if(u==-1) return;
		visit[u]=true;
		for(int i=0;i<adj[u].size();i++)
		{
			int to=adj[u][i].pointer;
			if(!visit[to]&&dis[to]>dis[u]+adj[u][i].distance)
			{
				dis[to]=dis[u]+adj[u][i].distance;
				route_num[to]=route_num[u];//路径数目继承至上一个点
				rescue_num[to]=rescue_num[u]+city_rescue[to];//总营救人员加上当前这个点（to）的人数
			}
			else if(dis[to]==dis[u]+adj[u][i].distance)//如果最短路径相等
			{
				if(rescue_num[to]<rescue_num[u]+city_rescue[to]) rescue_num[to]=rescue_num[u]+city_rescue[to];//考虑比较两条路径谁的营救人数更多
				route_num[to]+=route_num[u];//叠加上当前的路径数目
			}
		}
	}
}
void windows_cmd_support_utf8(void)
{
	#ifdef WIN32
	    system("chcp 65001 & cls");
	#endif
}
signed main(void)
{
	windows_cmd_support_utf8();//解决dev c++6.2不能输出中文的问题
	int c1,c2,l;
	scanf("%d%d%d%d",&n,&m,&sta,&fin);
	for(int i=0;i<n;i++) scanf("%d",&city_rescue[i]);//注意题目要求的顶点是从0~n-1还是1~n
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d%d",&c1,&c2,&l);//注意是双向建边
		adj[c1].push_back(node (c2,l));
		adj[c2].push_back(node (c1,l));
	}
	dj(fin);//也可以是dj(sta)，这样的话下边数组下标改成fin即可
	printf("%d %d",route_num[sta],rescue_num[sta]);
	return 0;
}

2.堆优化的点权做法

#include<bits/stdc++.h>
using std::cin;
using std::cout;
using std::endl;
using std::string;
using std::fill;
int num[1002],weight[1002],value[1002],dis[1002];
bool visit[1002];
int n,m,sta,fin;
struct node
{
	int distance;
	int point;
	node(int _distance,int _point):distance(_distance),point(_point) {}
	node() {}
	node(int _point):point(_point) {}
	friend bool operator < (const node &a,const node &b)//重载运算符<，使得优先队列从小到大排序
	{
		return a.distance<b.distance;
	}
};
std::priority_queue<node> q;
std::vector<node> adj[1002];
void dj(int root)
{
	std::fill(dis,dis+n+1,0x3fffffff);
	num[root]=1;
	dis[root]=0;
	value[root]=weight[root];//初始化，和朴素的一样
	q.push(node(0,root));
	while(!q.empty())
	{
		int to=q.top().point;
		q.pop();
		if(visit[to]) continue;
		visit[to]=true;
		for(int i=0;i<adj[to].size();i++)
		{		
			int v=adj[to][i].point;
			int w=adj[to][i].distance;
			if(dis[v]>dis[to]+w)
			{
				dis[v]=dis[to]+w;
				num[v]=num[to];//继承最短路径数目
				value[v]=value[to]+weight[v];//更新点权
				q.push(node(-dis[v],v));//保证其在优先队列中的优先级
			}
			else if(dis[v]==dis[to]+w)
			{
				if(value[v]<value[to]+weight[v]) value[v]=value[to]+weight[v];//如果最短距离相同，比较能否找到更大的点权
				num[v]+=num[to];//叠加最短路径数目，注意不能写在上面那个if里面去
			}
		}		
	}
}
signed main(void)
{
	int c1,c2,l;
	scanf("%d%d%d%d",&n,&m,&sta,&fin);
	for(int i=0;i<n;i++) scanf("%d",&weight[i]);
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d%d",&c1,&c2,&l);
		node temp1(l,c2);
		node temp2(l,c1);
		adj[c1].push_back(temp1);
		adj[c2].push_back(temp2);
	}
	dj(sta);
	printf("%d %d",num[fin],value[fin]);
	return 0;
}

3.dfs+迪杰斯特拉

#include<bits/stdc++.h>
using std::cin;
using std::cout;
using std::endl;
using std::string;
using std::fill;
using std::vector;
int n,m,sta,fin;
int optvalue;//第二标尺的最优值
int dis[505],city_rescue[505];
int edge[501][501];
bool visit[501];
int route_num[501];
struct node
{
	int point;
	int distance;
	node(int _point,int _distance):point(_point),distance(_distance) {}
	node() {}
};
vector<node> adj[505];
vector<int> pre[505];//存放结点的前驱结点，如果要查询一个结点是否在某条路径中，可以使用set<int>
vector<int> path,temppath;//最优路径与暂时路径(存储在temppath中的结点是逆序的)
void dj(int root)
{
	fill(dis,dis+n+1,0x3fffffff);
	dis[root]=0;
	route_num[root]=1;
	for(int i=0;i<n;i++)
	{
		int u=-1,minn=0x3fffffff;
		for(int j=0;j<n;j++)
		{
			if(visit[j]==false&&dis[j]<minn)
			{
				u=j;
				minn=dis[j];
			}
		}
		if(u==-1) return;
		visit[u]=true;
		for(int v=0;v<adj[u].size();v++)
		{
			int to=adj[u][v].point;
			int w=adj[u][v].distance;
			if(!visit[to]&&dis[to]>dis[u]+w)
			{
				dis[to]=dis[u]+w;
				pre[to].clear();//1.
				pre[to].push_back(u);//2.实际上就多了三步
				route_num[to]=route_num[u];
			}
			else if(dis[to]==dis[u]+adj[u][v].distance)
			{
				pre[to].push_back(u);//3.
				route_num[to]+=route_num[u];
			}
		}
	}
}
int calculate1()//计算边权
{
	int value=0;
	for(int i=temppath.size()-1;i>0;i--)//边的数量是点的数量-1，因此不包含i=0
	{
		int id=temppath[i];
		int id_next=temppath[i-1];
		value+=edge[id][id_next];
	}
	return value;
}
int calculate2()//计算点权
{
	int value=0;
	for(int i=temppath.size()-1;i>=0;i--)
	{
		int id=temppath[i];
		value+=city_rescue[id];
	}
	return value;
}
void dfs(int v)
{
	if(v==sta)//递归边界
	{
		temppath.push_back(v);
		int value=calculate2();//如果是边权的话，看calculate1
		if(value>optvalue)//if中的条件是value优于当前optvalue；如果是边劝最小的话这里就应该是value<optvalue
		{
			optvalue=value;//更新最优值
			path=temppath;//更新最优路径
		}
		temppath.pop_back();
		return;
	}
	temppath.push_back(v);
	for(int i=0;i<pre[v].size();i++)
	{
		dfs(pre[v][i]);//递归
	}
	temppath.pop_back();
}
void windows_cmd_support_utf8(void)
{
	#ifdef WIN32
	    system("chcp 65001 & cls");
	#endif
}
signed main(void)
{
	windows_cmd_support_utf8();
	int c1,c2,l;
	scanf("%d%d%d%d",&n,&m,&sta,&fin);
	for(int i=0;i<n;i++) scanf("%d",&city_rescue[i]);
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d%d",&c1,&c2,&l);
		node temp1(c2,l);
		node temp2(c1,l);
		adj[c1].push_back(temp1);
		adj[c2].push_back(temp2);
	}
	dj(sta);
	dfs(fin);
	printf("%d %d\n",route_num[fin],optvalue);
	return 0;
}