poj 3592 强连通+缩点+dfs

传送们

题意：给你一个m*n的矩阵（抱歉，由于习惯，我这里mn和题中的mn是反的，我的m是行，n是列），只可以向下或向右走，如果是数字，代表矿石数，*代表传送点，#代表不能到的地方。有多少个*之后就有多少个坐标，代表被传送到的地方，传送是主动选择的，可以选择不传送。问你最多能采多少矿石。

思路：根据矩阵建立有向图，#的地方就不用管啦，*的地方多建一条边就是了。然后根据有向图求强连通分量，然后把每个分量当作单独的点来处理，它的矿石数等于它包含的点的矿石数之和，再根据原来点的边给新点建边，注意去重，最后我是从原来0 0坐标所在的强连通分量的那个新点开始dfs，从而找到最大的矿石数。

吐槽：*可能传送到#，这里要注意，我就是被这个给坑啦。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<stack>
using namespace std;
int t,m,n;
char a[50][50];
int scnum,snum,ti,d[2000],dfn[2000],low[2000],k[2000],s[4000];
int num,fst[4000],next[8000000],node[8000000],ans;
bool inst[2000],vis[4000],p[4000][4000];
stack<int>st;
void init()
{
    int x,y;
    scnum=0;
    num=0;
    ti=0;
    ans=0;
    memset(fst,-1,sizeof(fst));
    memset(s,0,sizeof(s));
    memset(dfn,0,sizeof(dfn));
    memset(vis,0,sizeof(vis));
    memset(inst,0,sizeof(inst));
    memset(p,0,sizeof(p));
    memset(d,0,sizeof(d));
    scanf("%d%d",&m,&n);
    snum=m*n;
    for(int i=0; i<m; i++)scanf("%s",a[i]);
    for(int i=0; i<m; i++)
    {
        for(int j=0; j<n; j++)
        {
            if(a[i][j]!='#')
            {
                if(a[i][j]!='*')k[i*n+j]=(a[i][j]-'0');
                else
                {
                    k[i*n+j]=0;
                    scanf("%d%d",&x,&y);
                    if(a[x][y]!='#')
                    {
                        next[++num]=fst[i*n+j];
                        fst[i*n+j]=num;
                        node[num]=x*n+y;
                    }
                }
                if(i+1<m&&a[i+1][j]!='#')
                {
                    next[++num]=fst[i*n+j];
                    fst[i*n+j]=num;
                    node[num]=(i+1)*n+j;
                }
                if(j+1<n&&a[i][j+1]!='#')
                {
                    next[++num]=fst[i*n+j];
                    fst[i*n+j]=num;
                    node[num]=i*n+j+1;
                }
            }
        }
    }
}
void tarjan(int u)
{
    int v;
    dfn[u]=low[u]=++ti;
    st.push(u);
    inst[u]=1;
    for(int i=fst[u]; i!=-1; i=next[i])
    {
        v=node[i];
        if(!dfn[v])
        {
            tarjan(v);
            low[u]=min(low[u],low[v]);
        }
        else if(inst[v])low[u]=min(low[u],dfn[v]);
    }
    if(low[u]==dfn[u])
    {
        scnum++;
        do
        {
            v=st.top();
            st.pop();
            inst[v]=0;
            d[v]=snum+scnum;
            s[d[v]]+=k[v];
        }
        while(u!=v);
    }
}
void dfs(int u,int y)
{
    vis[u]=1;
    if(y+s[u]>ans)ans=y+s[u];
    for(int i=fst[u]; i!=-1; i=next[i])
    {
        int v=node[i];
        if(!vis[v])dfs(v,y+s[u]);
    }
    vis[u]=0;
}
void solve()
{
    int v;
    for(int i=0; i<snum; i++)
    {
        if(!dfn[i]&&a[i/n][i%n]!='#')tarjan(i);
    }
    for(int i=0; i<snum; i++)
    {
        if(d[i])
        {
            for(int j=fst[i]; j!=-1; j=next[j])
            {
                v=node[j];
                if(d[i]!=d[v]&&!p[d[i]][d[v]])
                {
                    p[d[i]][d[v]]=1;
                    next[++num]=fst[d[i]];
                    fst[d[i]]=num;
                    node[num]=d[v];
                }
            }
        }
    }
    dfs(d[0],0);
    printf("%d\n",ans);
}
int main()
{
    scanf("%d",&t);
    while(t--)
    {
        init();
        solve();
    }
    return 0;
}