【模拟】bzoj1686: [Usaco2005 Open]Waves 波纹

打完模拟题来庆祝一波：)；感觉最近陷入一种“口胡五分钟打题两小时”的巨坑之中……

Description

Input

    第1行：四个用空格隔开的整数Pj Bi，B2，R. P(1≤P≤5)表示石子的个数，Bi(-5×100000≤Bi≤5×100000)和B2（-5×100000≤B2≤5×100000）表示两个河堤的z坐标，R(1≤R≤5×100000)表示你要描述湖面多少秒．没有两个石子会在同一时间砸到同一地点，两个河堤一定有不同的坐标，没有石子会砸到河堤上去．

    第2到P+1行：每行有三个用空格隔开的整数描述了一颗石子，X，K T(-5×100000≤X，K T≤5×100000).X．y表示石子砸的地点的坐标，T表示石子是什么时候砸下去的．

Output

    输出是一个9半9的矩阵，中心在(0，0)点．左下点的坐标为(-4，-4)，右上点的坐标为(4，4)．这个矩阵表现的是R秒时湖面状态．

---

题目分析

有用的信息只有9x9的这一部分，所以我们来对于每一颗石子分别处理。

如果只是手算前两三个石头，看上去好像状态是指数阶的感觉很吓人，但实际上显然状态是与$n$同阶的（为什么我第一眼没看出来呢……）。

可以发现每一个水波到中心的曼哈顿距离相同，那么枚举水波就变得很容易了，我们来考虑如何计算反射。

注意到反射时$y$坐标是不会改变的，那只需要考虑起点和终点的$x$坐标。

有如下三种情况：

1.大坝包含了两点，所以传递没有受到任何影响。

2.大坝包住了起点，那么比较麻烦，需要考虑多次反射。

3.大坝包住了终点，那么只用考虑一次反射。

应该是有更优的写法的，不过分类讨论更稳（？）一些吧。

之后嘛，注意不要打挂就好了。

 #include<bits/stdc++.h>
 
 struct node
 {
     int a[][];
     int *const operator[](int x)
     {
         return a[x+];
     }
     node() {}
 }f;
 int stones,lRes,rRes,times;
 int s;
 
 int read()
 {
     char ch = getchar();
     int num = ;
     bool fl = ;
     for (; !isdigit(ch); ch = getchar())
         if (ch=='-') fl = ;
     for (; isdigit(ch); ch = getchar())
         num = (num<<)+(num<<)+ch-;
     if (fl) num = -num;
     return num;
 }
 inline bool limit(int x, int y){return x>=-&&x<=&&y>=-&&y<=;}
 inline bool illegal(int x)
 {
     if (x<-) return (x<lRes&&lRes<=-)||(x<rRes&&rRes<=-);
     if (x>) return (x>lRes&&lRes>=)||(x>rRes&&rRes>=);
     return ;
 }
 void add(int x, int y, int c)
 {
     register int bar,cnt;
     if (y < -||y > ) return;
     if (s < x){
         if ((s<=lRes&&lRes<=x)||(s<=rRes&&rRes<=x)){
             if (s <= lRes){
                 bar = lRes, x = *bar-x-;
             }else{
                 cnt = ;
                 while (x <= lRes||x >= rRes)
                 {
                     if (cnt) bar = rRes;
                     else bar = lRes;
                     x = *bar-x, cnt = -cnt;
                     if (cnt) x++;
                     else x--;
                 }
             }
         }
     }
     if (s > x){
         if ((x<=lRes&&lRes<=s)||(x<=rRes&&rRes<=s)){
             if (s >= rRes){
                 bar = rRes, x = *bar-x+;
             }else{
                 cnt = ;
                 while (x <= lRes||x >= rRes)
                 {
                     if (cnt) bar = lRes;
                     else bar = rRes;
                     x = *bar-x, cnt = -cnt;
                     if (cnt) x--;
                     else x++;
                 }
             }
         }
     }
     if (limit(x, y)) f[x][y] += c;
 }
 void deal(int x, int y, int t)
 {
     if (t==){
         if (limit(x, y)) f[x][y]++;
         return;
     }
     s = x;
     for (int delta=; delta<=t; delta++)
     {
         add(x+delta, y+t-delta, ), add(x+delta, y-t+delta, );
         add(x-delta, y+t-delta, ), add(x-delta, y-t+delta, );
     }
     add(x, y+t, -), add(x, y-t, -);
     add(x+t, y, -), add(x-t, y, -);
     t -= ;
     if (t==){
         if (limit(x, y)) f[x][y]--;
         return;
     }else if (t < ) return;
     for (int delta=; delta<=t; delta++)
     {
         add(x+delta, y+t-delta, -), add(x+delta, y-t+delta, -);
         add(x-delta, y+t-delta, -), add(x-delta, y-t+delta, -);
     }
     add(x, y+t, ), add(x, y-t, );
     add(x+t, y, ), add(x-t, y, );
 }
 int main()
 {
     stones = read(), lRes = read(), rRes = read(), times = read();
     if (lRes > rRes) std::swap(lRes, rRes);
     for (int i=; i<=stones; i++)
     {
         int x = read(), y = read(), t = read();
         if (t > times||illegal(x)) continue;
         deal(x, y, times-t);
     }
     for (int j=; j>=-; j--)
     {
         for (int i=-; i<=; i++)
         {
             if (i==lRes||i==rRes) putchar('X');
             else if (f[i][j] < ) putchar('o');
             else if (f[i][j] > ) putchar('*');
             else putchar('-');
         }
         putchar('\n');
     }
     return ;
 }

后记

逛了一圈发现我的代码又长又慢？？？

END