BZOJ2548：[CTSC2002]灭鼠行动
我对模拟的理解：https://www.cnblogs.com/AKMer/p/9064018.html
题目传送门：https://www.lydsy.com/JudgeOnline/problem.php?id=2548
原本计划先在三国杀赢了那几头猪再来灭这群老鼠的……
然后因为被猪吊锤了一周……然后特工队也等不下去了，我就被调来先灭老鼠了……
灭完老鼠继续跟猪打三国杀去\(emmm\)……
\([SDOI2010]\)猪国杀和\([CTSC2002]\)灭鼠行动就是我模拟旅程的最后两站了，弄完就不专门搞模拟了。
由于有[ZJOI2008]杀蚂蚁这题作为前提，我对于这种在地图上杀什么玩意儿的题理解深了很多，所以写灭鼠行动的时候比杀蚂蚁轻松多了。很多人说灭鼠行动比杀蚂蚁难，我觉得并不然。
首先，让我们把研究的重点放在炸弹上：
炸弹放下即爆，或者直接晕和变性（听起来好变态）
定时炸弹相当于一个即爆炸弹，我们可以假装它其实是t+3才放下来的。
然后就没什么别的了。
然后我们再把目光放到老鼠身上来：
老鼠移动算法比蚂蚁移动算法简单多了，能往前就往前，不然就拐弯，拐弯判一判就没了。
然后繁殖的时候要注意，只有在某个格子上只有一雄一雌两只成年老鼠才会生仔（他们需要私人空间，有其它老鼠在是不会繁殖的）
然后繁殖以及休息相当于眩晕三秒，第二秒后小老鼠就可以出来了。
如果在生产过程中被眩晕，那么小老鼠也要推迟时间出来。
怎么都觉得比杀蚂蚁简单多了……唉，不说了，我跟猪打三国杀去了。
时间复杂度：\(O(Tnm)\)
空间复杂度：\(O(不会炸)\)
代码如下：
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define sqr(x) ((x)*(x))
const int inf=2e9;
char s[10];//读入字符串用
bool explode[55][55];//explode[x][y]表示(x,y)是否在爆炸范围内
int dx[4]={-1,0,1,0};
int dy[4]={0,1,0,-1};//四个方向，分别是北，东，南，西，对应编号0,1,2,3，刚好对应读入代表方向的数字2^0,2^1,2^2,2^3
int left[4]={3,0,1,2};//left[x]表示x方向的左边方向编号
int right[4]={1,2,3,0};//right[x]表示x方向的右边编号
int fake[4]={1,2,4,8};//fake[x]表示2^x，也就是读入的网格信息里表示方向的数字
int now_time,now_bomb=1;//now_time表示现在时间，now_bomb表示现在轮到哪一个炸弹了
int con[55][55],mice_sum[55][55];//con[x][y]记录(x,y)与四周的连通性，mice_sum[x][y]记录(x,y)上一共多少老鼠
int len,radius,n,m,tot,limit,Time,cnt;//len记录1号炸弹范围，radius记录2号炸弹范围半径，n，m记录地图大小，tot实时更新记录老鼠个数，limit为评判鼠疫是否爆发的标准，Time是要求模拟的时间，cnt是炸弹个数
int read() {
int x=0,f=1;char ch=getchar();
for(;ch<'0'||ch>'9';ch=getchar())if(ch=='-')f=-1;
for(;ch>='0'&&ch<='9';ch=getchar())x=(x<<1)+(x<<3)+ch-'0';
return x*f;
}//快读
struct mouse {
bool dead,turn_left;
int x,y,dir,sex,ver,bre,grow;
//dead记录老鼠是否死亡，turn_left记录老鼠碰到前方无路左右都有路是否要往左转，x，y记录老鼠坐标，dir记录方向，sex记录性别，ver记录眩晕时间，bre记录小老鼠出生的时间，grow记录还有多久成年
void build() {
x=read(),y=read();
scanf("%s",s);
if(s[0]=='N')dir=0;
if(s[0]=='E')dir=1;
if(s[0]=='S')dir=2;
if(s[0]=='W')dir=3;
turn_left=1;//奇数次往左拐，所以初始为true
scanf("%s",s);
if(s[0]=='X')sex=1;
bre=-1;
}//读入老鼠信息
void move() {
bre=-1;//能移动了说明没有预产期，把预产期恢复成-1
if(con[x][y]&fake[dir]) {
x+=dx[dir],y+=dy[dir];
return;//如果能在原方向上向前走就向前走
}
if(fake[left[dir]]&con[x][y]) {//如果能向左走
if(fake[right[dir]]&con[x][y]) {//如果能向右走
if(turn_left)dir=left[dir];
else dir=right[dir];//奇数次向左，偶数次向右
turn_left^=1;//改变下一次拐弯的方向
}
else dir=left[dir];//只能向左走，花一单位时间转向
return;
}
if(fake[right[dir]]&con[x][y]) {
dir=right[dir];
return;//只能向右走，花一单位时间转向
}
dir=left[dir];//掉头相当于往左转两次
}
}M[500];
struct bomb {
int type,tim,x,y;
//type记录炸弹类型，tim记录炸弹爆炸时间，x，y记录炸弹坐标
bool operator<(const bomb &a)const {
return tim<a.tim;
}//按时间排序
void build() {
type=read(),tim=read();
x=read(),y=read();
if(type==3)tim+=3;//如果是三号炸弹那么时间推迟3个单位
}
void mark_1(int x,int y) {//标记1号炸弹爆炸范围
explode[x][y]=1;//炸弹所在点标记
for(int i=0;i<4;i++) {//往四个方向走
int now_x=x,now_y=y,step_sum=0;
//step_sum表示从(x,y)到(now_x,now_y)距离多少
while(step_sum<len&&(fake[i]&con[now_x][now_y]))
now_x+=dx[i],now_y+=dy[i],explode[now_x][now_y]=1,step_sum++;//如果能继续深入就继续深入
}
}
void mark_2(int x,int y) {//2号炸弹标记范围
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
if(sqr(i-x)+sqr(j-y)<=radius)
explode[i][j]=1;//勾股定理
}
void blast() {
memset(explode,0,sizeof(explode));//清空标记
if(type==1)mark_1(x,y);
if(type==2)mark_2(x,y);
if(type==3)explode[x][y]=1;
if(type==4)explode[x][y]=1;//标记范围
for(int i=1;i<=tot;i++)
if(explode[M[i].x][M[i].y]) {//如果在范围内
if(type==1||type==3)M[i].dead=1;//如果是1号或者3号炸弹直接秒了老鼠
if(type==2) {
M[i].ver+=3;
if(M[i].bre>now_time)M[i].bre+=3;//2号炸弹就加眩晕，如果有预产期预产期也跟着推后
}
if(type==4)M[i].sex^=1;//4号炸弹变性
}
}
}B[105];
void check_dead() {
int top=0;
for(int i=1;i<=tot;i++)
if(!M[i].dead)
M[++top]=M[i];
tot=top;//清理死亡老鼠尸体
}
void breed() {
memset(mice_sum,0,sizeof(mice_sum));//清空mice_sum[][]
for(int i=1;i<=tot;i++)
mice_sum[M[i].x][M[i].y]++;//更新mice_sum[][]
for(int i=1;i<=tot;i++)
if(!M[i].ver&&M[i].bre==-1&&!M[i].grow&&M[i].sex)//如何i号老鼠没有眩晕时间并且可以繁衍，也是成年的雄老鼠
for(int j=1;j<=tot;j++)
if(!M[j].ver&&M[j].bre==-1&&!M[j].grow&&!M[j].sex)
//j号老鼠也是没有眩晕的成年可生育老鼠，并且是雌的 if(mice_sum[M[i].x][M[i].y]==2&&M[i].x==M[j].x&&M[i].y==M[j].y){//如果他们在同一个格子上并且没有其它老鼠在场就开始造老鼠
M[i].ver+=3;M[j].ver+=3;//眩晕时间+3
M[i].bre=M[j].bre=now_time+2;//2单位时间后生出小老鼠
break;
}
for(int i=1;i<=tot;i++)
if(M[i].bre==now_time&&M[i].sex)//如果现在的时间这只老鼠会生育并且他是雄的，每一对老鼠生育只算一次
for(int j=0;j<4;j++)
if(fake[j]&con[M[i].x][M[i].y]) {//如果该点有该方向的管道那么生成一只小老鼠
tot++;
M[tot].x=M[i].x;
M[tot].y=M[i].y;
M[tot].ver=0;
M[tot].dir=j;
M[tot].grow=5;
M[tot].bre=-1;
M[tot].dead=0;
M[tot].turn_left=1;
if(j==0||j==2)M[tot].sex=1;
else M[tot].sex=0;
//记得所有数据全部都要初始化，因为tot+1保不定存了某只死亡的老鼠的信息，不全部初始化就GG了
}
}
int main() {
len=read(),radius=read(),n=read(),m=read();
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
con[i][j]=read();
tot=read();radius*=radius;//读入，半径先平方一遍，到时候直接上勾股定理
for(int i=1;i<=tot;i++)
M[i].build();//读入老鼠信息
cnt=read(),limit=read();
for(int i=1;i<=cnt;i++)
B[i].build();//读入炸弹
sort(B+1,B+cnt+1);B[cnt+1].tim=inf;//炸弹排序
Time=read();//读入时间
for(now_time=0;now_time<=Time;now_time++) {
while(B[now_bomb].tim==now_time)
B[now_bomb].blast(),now_bomb++;//把现在这个时间会炸的炸弹点爆。每一秒最先发生的事情是炸炸弹而不是老鼠繁殖和移动，这一点我改了好久……把它当杀蚂蚁了……
check_dead();//清理尸体
breed();//繁殖
if(tot>limit) {puts("-1");return 0;}//如果鼠疫爆发就直接没得玩了
if(now_time==Time)break;//时间过完了就不移动了。
for(int i=1;i<=tot;i++) {
if(M[i].ver)
M[i].ver--;//如果眩晕就减一秒眩晕时间
else {
M[i].move();//移动
if(M[i].grow)M[i].grow--;//小老鼠长大一秒
}
}
}
printf("%d\n",tot);//输出老鼠总数
return 0;
}
