hdoj4859海岸线

1.地图周围再加一圈海

2.周长最多为sum=n*(m+1)+m*(n+1)。如果有邻接相同，要减1.最小割使相同最少。结果为sum-最小割

3.但是有E海滩，两边都能选。让E到S,T都连一个很大的数（10000，比最小割大就行）。这样最小割时两个10000一定要选一个。那么在sum中也加上10000就能消除。这样一定割一边而且不都割。而且连S,和T都为一样的数能确保消去E的决策的干扰。

剩下的就和套路一样了，奇偶建图，偶数时,'.'连S     INF,'D'连T   INF,E连S,T     10000。邻接互相连1。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <vector>
#include <iomanip>
#include <cstring>
#include <map>
#include <queue>
#include <set>
#include <cassert>
#include <stack>
#include <bitset>
#define mkp make_pair
using namespace std;
const double EPS=1e-;
typedef long long lon;
const lon SZ=,SSZ=*SZ*SZ,APB=,one=,INF=0x7FFFFFFF,mod=;
lon n,m,S=,T=,sum,cnt,head[SZ*SZ];
lon nex[SSZ],to[SSZ],wt[SSZ],dep[SSZ];
int dx[]={,},dy[]={,};
char ch[SZ][SZ];
 
void add(int u,int v,int w)
{
    ++cnt;
    nex[cnt]=head[u];
    head[u]=cnt;
    to[cnt]=v,wt[cnt]=w;
}
 
int getid(int i,int j)
{
    return (i-)*(m+)+j;
}
 
void init()
{
    cin>>n>>m;
    sum=n*(m+)+m*(n+),cnt=-;
    memset(head,-,sizeof(head));
    for(int i=;i<=n+;++i)
    {
        for(int j=;j<=m+;++j)
        {
            if(i==||j==||i==n+||j==m+)
            {
                ch[i][j]='D';
            }
            else
            {
                cin>>ch[i][j];
                if(ch[i][j]=='E')sum+=APB;
            }
        }
    }
    for(int i=;i<=n+;++i)
    {
        for(int j=;j<=m+;++j)
        {
            if(!((i+j)&))
            {
                if(ch[i][j]=='.')
                {
                    add(S,getid(i,j),INF);
                    add(getid(i,j),S,);
                }
                else if(ch[i][j]=='D')
                {
                    add(getid(i,j),T,INF);
                    add(T,getid(i,j),);
                }
                else
                {
                    add(S,getid(i,j),APB);
                    add(getid(i,j),S,);
                    add(getid(i,j),T,APB);
                    add(T,getid(i,j),);
                }
            }
            else
            {
                if(ch[i][j]=='.')
                {
                    add(getid(i,j),T,INF);
                    add(T,getid(i,j),);
                }
                else if(ch[i][j]=='D')
                {
                    add(S,getid(i,j),INF);
                    add(getid(i,j),S,);
                }
                else
                {
                    add(S,getid(i,j),APB);
                    add(getid(i,j),S,);
                    add(getid(i,j),T,APB);
                    add(T,getid(i,j),);
                }
            }
        }
    }
    for(int i=;i<=n+;++i)
    {
        for(int j=;j<=m+;++j)
        {
            for(int k=;k<;++k)
            {
                int nx=i+dx[k],ny=j+dy[k];
                if(i==&&nx==||j==&&ny==)continue;
                add(getid(i,j),getid(nx,ny),);
                add(getid(nx,ny),getid(i,j),);
                add(getid(nx,ny),getid(i,j),);
                add(getid(i,j),getid(nx,ny),);
            }
        }
    }
//    for(int i=1;i<=n+2;++i)
//    {
//        for(int j=1;j<=m+2;++j)
//        {
//            for(int k=head[getid(i,j)];k!=-1;k=nex[k])
//            {
//                cout<<i<<" "<<j<<" "<<to[k]<<" "<<wt[k]<<endl;
//            }
//        }
//    }
}
 
bool bfs()
{
    memset(dep,,sizeof(dep));
    dep[S]=;
    queue<int> q;
    q.push(S);
    for(;q.size();)
    {
        int fr=q.front();
        q.pop();
        for(int i=head[fr];i!=-;i=nex[i])
        {
            int t=to[i],w=wt[i];
            if(!dep[t]&&w)
            {
                dep[t]=dep[fr]+;
                q.push(t);
                if(t==T)return ;
            }
        }
    }
    return ;
}
 
int dinic(int x,int flow)
{
    if(x==T)return flow;
    else
    {
        int rem=flow;
        for(int i=head[x];i!=-&&rem;i=nex[i])
        {
            int t=to[i],w=wt[i];
            if(dep[t]==dep[x]+&&w)
            {
                int tmp=dinic(t,min(w,rem));
                if(!tmp)dep[t]=;
                rem-=tmp;
                wt[i]-=tmp;
                wt[i^]+=tmp;
            }
        }
        return flow-rem;
    }
}
 
void work()
{
    int res=;
    for(;bfs();)
    {
        res+=dinic(S,INF);
    }
    cout<<sum-res<<endl;
}
 
void release()
{
 
}
 
int main()
{
    std::ios::sync_with_stdio();
    //freopen("d:\\1.txt","r",stdin);
    lon casenum;
    cin>>casenum;
    //cout<<casenum<<endl;
    for(int time=;time<=casenum;++time)
    //for(int time=1;cin>>n>>m,n;++time)
    {
        cout<<"Case "<<time<<": ";
        init();
        work();
        release();
    }
    return ;
}