【NOI2008】 糖果雨

糖果雨
【问题描述】
有一个美丽的童话:在天空的尽头有一个"糖果国",这里大到摩天大厦,小
到小花小草都是用糖果建造而成的。更加神奇的是,天空中飘满了五颜六色的糖
果云,很快糖果雨密密麻麻从天而落,红色的是草莓糖,黄色的是柠檬糖,绿色
的是薄荷糖,
黑色的是巧克力糖......这时糖果国的小朋友们便会拿出大大小小的
口袋来接天空中落下的糖果,拿回去与朋友们一起分享。
对糖果情有独钟的小 Z 憧憬着能够来到这样一个童话的国度。所谓日有所
思,夜有所梦,这天晚上小 Z 梦见自己来到了"糖果国"。他惊喜地发现,任何时
候天空中所有的云朵颜色都不相同,
不同颜色的云朵在不断地落下相应颜色的糖
果。
更加有趣的是所有的云朵都在做着匀速往返运动,
不妨想象天空是有边界的,
而所有的云朵恰好在两个边界之间做着往返运动。
每一个单位时间云朵向左或向
右运动一个单位,当云朵的左界碰到天空的左界,它会改变方向向右运动;当云
朵完全移出了天空的右界,它会改变方向向左运动。
我们不妨把天空想象为一个平面直角坐标系,而云朵则抽象为线段(线段可
能退化为点):
如上图,不妨设天空的左界为 0,右界为 len。图中共有 5 片云朵,其中标
号为 1 的云朵恰好改变方向向右运动,标号为 2 的云朵恰好改变方向向左运动。
忽略云朵的纵坐标,它们在运动过程中不会相互影响。
小 Z 发现天空中会不断出现一些云朵(某个时刻从某个初始位置开始朝某个
方向运动),而有的云朵运动到一定时刻就会从天空中消失,而在运动的过程中
糖果在不断地下落。小 Z 决定拿很多口袋来接糖果,口袋容量是无限的,但袋
口大小却是有限的。例如在时刻 T 小 Z 拿一个横坐标范围为[L, R]的口袋来接糖
果,如果[L, R]存在一个位置 x,该位置有某种颜色的糖果落下,则认为该口袋可
接到此种颜色的糖果。极端情况下,袋口区间可能是一个点,譬如[0,0]、[1,1],
但仍然可以接到相应位置的糖果。通常可以接到的糖果总数会很大,因而小 Z
想知道每一次(即拿出口袋的一瞬间)他的口袋可以接到多少种不同颜色的糖果。
糖果下落的时间忽略不计。
【输入格式】
输入文件 candy.in 的第一行有两个正整数 n, len,
分别表示事件总数以及天空
的“边界”
。
接下来 n 行每行描述一个事件,
所有的事件按照输入顺序依次发生。
每行的
第一个数 k(k = 1,2,3)分别表示事件的类型,分别对应三种事件:插入事件,
询问事件以及删除事件。输入格式如下:
事件类型
输入格式
1 Ti Ci Li Ri Di                  插入事件                         (天空中出现了一片云朵)
2 Ti Li Ri                           询问事件        (询问一个口袋可以接到多少种不同颜色的糖果)
3 Ti Ci                               删除事件                        (天空中一片云朵消失了)
说明
时刻 Ti,天空中出现了一片坐标范围为[Li, Ri],颜色为 Ci 的云朵,初始的时候云朵运动方向
为向左(Di = -1)或向右(Di = 1)。满足 0 ≤ Li ≤ Ri ≤ len,Di = -1 或
1。
数据保证任何时刻空中不会出现两片颜色相同的云朵。时刻 Ti,小 Z 用一个坐标范围为[Li, Ri]的大口袋去接糖果,询问可以接到多少种不同的糖果。满足 0 ≤ Li ≤ Ri ≤ len。
时刻 Ti,颜色为 Ci 的云朵从天空消失中。数据保证当前天空中一定存在一片颜色为 Ci 的云朵。
【输出格式】
对于每一个询问事件,输出文件 candy.out 中应包含相应的一行,为该次询
问的答案,即口袋可以接到多少种不同的糖果。
【输入样例】
10 10
1 0 10 1 3 -1
2100
2 11 0 10
2 11 0 9
1 11 13 4 7 1
2 13 9 9
2 13 10 10
3 100 13
3 1999999999 10
1 2000000000 10 0 1 1
【输出样例】
1
1
0
2
1
【样例说明】
共 10 个事件,包括 3 个插入事件,5 个询问事件以及 2 个删除事件。
时刻 0,天空中出现一片颜色为 10 的云朵,初始位置为[1, 3],方向向左。
时刻 1,范围为[0, 0]的口袋可以接到颜色为 10 的糖果(云朵位置为[0, 2])。
时刻 11,范围为[0,10]的口袋可以接到颜色为 10 的糖果(云朵位置为[10, 12])。
时刻 11,范围为[0, 9]的口袋不能接到颜色为 10 的糖果(云朵位置为[10, 12])。
时刻 11, 天空中出现一片颜色为 13 的云朵, 初始位置为[4, 7], 方向向右。
时刻 13,范围为[9, 9]的口袋可以接到颜色为 10(云朵的位置为[8, 10])和颜色为 13(云朵的位置为[6, 9])两种不同的糖果。
时刻 13,范围为[10, 10]的口袋仅仅可以接到颜色为 10 的一种糖果(云朵的位
置为[8, 10]),而不可以接到颜色为 13 的糖果(云朵的位置为[6, 9]),。
时刻 100, 颜色为 13 的云朵从天空中消失。
时刻 1999999999,颜色为 10 的云朵从天空中消失。
时刻 2000000000,天空中又出现一片颜色为 10 的云朵,初始位置为[0, 1],
方向向右。
【数据规模和约定】
对于所有的数据,0 ≤ Ti ≤ 2000000000,1 ≤ Ci ≤ 1000000。
数据保证{Ti}为非递减序列即 T1 ≤ T2 ≤ ... ≤ Tn-1 ≤ Tn。
对于所有的插入事件,令 Pi = Ri – Li,即 Pi 表示每片云朵的长度。
数据编号           n                    len                  pi
1                       20                   10                 <=len
2                      200                 100                <=len
3                     2000               1000               <=len
4                   100000               10                 <=len
5                   100000              100                <=2
6                   150000             1000               <=3
7                   200000             1000               <=3
8                   100000             1000               <=len
9                   150000             1000               <=len
10                 200000             1000               <=len              

 

题解

 

(*
    *Problem:    NOI2008 糖果雨
    *Author :    Chen Yang
    *Time   :    2012.5.20
    *State  :    30分
    *Memo    :    模拟
*)
program candy;
const maxn=200200;
type
  ty1=record
    t,c,l,r,d:longint;
  end;
var
  n,len,tot,i,j,x,l,r,t,ans:longint;
  wind:array[0..maxn] of ty1;
  del:array[0..maxn] of boolean;
//=======================
function check1(var l,r:longint):boolean;
var
  t:longint;
begin
  t:=r-l;
  if r>len+t then
  begin
    r:=len+t-(r-len-t); l:=r-t;
    exit(true);
  end;
  exit(false);
end;
//=======================
function check2(var l,r:longint):boolean;
var
  t:longint;
begin
  t:=r-l;
  if l<0 then
  begin
    l:=0-l; r:=l+t;
    exit(true);
  end;
  exit(false);
end;
//=======================
function check(x:longint):boolean;
var
  lx,rx:longint;
  can:boolean;
begin
  lx:=wind[x].l+wind[x].d*((t-wind[x].t) mod (len shl 1));
  rx:=wind[x].r+wind[x].d*((t-wind[x].t) mod (len shl 1));
  //writeln(x,' ',lx,' ',rx,' ',wind[x].l,' ',wind[x].r);
  can:=true;
  while can do
  begin
    //writeln(x,' ',lx,' ',rx);
    can:=false;
    if check1(lx,rx) then can:=true;
    if check2(lx,rx) then can:=true;
  end;
  //writeln(t,' ',lx,' ',rx,' ',l,' ',r);
  exit(not ((lx>r)or(rx<l)));
end;
//=======================
begin
  assign(input,'candy.in'); reset(input);
  assign(output,'candy.out'); rewrite(output);
  read(n,len);
  for i:=1 to n do
  begin
    read(x);
    case x of
    1:begin
        inc(tot);
        read(wind[tot].t,wind[tot].c,wind[tot].l,wind[tot].r,wind[tot].d);
      end;
    2:begin
        read(t,l,r); ans:=0;
        for j:=1 to tot do
        if not del[j] and check(j) then inc(ans);
        writeln(ans);
      end;
    3:begin
        read(t,l);
        for j:=1 to tot do
        if wind[j].c=l then del[j]:=true;
      end;
    end;
  end;
  close(input); close(output);
end.

 

(*
        Problem:        NOI2008 糖果雨
    Author:        Chen Yang
    Time:        2012.5.25 8:29 pm
    State:        AC
    Memo:        二维树状数组、平行四边形变换统计
*)
program candy;
const maxc=1000100;
      maxm=4010;
      maxn=2010;
var
  n,len,m,i,p,t,l,r,c,d,q:longint;
  s:array[0..1,0..maxn,0..maxm] of longint;
  x,y:array[0..maxc] of longint;
//=================
procedure add(p,x,y,v:longint);
var
  i:longint;
begin
  inc(x); inc(y);
  while x<maxn do
  begin
    i:=y;
    while i<maxm do
    begin
      inc(s[p,x,i],v);
      inc(i,i and -i);
    end;
    inc(x,x and -x);
  end;
end;
//=================
function sum(p,x,y:longint):longint;
var
  i:longint;
begin
  if (x<0)or(y<0) then exit(0);
  inc(x); inc(y); sum:=0;
  if x>n then x:=n+1;
  if y>m then y:=m+1;
  while x>0 do
  begin
    i:=y;
    while i>0 do
    begin
      inc(sum,s[p,x,i]);
      dec(i,i and -i);
    end;
    dec(x,x and -x);
  end;
end;
//=================
procedure updata(t,c,l,r,d:longint);            inline;
begin
  x[c]:=(t-d*l+n) mod n; y[c]:=r-l;
  add(0,x[c],y[c]+x[c],1);
  add(1,x[c],y[c]-x[c]+n,1);
end;
//=================
procedure change(c:longint);                    inline;
begin
  add(0,x[c],y[c]+x[c],-1);
  add(1,x[c],y[c]-x[c]+n,-1);
end;
//=================
function area(p,x1,y1,x2,y2:longint):longint;   inline;
begin
  area:=sum(p,x2,y2)+sum(p,x1-1,y1-1)-sum(p,x1-1,y2)-sum(p,x2,y1-1);
end;
//=================
function ask(t,l,r:longint):longint;            inline;
var
  d:longint;
begin
  d:=longint(r=len);
  ask:=area(0,t,l+t,t+r,m)+area(0,0,l+t-n,t+r-n-d,m)+
       area(1,t-r+n+d,l-t,n,m)+area(1,t-r,l-t+n,t-1,m);
end;
//=================
begin
  assign(input,'candy.in'); reset(input);
  assign(output,'candy.out'); rewrite(output);
  read(q,len); n:=len << 1; m:=len << 2;
  for i:=1 to q do
  begin
    read(p);
    case p of
    1:begin read(t,c,l,r,d); updata(t,c,l,r,d);    end;
    2:begin read(t,l,r); writeln(ask(t mod n,l,r)); end;
    3:begin read(t,c); change(c); end;
    end;
  end;
  close(input); close(output);
end.

 

/*
        problem:        NOI2008 糖果雨
    Author:        Chen Yang
    Time:        2012.5.25 9:29 pm
    State:        AC
    Memo:        二维树状数组、平行四边形变换统计
*/
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <string>
#include <cstring>
using namespace std;
const int maxn=2020,maxm=4040,maxc=1000010;
int len,n,m,q;
int s[2][maxn][maxm],x[maxc],y[maxc];

void add(int p,int x,int y,int v)
{
    for (++x, ++y; x<maxn; x += x & -x)
    for (int i=y; i<maxm; i += i & -i) s[p][x][i] += v;
}

int sum(int p,int x,int y)
{
    if (x<0 || y<0) return 0; int k = 0;
    if (++x > n) x = n+1; if (++y > m) y = m+1;    
    for (; x; x -= x & -x)
    for (int i=y; i; i -= i & -i) k += s[p][x][i];
    return k;
}

inline void updata(int t,int c,int l,int r,int d)
{
    x[c] = (t - l*d + n) % n, y[c] = r - l;
    add(0,x[c],y[c]+x[c],1); add(1,x[c],y[c]-x[c]+n,1);
}

inline void del(int c)
{
    add(0,x[c],y[c]+x[c],-1); add(1,x[c],y[c]-x[c]+n,-1);
}

inline int area(int p,int x1,int y1,int x2,int y2)
{
    return sum(p,x2,y2)+sum(p,x1-1,y1-1)-sum(p,x1-1,y2)-sum(p,x2,y1-1);
}

inline int ask(int t,int l,int r)
{
    int d=r==len;
    return area(0,t,l+t,t+r,m)+area(0,0,l+t-n,t+r-n-d,m)+
           area(1,n-r+t+d,l-t,n,m)+area(1,t-r,l-t+n,t-1,m);
}

int main()
{
    freopen("candy.in","r",stdin);
    freopen("candy.out","w",stdout);
    scanf("%d%d", &q, &len); 
    n = len << 1; m = len << 2;
    for (int t,c,l,r,d,p; q--;)
    {
        scanf("%d%d", &p, &t);
        if (p==1)
        {
            scanf("%d%d%d%d", &c, &l, &r, &d);
            updata(t,c,l,r,d);
        } 
        else if (p==2)
        {
            scanf("%d%d", &l, &r);
            printf("%d\n",ask(t%n,l,r));
        } else
        {
            scanf("%d", &c);
            del(c);
        }
    }
}