习题：天降馅饼（tarjan缩点）

题目描述

今天是K公司成立百年的日子，在这个喜庆的日子里，K公司决定给所有到他们线下门店的顾客发红包。林老师得知这个消息后，立马出门抢红包去啦。在林老师居住的P市内，一共有n家K公司的门店，在n家门店之间，有m条单向连通的道路，每家门店所发放的红包金额大小不同，且林老师在同一家门店只能领取一次红包（即第二次到达同一家门店是没有红包领取的）。林老师从离家最近的门店s开始出发抢红包，一路沿着单向道路行走，可K公司的门店不光发红包，还有各种各样其他的优惠活动吸引着林老师，于是，林老师在出发到第一家门店的路上四处打听，了解到有q家门店有他想要的东西，他希望在这q家门店中的一家停下来购物。林老师想知道他最多能得到多少的红包并能够购买他心仪的商品（即在给出的q家门店中的一家停下）。

输入格式

第一行输入两个数n和m(1≤n,m≤500000)。 接下来m行，每行输入两个数字u和v(1≤u,v≤n)，表示从u到v有一条单向道路。 接下来n行，每行输入一个整数x(0≤x≤4000)，按顺序表示标号为1到n的门店发放红包的数额。 接下来输入两个数s，q，表示离林老师最近的门店的编号和有林老师心仪商品的门店数。 接下来输入q个数，表示这q个门店的编号。 输入保证至少有一条道路能从离林老师家最近的门店到有林老师心仪商品的门店。

输出格式

输出一个整数，表示林老师最多能得到的红包金额数。

样例

样例输入1

6 7
1 2
2 3
3 5
2 4
4 1
2 6
6 5
10
12
8
16
1
5
1 4
4 3 5 6

样例输出1

47

 

思路：可能有些人看到这道题第一个想到的是dij最短路算法，但是很可惜，这道并不满足dij的性质，同时如果对边进行取反，那更是错的，因为dij不能处理负权边。

接着便想到了广搜，但是如果对图直接进行广搜，时间一定会炸

考虑对广搜进行优化，因为路是可以重复走的，况且每个点只能领取一次，便想到了tarjan缩点，

喜闻(van)乐见的代码时间

#include<iostream>
#include<cstring>
#include<climits>
#include<cstdio>
#include<vector>
#include<stack>
#include<queue>
using namespace std;
void read(int &x)//快速读入
{
	x=0;
	int f=1;
	char c=getchar();
	while('0'>c||c>'9')
	{
		if(c=='-')
			f=-1;
		c=getchar();
	}
	while('0'<=c&&c<='9')
	{
		x=(x<<3)+(x<<1)+c-'0';
		c=getchar();
	}
	x*=f;
}
int n,m;
int s,q;
int tot;
int scc;
int ans=INT_MIN;
int money[500005];
int dfs_num[500005];
int low[500005];
int bel[500005];
int x[500005];
int want[500005];
int lov[500005];
int dis[500005];
bool vis[500005];
bool ins[500005];
stack<int> tar_s;
queue<int> que;
vector<int> g[500005];//原图
vector<int> g2[500005];//缩点之后的图
void tar(int u)//求每个点属于哪一个强连通分量
{
	dfs_num[u]=low[u]=++tot;
	tar_s.push(u);
	ins[u]=1;
	for(int i=0;i<g[u].size();i++)
	{
		int v=g[u][i];
		if(dfs_num[v]==0)
		{
			tar(v);
			low[u]=min(low[u],low[v]);
		}
		else if(ins[v]==1)
		{
			low[u]=min(low[u],low[v]);
		}
	}
	if(dfs_num[u]==low[u])
	{
		int v;
		scc++;
		do
		{
			v=tar_s.top();
			tar_s.pop();
			bel[v]=scc;
			ins[v]=0;
		}while(v!=u);
	}
}
int main()
{
	read(n);
	read(m);
	for(int i=1;i<=m;i++)
	{
		int s,e;
		read(s);
		read(e);
		g[s].push_back(e);
	}
	for(int i=1;i<=n;i++)
		read(x[i]);
	read(s);
	read(q);
	for(int i=1;i<=q;i++)
	{
		int t;
		read(t);
		want[t]=1;
	}
	for(int i=1;i<=n;i++)
		if(dfs_num[i]==0)
			tar(i);
	for(int i=1;i<=n;i++)
	{
		if(want[i]==1)//将原来的最后到的点在新图上进行映射
		{
			lov[bel[i]]=1;
		}
		money[bel[i]]+=x[i];
		for(int j=0;j<g[i].size();j++)
		{
			if(bel[i]!=bel[g[i][j]])
			{
				g2[bel[i]].push_back(bel[g[i][j]]);
			}
		}
	}
	s=bel[s];
	dis[s]=money[s];
	que.push(s);
	while(!que.empty())
	{
		int t1=que.front();
		que.pop();
		vis[t1]=1;
		for(int i=0;i<g2[t1].size();i++)
		{
			int v=g2[t1][i];
			if(dis[v]<dis[t1]+money[v])
			{
				dis[v]=dis[t1]+money[v];
				if(lov[v]==1)
					ans=max(ans,dis[v]);
				que.push(v);
			}
		}
	}
	cout<<ans;
	return 0;
}