bzoj 1066: [SCOI2007]蜥蜴

news/2024/8/26 14:40:28

Description

在一个r行c列的网格地图中有一些高度不同的石柱,一些石柱上站着一些蜥蜴,你的任务是让尽量多的蜥蜴逃到边界外。 每行每列中相邻石柱的距离为1,蜥蜴的跳跃距离是d,即蜥蜴可以跳到平面距离不超过d的任何一个石柱上。石柱都不稳定,每次当蜥蜴跳跃时,所离开的石柱高度减1(如果仍然落在地图内部,则到达的石柱高度不变),如果该石柱原来高度为1,则蜥蜴离开后消失。以后其他蜥蜴不能落脚。任何时刻不能有两只蜥蜴在同一个石柱上。

Input

输入第一行为三个整数r,c,d,即地图的规模与最大跳跃距离。以下r行为石竹的初始状态,0表示没有石柱,1~3表示石柱的初始高度。以下r行为蜥蜴位置,“L”表示蜥蜴,“.”表示没有蜥蜴。

Output

输出仅一行,包含一个整数,即无法逃离的蜥蜴总数的最小值。

Sample Input

5 8 2
00000000
02000000
00321100
02000000
00000000
........
........
..LLLL..
........
........

Sample Output

1

HINT

100%的数据满足:1<=r, c<=20, 1<=d<=3


很裸的拆点网络流。。每个点拆成u和v

v连可以到达的u,容量INF。如果当前点可以跳出边界则v连汇点,容量INF

有蜥蜴的地方从源点连u,容量为1。


T T拆点的序号弄错了结果找了半个多小时错。。


#include<cstdio>
#include<string>
#include<cstring>
using namespace std;
int head[100001];
struct map
{
	 int f;
	 int s,t;
     int next;
}a[400001];
int b[5001];
int edge;
int p;
int q[400001],d[400001];
inline void add(int s,int t,int f)
{
     a[edge].next=head[s];
     head[s]=edge;
     a[edge].s=s;
     a[edge].t=t;
     a[edge].f=f;
}
inline bool bfs()
{
     int l=0,r=0;
     memset(q,0,sizeof(q));
     r++;
     q[r]=0;
     memset(d,-1,sizeof(d));
     d[0]=0;
     int i,k;
     while(l<r)
     {
     	  l++;
          int k=q[l];
          for(i=head[k];i!=0;i=a[i].next)
          {
               if(a[i].f>0&&d[a[i].t]==-1)
               {
                    d[a[i].t]=d[k]+1;
                    r++;
                    q[r]=a[i].t;
               }
          }
     }
     if(d[p]>=0)
          return true;
     return false;
}
inline int dfs(int k,int s)
{
     if(k==p)
          return s;
     int t=s;
     int i;
     for(i=head[k];i!=0;i=a[i].next)
     {
          if(d[a[i].t]==d[k]+1&&a[i].f>0)
          {
               int xx=dfs(a[i].t,min(s,a[i].f));
               a[i].f-=xx;
               if(i%2==0)
                    a[i-1].f+=xx;
               else
                    a[i+1].f+=xx;
               s-=xx;
          }
     }
     return t-s;
}
inline int maxflow()
{
     int s=0;
     while(bfs())
          s+=dfs(0,2100000000);
     return s;
}
inline int abs(int x)
{
     if(x<0)
          x=-x;
     return x;
}
int map[21][21];
int main()
{
     int r,c,d;
     scanf("%d%d%d",&r,&c,&d);
     p=r*c*2+1;
     int pp=r*c;
     int i,j;
     char x;
     for(i=1;i<=r;i++)
     {
          scanf("%c",&x);
          for(j=1;j<=c;j++)
          {
               scanf("%c",&x);
               map[i][j]=x-'0';
               if(map[i][j]!=0)
               {
                    edge++;
                    add((i-1)*c+j,(i-1)*c+j+pp,map[i][j]);
                    edge++;
                    add((i-1)*c+j+pp,(i-1)*c+j,0);
               }
          }
     }
     int ii,jj;
     for(i=1;i<=r;i++)
     {
          for(j=1;j<=c;j++)
          {
          	   if(map[i][j]>0)
          	   {
                    for(ii=1;ii<=r;ii++)
                    {
                         for(jj=1;jj<=c;jj++)
                         {
                       	      if(map[ii][jj]>0)
                       	      {
                                   if((i-ii)*(i-ii)+(j-jj)*(j-jj)<=d*d)
                                   {
                                        edge++;
                                        add((i-1)*c+j+pp,(ii-1)*c+jj,2100000000);
                                        edge++;
                                        add((ii-1)*c+jj,(i-1)*c+j+pp,0);
                                   }
                              }
                         }
                    }
                    if(i<=d||r-i<d||j<=d||c-j<d)
                    {
                         edge++;
                         add((i-1)*c+j+pp,p,2100000000);
                         edge++;
                         add(p,2*((i-1)*c+j),0);
                    }
               }
          }
     }
     int sum=0;
     for(i=1;i<=r;i++)
     {
          scanf("%c",&x);
          for(j=1;j<=c;j++)
          {
               scanf("%c",&x);
               if(x=='L')
               {
                    sum++;
                    edge++;
                    add(0,(i-1)*c+j,1);
                    edge++;
                    add((i-1)*c+j,0,0);
               }
          }
     }
     printf("%d\n",sum-maxflow());
     return 0;
}



http://www.niftyadmin.cn/n/943253.html

相关文章

组合数的一些模板(好用)

一、 ll qpow(ll a,ll b,ll p) {ll ret1;a%p;while(b){if(b&1) retret*a%p;b/2;aa*a%p;}return ret; } 二、Lucas ll lucas(ll n,ll m,ll p) {if(m0) return 1;return C(n%p,m%p,p)*lucas(n/p,m/p,p)%p; } 三、C(n,m) 方案一&#xff1a;&#xff08;直接算&#xff0…

bzoj 1934: [Shoi2007]Vote 善意的投票

Description 幼儿园里有n个小朋友打算通过投票来决定睡不睡午觉。对他们来说&#xff0c;这个问题并不是很重要&#xff0c;于是他们决定发扬谦让精神。虽然每个人都有自己的主见&#xff0c;但是为了照顾一下自己朋友的想法&#xff0c;他们也可以投和自己本来意愿相反的票。我…

bzoj 1706: [usaco2007 Nov]relays 奶牛接力跑

Description FJ的N(2 < N < 1,000,000)头奶牛选择了接力跑作为她们的日常锻炼项目。至于进行接力跑的地点 自然是在牧场中现有的T(2 < T < 100)条跑道上。 农场上的跑道有一些交汇点&#xff0c;每条跑道都连结了两个不同的交汇点 I1_i和I2_i(1 < I1_i < 1,0…

杭电ACM——1247,Hat’s Words(Trie树)

这道题是有问题的&#xff0c;一个单词是可以被另外一个单词拼接两次得到的&#xff0c;如ab&#xff0c;abab&#xff0c;输出是abab&#xff0c;然而原题确指明&#xff0c;某个单词当且仅当是另外字典里的两个单词拼接的&#xff0c;才是一个hat’s word&#xff0c;实际上这…

bzoj 1001: [BeiJing2006]狼抓兔子

Description 现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到&#xff0c;但抓兔子还是比较在行的&#xff0c;而且现在的兔子还比较笨&#xff0c;它们只有两个窝&#xff0c;现在你做为狼王&#xff0c;面对下面这样一个网格的地形&#xff1a; 左上角点…

杭电ACM——2533,N皇后问题(搜索)

突破口&#xff1a; N*N的棋盘放置N个皇后&#xff0c;因为皇后的特性&#xff0c;合法的放置方法中每行&#xff08;每列&#xff09;必定有一个皇后&#xff0c;因此我们可以先搜索第1行可以放置的位置&#xff0c;再逐步搜索第2行&#xff0c;直到第N行。而不用去穷举第一个…

杭电ACM——2546,饭卡(DP)

一维背包问题。 突破口&#xff1a;将菜的价格从小到大排序&#xff0c;对前n-1个物品进行DP&#xff0c;找出消费金额不超过m-5的方案中最大的消费金额sum&#xff0c;在此基础上&#xff0c;m - (sumprice[n]) 即为卡内最小余额。 因此&#xff0c;我们DP应在余额为0~m-5之内…

bzoj 1797: [Ahoi2009]Mincut 最小割

Description A,B两个国家正在交战&#xff0c;其中A国的物资运输网中有N个中转站&#xff0c;M条单向道路。设其中第i (1≤i≤M)条道路连接了vi,ui两个中转站&#xff0c;那么中转站vi可以通过该道路到达ui中转站&#xff0c;如果切断这条道路&#xff0c;需要代价ci。现在B国想…