首页 > 代码库 > HDU5093【最大流】

HDU5093【最大流】

题意:在*上建设炮塔,炮塔会像上下左右4个方向发射炮弹。o为浮冰,炮弹会越过。#为冰山,会阻挡炮弹。建设的炮塔会相互攻击,问最多建设多少个不互相攻击的炮塔。

本来我以为是贪心的,就像http://acm.hdu.edu.cn/showproblem.php?pid=1045一样,结果WA了,不懂是写错还是非法不行。

后来改成用http://acm.hdu.edu.cn/showproblem.php?pid=1045的网络流方法做就AC了。

建图:先用#把整个地图围一圈,把所有点分成ab两种点,源点向a建边,b点向汇点建边。对地图上的所有*,找它左端第一个#设为A,上端第一个#设为B,A的a点向B的b点建边。所建边容量都为1.

#include<iostream>#include<cstdio>#include<queue>using namespace std;const int NO=57;const int INF=1000000000;struct X{    int x,y;    X(){x=y=0;}    X(int a,int b){x=a;y=b;}}map[NO][NO],u[NO*NO*4],v[NO*NO*4],st,ed,r[NO][NO],c[NO][NO];char MAP[NO][NO];int n,m;int first[NO*2][NO*2],next[NO*NO*4],w[NO*NO*4],num;int dis[NO*2][NO*2];int ans[NO][NO];void reset_first(){    int i,j;    num=0;    for(i=0;i<=ed.x;i++)        for(j=0;j<=ed.y;j++)            first[i][j]=-1;}void add(int x1,int y1,int x2,int y2,int c){    u[num].x=x1,u[num].y=y1;    v[num].x=x2,v[num].y=y2;    w[num]=c;    next[num]=first[x1][y1];    first[x1][y1]=num++;    u[num].x=x2,u[num].y=y2;    v[num].x=x1,v[num].y=y1;    w[num]=0;    next[num]=first[x2][y2];    first[x2][y2]=num++;}void reset_dis(){    int i,j;    for(i=0;i<=ed.x;i++)        for(j=0;j<=ed.y;j++)            dis[i][j]=-1;}bool bfs(){    int i;    X a;    reset_dis();    queue<X> t;    t.push(st);    dis[st.x][st.y]=0;    while(!t.empty())    {        a=t.front();        t.pop();        for(i=first[a.x][a.y];i!=-1;i=next[i])            if(w[i]&&dis[v[i].x][v[i].y]==-1)            {                dis[v[i].x][v[i].y]=dis[a.x][a.y]+1;                t.push(v[i]);            }    }    return dis[ed.x][ed.y]!=-1;}int min(int a,int b){return a<b?a:b;}int dfs(const X &k,int MIN){    if(k.x==ed.x&&k.y==ed.y)        return MIN;    int sum=0,a,i;    for(i=first[k.x][k.y];i!=-1&&sum<MIN;i=next[i])        if(w[i]&&dis[v[i].x][v[i].y]==dis[k.x][k.y]+1&&(a=dfs(v[i],min(MIN-sum,w[i]))))        {            w[i]-=a;            w[i^1]+=a;            sum+=a;        }    if(!sum)        dis[k.x][k.y]=-1;    return sum;}int DANIC(){    int ans=0,a;    while(bfs())        while(a=dfs(st,INF))            ans+=a;    return ans;}void read(){    for(int i=0;i<=n+1;i++)        MAP[i][0]=MAP[i][m+1]=#;    for(int j=0;j<=m+1;j++)        MAP[0][j]=MAP[n+1][j]=#;    for(int i=1;i<=n;i++)    {        getchar();        for(int j=1;j<=m;j++)            MAP[i][j]=getchar();    }    for(int i=0;i<=n+1;i++)        for(int j=0;j<=m+1;j++)            if(MAP[i][j]==#)            {                ans[i][j]=-1;                map[i][j].x=MAP[i+1][j]!=#;                map[i][j].y=MAP[i][j+1]!=#;            }            else                ans[i][j]=map[i][j].x=map[i][j].y=0;}void build(){    int i,j,k;    for(i=0;i<=n;i++)        for(j=0;j<=m;j++)        {            if(map[i][j].x)            {                for(k=0;i+k+1<=n;k++)                    if(ans[i+k+1][j]==-1)                        break;                    else                        c[i+k+1][j].x=i,c[i+k+1][j].y=j;                add(n+i,m+j,ed.x,ed.y,1);            }            if(map[i][j].y)            {                for(k=0;j+k+1<=m;k++)                    if(ans[i][j+k+1]==-1)                        break;                    else                        r[i][j+k+1].x=i,r[i][j+k+1].y=j;                add(st.x,st.y,i,j,1);            }        }    for(i=1;i<=n;i++)        for(j=1;j<=m;j++)            if(MAP[i][j]==*)                add(r[i][j].x,r[i][j].y,n+c[i][j].x,m+c[i][j].y,1);}int main(){    int ttt;    scanf("%d",&ttt);    while(ttt--)    {        scanf("%d%d",&n,&m);        ed.x=n*2+1;ed.y=m*2+1;        reset_first();        read();        build();        printf("%d\n",DANIC());    }    return 0;}
View Code

 

HDU5093【最大流】