poj 3592 强连通+缩点+dfs

来源:互联网 发布:域控推送软件 编辑:程序博客网 时间:2024/05/20 06:23

传送们

题意:给你一个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;}


原创粉丝点击