POJ 3592 Instantaneous Transference

题目大意：

        给出一个n*m的格子地图，每一格上面是0~9，“*”或“#”。如果格子上是数字代表这个格子上有当前数量的矿石。如果是“*” 代表着当前格子是一个传送阵可以传送到指定的地方。如果是“#”代表当前格子不可达。

        现在有一个矿车在坐标（0，0），也就是左上角。他只能向右和下行驶。当遇到传送阵时可以被传送到指定的位置。当他遇到数字时就可以得到那些数量的矿石，那个地方的矿石数量就变为“0”。问矿车最多可以采多少矿。

解题思路：

1、我们先要根据格子地图建图。（注1）

2、因为建立出的图可能会有环，我们要用Tarjan算法将环缩成点。（注2）

3、我们将缩点处理后的图重新建图。（注3）

4、对图求最长路。

注意：

1、要注意到“#”和从“#”出发的边是不建立的。当从“*”出发时需要建三条边：右，下，指定位置。

2、根据题意，经过的地方将不再有矿石。所以将环缩成点是最好的选择。因为如果可以到环中的任意一点，其他点也就可以到了。

3、重新建图的目的是为后来更容易的求最长路做准备。建图时对于每一条边（u,v），边的权值是点v所在的缩点的总权值。

4、求完最长路后得出的数值一定要再加上点0所在的缩点的总权值。

5、注意不要手残（手残毁一天啊~~~ T _ T ）.

下面是代码：

#include <stdio.h>#include <string.h>#include <queue>using namespace std;const int MAXN = 2005;struct node{    int u,to,next;} edge[MAXN*10];struct node1{    int v,w,next;} newedge[MAXN*10];int n,m,p;int dfn[MAXN],low[MAXN],vis[MAXN],head[MAXN],stack1[MAXN],top,cnt,time;int num[MAXN],numcnt,numw[MAXN],map1[MAXN];int newhead[MAXN],newcnt;int dis[MAXN];char s[45][45];void init(){    memset(dfn,0,sizeof(dfn));    memset(low,0,sizeof(low));    memset(vis,0,sizeof(vis));    memset(head,-1,sizeof(head));    memset(newhead,-1,sizeof(newhead));    memset(numw,0,sizeof(numw));    memset(stack1,-1,sizeof(stack1));    memset(num,-1,sizeof(num));    top=0;    cnt=0;    newcnt=0;    time=1;    numcnt=0;}int min(int a,int b){    if(a>b)a=b;    return a;}void addedge(int u,int v){    edge[cnt].u=u;    edge[cnt].to=v;    edge[cnt].next=head[u];    head[u]=cnt;    cnt++;}void Buildedge(){    int x,y;    for(int i=0; i<n; i++)    {        scanf("%s",s[i]);    }    for(int i=0; i<n; i++)    {        for(int j=0; j<m; j++)        {            if(s[i][j]=='*')            {                scanf("%d%d",&x,&y);                if(s[x][y]!='#')                    addedge(i*m+j,x*m+y);                map1[i*m+j]=0;                if(i+1<n&&s[i+1][j]!='#')                {                    addedge(i*m+j,(i+1)*m+j);                }                if(j+1<m&&s[i][j+1]!='#')                {                    addedge(i*m+j,i*m+j+1);                }            }            else if(s[i][j]!='#')            {                map1[i*m+j]=s[i][j]-'0';                if(i+1<n&&s[i+1][j]!='#')                {                    addedge(i*m+j,(i+1)*m+j);                }                if(j+1<m&&s[i][j+1]!='#')                {                    addedge(i*m+j,i*m+j+1);                }            }            else            {                map1[i*m+j]=0;            }        }    }}void dfs(int u,int fa){    dfn[u]=time;    low[u]=time;    time++;    vis[u]=1;    stack1[top]=u;    top++;    for(int i=head[u]; i!=-1; i=edge[i].next)    {        int v=edge[i].to;        if(!vis[v])        {            dfs(v,u);            low[u]=min(low[u],low[v]);        }        else if(vis[v]==1)        {            low[u]=min(low[u],dfn[v]);        }    }    if(low[u]==dfn[u])    {        while(stack1[top]!=u&&top>0)        {            top--;            num[stack1[top]]=numcnt;            vis[stack1[top]]=2;            numw[numcnt]+=map1[stack1[top]];        }        numcnt++;    }}void addedgenew(int u,int v){    for(int i=newhead[u]; i!=-1; i=newedge[i].next)    {        if(newedge[i].v==v)return;    }    newedge[newcnt].v=v;    newedge[newcnt].w=numw[v];    newedge[newcnt].next=newhead[u];    newhead[u]=newcnt;    newcnt++;}int spfa(int src){    queue <int> q;    q.push(src);    memset(vis,0,sizeof(vis));    vis[src]=true;    for(int i=0; i<=numcnt; i++)    {        dis[i]=-1;    }    dis[src]=0;    while(!q.empty())    {        int p,t=q.front();        q.pop();        p=newhead[t];        vis[t]=false;        while(p!=-1)        {            if(dis[newedge[p].v]<dis[t]+newedge[p].w)            {                dis[newedge[p].v]=dis[t]+newedge[p].w;                if(!vis[newedge[p].v])                {                    vis[newedge[p].v]=true;                    q.push(newedge[p].v);                }            }            p=newedge[p].next;        }    }    int max1=0;    for(int i=0; i<numcnt; i++)    {        if(dis[i]>max1)        {            max1=dis[i];        }    }    return max1+numw[src];}void Tarjan(){    for(int i=0; i<p; i++)    {        if(!vis[i])dfs(i,-1);    }}void BuildAgain(){    for(int i=0; i<p; i++)    {        for(int j=head[i]; j!=-1; j=edge[j].next)        {            if(num[i]!=num[edge[j].to]&&num[i]!=-1&&num[edge[j].to]!=-1)            {                addedgenew(num[i],num[edge[j].to]);            }        }    }}int main(){    int t;    scanf("%d",&t);    while(t--)    {        init();        scanf("%d%d",&n,&m);        Buildedge();//建图        p=n*m; //总点数        Tarjan();//求强连通分量用来缩点        BuildAgain();//重新建图        printf("%d\n", spfa(num[0]));//求最长路    }    return 0;}