poj 3592 Instantaneous Transference tarjan缩点 最长路 ++tarjan模版 && spfa最长路模版

/*
题意：给定一个矩阵，西南角为起点，每个单元都有一定价值的金矿（#表示岩石，不可达，*表示时空门，可以到达指定单元）
      现在要求得最多可以获得多大利益
题解：强联通分量，最长路；如果没有时空门，就是纯粹的有向无环图的最长路了，现在出现时空门了，只要求强联通分量进行       缩点，对缩点后的图建立有向无环图，然后求最长路就ok了
*/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<climits>
#include<queue>
#include<vector>
using namespace std;
const int MAXV=49*49;
int n,m;

int value[MAXV];
char graph[41][41];

struct spfa//  spfa求最长路
{
	static const int INF =-(1<<30);
	struct Edge {
		int v,next;
		int w;
	} ep[MAXV*MAXV];
	int first[MAXV];
	int d[MAXV];
	int e;
	int V;
	inline void init(int n) 
	{
		memset(first,-1,sizeof(first));
		e=0;
		V=n;
	}
	void addEdge(int u, int v, int dd) {
		ep[e].v=v;
		ep[e].w=dd;
		ep[e].next=first[u];
		first[u]=e++;
	}
	void solve(int source) {
		static bool ok[MAXV];
		memset(ok,false,sizeof(ok));
		queue<int> q;
		q.push(source);
		for (int i = 0; i < V; i++) d[i] = INF;
		d[source] = 0;
		while(!q.empty())
		{
			int x=q.front();q.pop();
			ok[x]=false;
			for(int k=first[x];k!=-1;k=ep[k].next)
			{
				int v=ep[k].v;
				if(d[v]<d[x]+ep[k].w)
				{
					d[v]=d[x]+ep[k].w;
					if(ok[v]==false)
					{
						ok[v]=true;
						q.push(v);
					}
				}
			}
		}
	}
} g;
struct TarJan
{
	struct Edge
	{
	int v,next;
	}edge[4*MAXV];
	int first[MAXV],ins[MAXV],dfn[MAXV],low[MAXV],Stack[MAXV],scc[MAXV];
	int top,index,sccnum,e,V;
	void init (int n)
	{
	e=0;V=n;
	memset(first,-1,sizeof(first[0])*n);memset(ins,0,sizeof(ins[0])*n);
	memset(dfn,0,sizeof(dfn[0])*n);  memset(low,0,sizeof(dfn[0])*n);
	}
	void addEdge(int u,int v)
	{ edge[e].v=v;edge[e].next=first[u];first[u]=e++;}
	void solve()
	{
		index=1;top=-1;sccnum=-1;
		for(int i=0;i<V;i++)if(dfn[i]==0) tarjan(i);
	}
	void tarjan(int u)  
	{  
		int v;   low[u] = dfn[u] = index++;   Stack[++top]=u; ins[u] = true;  
		for (int k=first[u]; k!=-1; k=edge[k].next)  
		{  
			v = edge[k].v;  
			if (dfn[v] == 0) {  tarjan(v);  low[u]=min(low[u],low[v]); }  
			else if (ins[v]) {  low[u]=min(low[u],dfn[v]);  }  
		}  
		if (dfn[u] == low[u])  {  sccnum++;  
		do{v = Stack[top--];ins[v] = false; scc[v] = sccnum;}while(u != v); }  
	}  
}my;
void init()
{
	int x,y;
	memset(value,0,sizeof(value));
	my.init(n*m);
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<m;j++)
		if(graph[i][j]!='#')
		{
			if(i!=0)
			my.addEdge((i-1)*m+j,i*m+j);
			if(j!=0)
			my.addEdge(i*m+j-1,i*m+j);
			if(graph[i][j]=='*')
			{
				scanf("%d%d",&x,&y);
				if(graph[x][y]!='#')
				my.addEdge(i*m+j,x*m+y);
			}
		}
	}
	my.solve();
}
void solve()
{	
	g.init(my.sccnum+2);
	for(int i=0;i<n;i++)
	for(int j=0;j<m;j++)
	{
		if(graph[i][j]!='*'&&graph[i][j]!='#')
		{
			value[my.scc[i*m+j]]+=graph[i][j]-'0';
		}
	}
	
	for(int i=0;i<n*m;i++)
	{
		for(int k=my.first[i];k!=-1;k=my.edge[k].next)
		{
			int v=my.edge[k].v;
			if(my.scc[i]!=my.scc[v])
			{
				g.addEdge(my.scc[i],my.scc[v],value[my.scc[v]]);
			}
		}
	}
	int ans=0;
	g.solve(my.scc[0]);
	for(int i=0;i<=my.sccnum;i++)
	{
		if(g.d[i]>ans)
		ans=g.d[i];
	}
	printf("%d\n",ans+value[my.scc[0]]);
}
int main()
{
	int ca;
	scanf("%d",&ca);
	while(ca--)
	{
		scanf("%d%d",&n,&m);
		for(int i=0;i<n;i++)
		{
			scanf("%s",graph[i]);
		}
		init();
		solve();
	}
return 0;
}