正解:整体二分

解题报告:

传送门! 还有个双倍经验!(明明是一样的题目为什么你们一个紫一个黑啊喂!

这题首先要想到可以二分嘛,然后看到多组询问肯定就整体二分鸭

那就是基本套路啊,发现是区间修改单点查询,于是就树状数组前缀和维护一波,就差不多辣

其实是个板子题,只是刚学整体二分所以还是写下题解QwQ

然后有几个细节分别港下

第一个是注意到它是个环形的,所以在update的时候记得分情况讨论一下

第二个是关于判断NIE,可以在最后增加一次流星雨,保证这次流星雨之后一定所有国家都能得到足够的流星,输出的时候判断一下就好

最后一个是比较重要的,我发现好几篇题解都被hack辣,,,就是第11个点,如果WA辣,应该就是这个问题,大概港下

首先从数据范围可以发现其实它是可以溢出的,就是假如有3e5场流星雨,每场提供1e9的流星,那在查询的时候就会炸longlong

但是看到它的需求也是<=1e9的,所以可以在算的时候判断一下,如果得到的已经大于1e9就可以直接break辣,只要加一句话就好QwQ

#include<bits/stdc++.h>
using namespace std;
#define il inline
#define rg register
#define gc getchar()
#define ll long long
#define lowbit(x) (x&(-x))
#define rp(i,x,y) for(rg ll i=x;i<=y;++i)

const ll N=3e5+10,inf=1e9+10;
ll n,m,q,as[N],bst[N];
vector<ll>bl[N];
struct node{ll l,r,dat;}quest[N];
struct nd{ll nd,id;}p[N],tmp_l[N],tmp_r[N];

il ll read()
{
    rg char ch=gc;rg ll x=0;rg bool y=1;
    while(ch!='-' && (ch>'9' || ch<'0'))ch=gc;
    if(ch=='-')ch=gc,y=0;
    while(ch>='0' && ch<='9')x=(x<<1)+(x<<3)+(ch^'0'),ch=gc;
    return y?x:-x;
}
il void updat(ll x,ll dat){while(x<=m)bst[x]+=dat,x+=lowbit(x);}
il ll query(ll x){ll ret=0;while(x)ret+=bst[x],x-=lowbit(x);return ret;}
void solv(ll ct_l,ll ct_r,ll as_l,ll as_r)
{
    if(ct_l>ct_r)return;
    if(as_l==as_r){rp(i,ct_l,ct_r)as[p[i].id]=as_l;return;}
    ll mid=(as_l+as_r)>>1,cnt_l=0,cnt_r=0;
    rp(i,as_l,mid)
        if(quest[i].l<=quest[i].r)updat(quest[i].l,quest[i].dat),updat(quest[i].r+1,-quest[i].dat);
        else updat(1,quest[i].dat),updat(quest[i].r+1,-quest[i].dat),updat(quest[i].l,quest[i].dat);
    rp(i,ct_l,ct_r)
    {
        ll ret=0,sz=bl[p[i].id].size();
        rp(j,0,sz-1){ret+=query(bl[p[i].id][j]);if(ret>inf)break;}
        if(ret>=p[i].nd)tmp_l[++cnt_l]=p[i];else p[i].nd-=ret,tmp_r[++cnt_r]=p[i];
    }
    rp(i,as_l,mid)
        if(quest[i].l<=quest[i].r)updat(quest[i].l,-quest[i].dat),updat(quest[i].r+1,quest[i].dat);
        else updat(1,-quest[i].dat),updat(quest[i].r+1,quest[i].dat),updat(quest[i].l,-quest[i].dat);
    rp(i,ct_l,ct_l+cnt_l-1)p[i]=tmp_l[i-ct_l+1];rp(i,ct_l+cnt_l,ct_r)p[i]=tmp_r[i-cnt_l-ct_l+1];
    solv(ct_l,ct_l+cnt_l-1,as_l,mid);solv(ct_l+cnt_l,ct_r,mid+1,as_r);
}

int main()
{
    n=read();m=read();rp(i,1,m)bl[read()].push_back(i);rp(i,1,n)p[i]=(nd){read(),i};q=read();rp(i,1,q)quest[i].l=read(),quest[i].r=read(),quest[i].dat=read();quest[q+1]=(node){1,m,1e9};
    solv(1,n,1,q+1);rp(i,1,n)if(as[i]!=q+1)printf("%lld\n",as[i]);else printf("NIE\n");
    return 0;
}

View Code

洛谷P3527 MET-Meteors [POI2011] 整体二分

正解:整体二分

解题报告:

相关推荐