Joy of Handcraft（线段树）

题意

有$n$个灯泡，每个灯泡有两个参数$t_i$和$x_i$，分别指运行周期和亮度。第$i$个灯泡在$2k{t}_i+1$秒到$2k{t}_i+t_i$秒是亮的，在$2k{t}_i+t_i+1$秒到$2k{t}_i+2t_i$秒是灭的。问第$1$秒到第$m$秒，每一秒亮度最大值是多少。

数据范围

$1\le T\le 100$
$1\le n,m\le 10^5$
$1\le t_i,x_i\le 10^5$

思路

其实这道题的关键是对区间数量的估计。这里有一个性质，如果两个灯泡的$t_i$相等，那么我们就可以忽略掉$x_i$较小的那个灯泡。因此所有的区间数量就是${\sum }_{i=1}^n\frac{m}{t_i}=mlogn$。所以，我们可以预处理出这些区间。
然后就是个区间染色问题，用线段树维护即可。

代码

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;

typedef pair<int,int> pii;

const int N = 100010;

int n, m;
pii bulbs[N];
bool is_ap[N];
struct Line
{
    int st, ed, w;
}line[20*N];

struct Node
{
    int l, r;
    int color, flag;
}tr[N * 4];

void pushup(int u)
{
    if(tr[u<<1].color==tr[u<<1|1].color) tr[u].color = tr[u<<1].color;
    else tr[u].color = 0;
}

void pushdown(int u)
{
    auto &root = tr[u], &left = tr[u<<1], &right = tr[u<<1|1];
    if(root.flag){
        left.flag = root.flag, left.color = root.flag;
        right.flag = root.flag, right.color = root.flag;
        root.flag = 0;
    }
}

void build(int u,int l,int r)
{
    tr[u] = {l,r,0,0};
    if(l==r) return;
    int mid = l + r >> 1;
    build(u<<1,l,mid), build(u<<1|1,mid+1,r);
    pushup(u);
}

void modify(int u,int l,int r,int c)
{
    if(tr[u].l>=l&&tr[u].r<=r){
        tr[u].flag = c;
        tr[u].color = c;
        return;
    }
    pushdown(u);
    int mid = tr[u].l + tr[u].r >> 1;
    if(l<=mid) modify(u<<1,l,r,c);
    if(r>mid) modify(u<<1|1,l,r,c);
    pushup(u);
}

int query(int u,int l,int r)
{
    if(tr[u].l>=l&&tr[u].r<=r) return tr[u].color;
    pushdown(u);
    int mid = tr[u].l + tr[u].r >> 1;
    if(l<=mid) return query(u<<1,l,r);
    else return query(u<<1|1,l,r);
}

int main()
{
    int T;
    scanf("%d",&T);
    for(int cas=1;cas<=T;cas++){
        scanf("%d%d",&n,&m);
        memset(is_ap,0,sizeof(is_ap));
        int cnt = 1;
        for(int i=1;i<=n;i++) scanf("%d%d",&bulbs[i].second,&bulbs[i].first);
        sort(bulbs+1,bulbs+n+1);
        for(int i=n;i>=1;i--){
            int x = bulbs[i].first, t = bulbs[i].second;
            if(is_ap[t]) continue;
            //for(int k=2;k<t/k;k++) is_ap[t/k] = true;
            is_ap[t] = true;
            for(int k=0;;k++){
                int st = 2*k*t+1, ed = 2*k*t+t;
                if(st>m) break;
                if(ed>m) ed = m;
                line[cnt++] = {st,ed,x};
            }
        }
        build(1,1,m);
        printf("Case #%d:",cas);
        for(int i=cnt-1;i>0;i--){
            int st = line[i].st, ed = line[i].ed, w = line[i].w;
            modify(1,st,ed,w);
        }
        for(int i=1;i<=m;i++){
            int ans = query(1,i,i);
            printf(" %d",ans);
        }
        printf("\n");
        //for(int i=1;i<=m*4;i++) tr[i] = {0,0,0,0};
    }
    return 0;
}
/*
3
2 3
1 1
2 2
2 5
1 2
2 3
3 3
1 1
1 2
1 3
*/