【模拟】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.大坝包住了终点,那么只用考虑一次反射。

应该是有更优的写法的,不过分类讨论更稳(?)一些吧。

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

 

  1 #include<bits/stdc++.h>
  2 
  3 struct node
  4 {
  5     int a[103][103];
  6     int *const operator[](int x)
  7     {
  8         return a[x+50];
  9     }
 10     node() {}
 11 }f;
 12 int stones,lRes,rRes,times;
 13 int s;
 14 
 15 int read()
 16 {
 17     char ch = getchar();
 18     int num = 0;
 19     bool fl = 0;
 20     for (; !isdigit(ch); ch = getchar())
 21         if (ch=='-') fl = 1;
 22     for (; isdigit(ch); ch = getchar())
 23         num = (num<<1)+(num<<3)+ch-48;
 24     if (fl) num = -num;
 25     return num;
 26 }
 27 inline bool limit(int x, int y){return x>=-4&&x<=4&&y>=-4&&y<=4;}
 28 inline bool illegal(int x)
 29 {
 30     if (x<-4) return (x<lRes&&lRes<=-4)||(x<rRes&&rRes<=-4);
 31     if (x>4) return (x>lRes&&lRes>=4)||(x>rRes&&rRes>=4);
 32     return 0;
 33 }
 34 void add(int x, int y, int c)
 35 {
 36     register int bar,cnt;
 37     if (y < -4||y > 4) return;
 38     if (s < x){
 39         if ((s<=lRes&&lRes<=x)||(s<=rRes&&rRes<=x)){
 40             if (s <= lRes){
 41                 bar = lRes, x = 2*bar-x-1;
 42             }else{
 43                 cnt = 1;
 44                 while (x <= lRes||x >= rRes)
 45                 {
 46                     if (cnt) bar = rRes;
 47                     else bar = lRes;
 48                     x = 2*bar-x, cnt = 1-cnt;
 49                     if (cnt) x++;
 50                     else x--;
 51                 }
 52             }
 53         }
 54     }
 55     if (s > x){
 56         if ((x<=lRes&&lRes<=s)||(x<=rRes&&rRes<=s)){
 57             if (s >= rRes){
 58                 bar = rRes, x = 2*bar-x+1;
 59             }else{
 60                 cnt = 1;
 61                 while (x <= lRes||x >= rRes)
 62                 {
 63                     if (cnt) bar = lRes;
 64                     else bar = rRes;
 65                     x = 2*bar-x, cnt = 1-cnt;
 66                     if (cnt) x--;
 67                     else x++;
 68                 }
 69             }
 70         }
 71     }
 72     if (limit(x, y)) f[x][y] += c;
 73 }
 74 void deal(int x, int y, int t)
 75 {
 76     if (t==0){
 77         if (limit(x, y)) f[x][y]++;
 78         return;
 79     }
 80     s = x;
 81     for (int delta=0; delta<=t; delta++)
 82     {
 83         add(x+delta, y+t-delta, 1), add(x+delta, y-t+delta, 1);
 84         add(x-delta, y+t-delta, 1), add(x-delta, y-t+delta, 1);
 85     }
 86     add(x, y+t, -1), add(x, y-t, -1);
 87     add(x+t, y, -1), add(x-t, y, -1);
 88     t -= 2;
 89     if (t==0){
 90         if (limit(x, y)) f[x][y]--;
 91         return;
 92     }else if (t < 0) return;
 93     for (int delta=0; delta<=t; delta++)
 94     {
 95         add(x+delta, y+t-delta, -1), add(x+delta, y-t+delta, -1);
 96         add(x-delta, y+t-delta, -1), add(x-delta, y-t+delta, -1);
 97     }
 98     add(x, y+t, 1), add(x, y-t, 1);
 99     add(x+t, y, 1), add(x-t, y, 1);
100 }
101 int main()
102 {
103     stones = read(), lRes = read(), rRes = read(), times = read();
104     if (lRes > rRes) std::swap(lRes, rRes);
105     for (int i=1; i<=stones; i++)
106     {
107         int x = read(), y = read(), t = read();
108         if (t > times||illegal(x)) continue;
109         deal(x, y, times-t);
110     }
111     for (int j=4; j>=-4; j--)
112     {
113         for (int i=-4; i<=4; i++)
114         {
115             if (i==lRes||i==rRes) putchar('X');
116             else if (f[i][j] < 0) putchar('o');
117             else if (f[i][j] > 0) putchar('*');
118             else putchar('-');
119         }
120         putchar('\n');
121     }
122     return 0;
123 }

后记

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

 

END

posted @ 2018-07-10 14:30  AntiQuality  阅读(212)  评论(0编辑  收藏  举报