Time:2016.05.08
Author:xiaoyimi
转载注明出处谢谢
传送门
思路:
说真的我一开始做这个题真不知道怎么建树,只知道应该以时间为根建树,记录优先值并维护前缀和,但想不出怎么用到主席树的前缀建树思想,暴力建树+修改肯定是不科学的……
后来去看题解,发现又是差分。
回想起被运输计划支配的恐惧……(还是自己太弱的缘故)
不过这个差分比较简单,因为它正好利用了第i棵主席树以第i-1棵树为基础建树的思路,比如某一任务的时间区间为[l,r],优先度为p,那么我们在root[l]上建树时+p,root[r+1]上建树时-p,这样的话就可以在枚举每个时间点的操作来建树了,因为有m个任务,每个任务有两个时间点的操作,所以最多建2m棵树,每个时间点所代表的主席树就是它最后一次修改所建的树,记录每个时间点最后被操作时所建的树的编号pos[i],查询操作就是空树与pos[i]相减得到了,这是主席树的基本操作,不过多赘述
注意:
1.简单计算可以得到这里的sum最多可以是100000*100000,超过了int,所以sum用long long计算
2.对于每个时间点的修改,有+有-,所以在原先的建树代码上稍作修改即可,记录系数什么的
代码:
#include<bits/stdc++.h>
#define M 100004
#define LL long long
#define ls(x) a[x].ch[0]
#define rs(x) a[x].ch[1]
using namespace std;
int n,m,cnt,k;
int ID[M],S[M],E[M],P[M],pos[M];
LL lastans=1;
vector <int> t[M];
struct disc
{
int data,id;
bool operator <(const disc other)const
{
return data<other.data;
}
}b[M];
struct Chairman_tree
{
int siz,ch[2];
LL sum;
}a[M<<6];
int in()
{
char ch=getchar();
int t=0;
while (!isdigit(ch)) ch=getchar();
while (isdigit(ch)) t=(t<<3)+(t<<1)+ch-48,ch=getchar();
return t;
}
void build(int now,int L,int R,int rt,int k,int val)
{
a[rt].siz=a[now].siz+k;
a[rt].sum=a[now].sum+k*ID[val];
if (L==R) return;
int mid=(L+R)>>1;
if (mid>=val)
rs(rt)=rs(now),
ls(rt)=++cnt,
build(ls(now),L,mid,ls(rt),k,val);
else
ls(rt)=ls(now),
rs(rt)=++cnt,
build(rs(now),mid+1,R,rs(rt),k,val);
}
LL get(int begin,int end,int L,int R,int k)
{
if (begin==end) return a[R].sum-a[L].sum;
int mid=(begin+end)>>1,t=a[ls(R)].siz-a[ls(L)].siz;
if (k<=t)
return get(begin,mid,ls(L),ls(R),k);
else
return a[ls(R)].sum-a[ls(L)].sum+get(mid+1,end,rs(L),rs(R),k-t);
}
main()
{
m=in();n=in();
for (int i=1;i<=m;i++)
S[i]=in(),E[i]=in(),P[i]=in(),
b[i].id=i,b[i].data=P[i];
sort(b+1,b+m+1);
for (int i=1;i<=m;i++)
ID[i]=b[i].data,
P[b[i].id]=i;
cnt=1;
for (int i=1;i<=m;i++)
t[S[i]].push_back(P[i]),
t[E[i]+1].push_back(-P[i]),
cnt+=1+(E[i]+1<=n);
for (int rt=1,i=1;i<=n;i++)
{
for (int j=0;j<t[i].size();j++)
build(rt,1,n,rt+1,t[i][j]/abs(t[i][j]),abs(t[i][j])),
rt++;
pos[i]=rt;
}
int x;
for (int i=1;i<=n;i++)
x=in(),
k=(lastans*in()+in())%in()+1,
printf("%lld\n",lastans=get(1,n,1,pos[x],k));
}