BZOJ 2150 cogs 1861 [国家集训队2011]部落战争

简介: 题目描述 lanzerb的部落在A国的上部,他们不满天寒地冻的环境,于是准备向A国的下部征战来获得更大的领土。 A国是一个M*N的矩阵,其中某些地方是城镇,某些地方是高山深涧无人居住。lanzerb把自己的部落分成若干支军队,他们约定: 1. 每支军队可以从任意一个城镇出发,并只能从上往向下征战,不能回头。

题目描述

lanzerb的部落在A国的上部,他们不满天寒地冻的环境,于是准备向A国的下部征战来获得更大的领土。 A国是一个M*N的矩阵,其中某些地方是城镇,某些地方是高山深涧无人居住。lanzerb把自己的部落分成若干支军队,他们约定: 1. 每支军队可以从任意一个城镇出发,并只能从上往向下征战,不能回头。途中只能经过城镇,不能经过高山深涧。 2. 如果某个城镇被某支军队到过,则其他军队不能再去那个城镇了。 3. 每支军队都可以在任意一个城镇停止征战。 4. 所有军队都很奇怪,他们走的方法有点像国际象棋中的马。不过马每次只能走1*2的路线,而他们只能走R*C的路线。 lanzerb的野心使得他的目标是统一全国,但是兵力的限制使得他们在配备人手时力不从心。假设他们每支军队都能顺利占领这支军队经过的所有城镇,请你帮lanzerb算算至少要多少支军队才能完成统一全国的大业。

输入

第一行包含4个整数M、N、R、C,意义见问题描述。接下来M行每行一个长度为N的字符串。如果某个字符是'.',表示这个地方是城镇;如果这个字符时'x',表示这个地方是高山深涧。

输出

输出一个整数,表示最少的军队个数。

样例输入

【样例输入一】
3 3 1 2
...
.x.
...
【样例输入二】
5 4 1 1
....
..x.
...x
....
x...

样例输出

【样例输出一】
4

【样例输出二】
5
【样例说明】

【数据范围】
100%的数据中,1<=M,N<=50,1<=R,C<=10。

提示

来源

  //BZOJ没有,我来加个:国家集训队2011

吐槽

  n和m都小于50,空间至少要开到50*50=2500!

解题思路

  一道最小路径覆盖的题目,可以按照题意建图,对于点$(i,j)$,向$(i+r,j+c)$、$(i+r,j-c)$、$(i+c,j+r)$、$(i+c,j-r)$连单向边(前提是点在地图上且是城市),然后跑最小路径覆盖即可,我依然跑匈牙利,这提数据范围用不到dinic,但是dinic跑得快。

源代码

#include<cstdio>
#include<vector>
#include<cstring>
int m,n,r,c,ans=0,N=0;
char ch;
bool g[2510][2510]={0};//(x-1)*n+y
bool vis[2510]={0};
struct men{
    int lover;
    std::vector<int> t;
}man[2510];

int woman[2510]={0};

bool dfs(int u)
{
    int sz=man[u].t.size();
    for(int i=0;i<sz;i++)
    {
        int w=man[u].t[i];
        if(!vis[w])
        {
            vis[w]=1;
            if(!woman[w]||dfs(woman[w]))
            {
                woman[w]=u;
                man[u].lover=w;
                return 1;
            }
        }
    }
    return 0;
}
int main()
{
    //freopen("lanzerb.in","r",stdin);
    //freopen("lanzerb.out","w",stdout);
    scanf("%d %d %d %d",&m,&n,&r,&c);
    for(int i=1;i<=m;i++)
    {
        for(int j=1;j<=n;j++)
        {
            scanf("\n%c",&ch);
            if(ch=='.')
                g[i][j]=1,N++;
        }
    }
    int bh[4][2]={{r,c},{c,-r},{c,r},{r,-c}};
    for(int i=1;i<=m;i++)
    {
        for(int j=1;j<=n;j++)
        {
            if(!g[i][j]) continue;
            int id=(i-1)*n+j;
            for(int k=0;k<4;k++)
            {
                if(r==c&&k>1) break;
                int x=i+bh[k][0],y=j+bh[k][1];
                if(x>=1&&x<=m&&y>=1&&y<=n&&g[x][y])
                {
                    man[id].t.push_back((x-1)*n+y);
                }
            }
        }
    }
    memset(woman,0,sizeof(woman));
    for(int i=1;i<=n*m;i++)
    {
        memset(vis,0,sizeof(vis));
        if(dfs(i))
            ans++;
    }
    printf("%d\n",N-ans);
    return 0;
}

 

目录
相关文章
|
3月前
|
Shell
[SWPUCTF 2021 新生赛]gift_pwn-入土为安的第十五天
[SWPUCTF 2021 新生赛]gift_pwn-入土为安的第十五天
128 0
|
3月前
日更[ACTF新生赛2020]Oruga wp
日更[ACTF新生赛2020]Oruga wp
31 0
|
6月前
|
算法
【随想】每日两题Day.9
【随想】每日两题
36 1
|
6月前
【随想】每日两题Day.8
【随想】每日两题
35 0
|
6月前
【随想】每日两题Day.17
【随想】每日两题Day.17
39 0
|
6月前
|
算法
【随想】每日两题Day.1
【随想】每日两题Day.1
33 0
|
6月前
【随想】每日两题Day.22
【随想】每日两题Day.22
33 0
|
6月前
【随想】每日两题Day.19
【随想】每日两题Day.19
36 0
|
6月前
|
数据安全/隐私保护
【随想】每日两题Day.14
【随想】每日两题Day.14
33 0
|
6月前
【随想】每日两题Day.4
【随想】每日两题Day.4
28 0