poj 3592 Instantaneous Transference

http://poj.org/problem?id=3592

这道题先有tarjan算法缩点，然后再建图，再有spfa算法求最长路。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#define maxn 300000
using namespace std;

const int inf=-1<<30;
int e,head[maxn],belong[maxn],stack1[maxn],top,dfn[maxn],low[maxn],bcc_clock,bcnt,n,m,num1[100][100],num,gg[maxn],point[maxn],cc[maxn],ee,head1[maxn],dis[maxn],cnt[maxn],temp;
bool vis[maxn],visi[maxn];
int dir[2][2]= {{0,1},{1,0}};
struct node
{
    int u,v,next;
} p[maxn];
struct node1
{
    int u,v,w,next;
} pp[maxn];
char g[50][50];

void add(int u,int v)
{
    p[e].u=u;
    p[e].v=v;
    p[e].next=head[u];
    head[u]=e++;
}

void addnode(int u,int v,int w)
{
    pp[ee].v=v;
    pp[ee].u=u;
    pp[ee].w=w;
    pp[ee].next=head1[u];
    head1[u]=ee++;
}

void tarjan(int u)
{
    vis[u]=true;
    dfn[u]=low[u]=++bcc_clock;
    stack1[++top]=u;
    for(int i=head[u]; i!=-1; i=p[i].next)
    {
        int v=p[i].v;
        if(!dfn[v])
        {
            tarjan(v);
            low[u]=min(low[u],low[v]);
        }
        else if(vis[v])
            low[u]=min(low[u],dfn[v]);
    }
    if(dfn[u]==low[u])
    {
        bcnt++;
        int j;
        do
        {
            j=stack1[top--];
            vis[j]=false;
            belong[j]=bcnt;
        }
        while(j!=u);
    }
}

void init()
{

    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    memset(belong,0,sizeof(belong));
    memset(cc,0,sizeof(cc));
    memset(vis,false,sizeof(vis));
}

void del()
{
    init();
    for(int i=0; i<n*m; i++)
    {
        if(!dfn[i])
        {
            tarjan(i);
        }
    }
}

bool ralex(int u,int v,int c)
{
    if(dis[v]<dis[u]+c)
    {
        dis[v]=dis[u]+c;
        return true;
    }
    return false;
}



bool spfa(int src)
{
    memset(visi,false,sizeof(visi));
    memset(cnt,0,sizeof(cnt));
    visi[src]=true;
    for(int i=0; i<=bcnt; i++)
    {
        dis[i]=inf;
    }
    queue<int>q;
    q.push(src);
    dis[src]=0;
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        visi[u]=false;
        for(int i=head1[u]; i!=-1; i=pp[i].next)
        {
            if(ralex(u,pp[i].v,pp[i].w)&&!visi[pp[i].v])
            {
                if((++cnt[pp[i].v])>n*m) return false;
                q.push(pp[i].v);
                visi[pp[i].v]=true;
            }
        }
    }
    temp=0;
    for(int i=1; i<bcnt+1; i++)
    {
        temp=max(temp,dis[i]);
    }
    return true;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&m);
        num=0;
        int cn=0;
        memset(head1,-1,sizeof(head1));
        memset(head,-1,sizeof(head));
        e=0,top=0,bcnt=0,bcc_clock=0,ee=0;
        getchar();
        for(int i=0; i<n; i++)
        {
            scanf("%s",g[i]);
        }
        for(int i=0; i<n; i++)
        {
            for(int j=0; j<m; j++)
            {
                int k=i*m+j;
                if(g[i][j]=='#')
                {
                    gg[k]=-1;
                    continue;
                }
                else
                {
                    if(g[i][j]=='*')
                    {
                        point[cn++]=k;
                        gg[k]=0;
                    }
                    else if(g[i][j]>='0'&&g[i][j]<='9')
                    {
                        gg[k]=g[i][j]-'0';
                    }
                    for(int c=0; c<2; c++)
                    {
                        int xx=i+dir[c][0];
                        int yy=j+dir[c][1];
                        if(xx<n&&yy<m)
                        {
                            if(g[xx][yy]!='#')
                            {
                                add(k,xx*m+yy);
                            }
                        }
                    }

                }
            }
        }
        for(int i=0; i<cn; i++)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            if(g[x][y]!='#')
            {
                add(point[i],x*m+y);
            }
        }
        del();
        for(int i=0; i<n*m; i++)
        {
            cc[belong[i]]+=gg[i];
        }
        addnode(0,belong[0],cc[belong[0]]);
        for(int i=0; i<n*m; i++)
        {
            for(int j=head[i]; j!=-1; j=p[j].next)
            {
                int v=p[j].v;
                if(belong[i]!=belong[v])
                {
                    addnode(belong[i],belong[v],cc[belong[v]]);
                }
            }
        }
        spfa(0);
        printf("%d\n",temp);
    }
    return 0;
}