[Ynoi2015]盼君勿忘

题目大意：

给定一个序列，每次查询一个区间\([l,r]\)中所有子序列分别去重后的和\(\bmod p\)（每次询问模数不同）。

解题思路：

在太阳西斜的这个世界里，置身天上之森。等这场战争结束之后，不归之人与望眼欲穿的众人， 人人本着正义之名，长存不灭的过去、逐渐消逝的未来。我回来了，纵使日薄西山，即便看不到未来，此时此刻的光辉，盼君勿忘。————世界上最幸福的女孩

珂朵莉，要永远幸福哦。

---

我们考虑每个数的贡献。即该区间内含有这个数的子序列个数。用补集转化为不含这个数的子序列个数。

那么，假设这个数在\([l,r]\)内出现了\(k\)次，则一共有\(2^{r-l+1}-2^{r-l+1-k}\)个子序列**包含**这个数。

所以，我们莫队，维护所有出现\(k\)次的不同数的和，然后包含这些数的子序列个数都相同，所以直接用和来计算贡献即可。

时间复杂度\(O(nm)\)，非常优秀。

考虑根号内分类讨论。

对于出现次数大于\(\sqrt n\)的数，这种数最多不到\(\sqrt n\)个。

所以我们莫队，统计出现次数小于等于\(\sqrt n\)的所有出现次数为\(i\)的和\(s_i\)，然后用一个东西来维护所有出现次数大于\(\sqrt n\)的数。

这样的话，每次询问的复杂度就是\(O(\sqrt n)\)了。

看上去好像没什么问题，实际上你会发现求2的幂次还有个\(\log\)，而且模数会改，不好预处理。

有一种分块\(O(\sqrt n)\)预处理，\(O(1)\)求幂次的方法，具体就是求出\(2^0,2^1,2^2,\dots,2^{{\sqrt n}-1}\)和\(2^{\sqrt n},2^{2\sqrt n},\dots,2^n\)（模意义），然后每个幂都分成两部分相乘即可。

每次询问重新计算一遍即可。

然后我们考虑用什么东西来维护出现次数大于\(\sqrt n\)的数。我用了unordered_set，理论上是线性的（当然这个理论有多可靠就不知道了）。

这样总时间复杂度\(O((n+m)\sqrt n)\)，空间复杂度\(O(n)\)。

C++ Code：

#include<cstdio>
#include<cctype>
#include<algorithm>
#include<vector>
#include<unordered_set>
#ifdef ONLINE_JUDGE
struct istream{
    char buf[23333333],*s;
    inline istream(){
        buf[fread(s=buf,1,23333330,stdin)]='\n';
        fclose(stdin);
    }
    inline istream&operator>>(int&d){
        d=0;
        for(;!isdigit(*s);++s);
        while(isdigit(*s))
        d=(d<<3)+(d<<1)+(*s++^'0');
        return*this;
    }
}cin;
struct ostream{
    char buf[8000005],*s;
    inline ostream(){s=buf;}
    inline ostream&operator<<(int d){
        if(!d){
            *s++='0';
        }else{
            static int w;
            for(w=1;w<=d;w*=10);
            for(;w/=10;d%=w)*s++=d/w^'0';
        }
        return*this;
    }
    inline ostream&operator<<(const char&c){*s++=c;return*this;}
    inline void flush(){
        fwrite(buf,1,s-buf,stdout);
        s=buf;
    }
    inline~ostream(){flush();}
}cout;
#else
#include<iostream>
using std::cin;
using std::cout;
#endif
#define siz 317
#define N 100005
int n,m,a[N],buc[N],_2[siz+2],__2[siz+2],out[N];
long long bvc[siz+2];
std::vector<int>lr;
std::unordered_set<int>s;
inline int pow2(int b,const int&md){
	int ret=1,a=2;
	for(;b;b>>=1,a=1LL*a*a%md)
	if(b&1)ret=1LL*ret*a%md;
	return ret;
}
inline int pw2(int m,const int&md){
	return 1LL*_2[m%siz]*__2[m/siz]%md;
}
struct que{
	int l,r,id,md;
	inline bool operator<(const que&rhs)const{
		return(l/siz!=rhs.l/siz)?(l<rhs.l):(r<rhs.r);
	}
}q[N];
inline void add(int pos){
	const int val=a[pos];
	if(buc[val]>siz)++buc[val];else
	if(buc[val]==siz)bvc[buc[val]++]-=lr[val],s.insert(val);else
	bvc[buc[val]]-=lr[val],bvc[++buc[val]]+=lr[val];
}
inline void del(int pos){
	const int val=a[pos];
	if(buc[val]>siz+1)--buc[val];else
	if(buc[val]==siz+1)s.erase(val),bvc[--buc[val]]+=lr[val];else
	bvc[buc[val]]-=lr[val],bvc[--buc[val]]+=lr[val];
}
int main(){
	cin>>n>>m;
	for(int i=1;i<=n;lr.push_back(a[i++]))cin>>a[i];
	std::sort(lr.begin(),lr.end());
	for(int i=1;i<=n;++i)a[i]=lower_bound(lr.begin(),lr.end(),a[i])-lr.begin();
	for(int i=1;i<=m;++i)
	cin>>q[q[i].id=i].l>>q[i].r>>q[i].md;
	std::sort(q+1,q+m+1);
	int l=1,r=0;
	for(int i=1;i<=n;++i)*bvc+=a[i];
	*_2=*__2=1;
	for(int i=1;i<=m;++i){
		while(r<q[i].r)add(++r);
		while(l>q[i].l)add(--l);
		while(r>q[i].r)del(r--);
		while(l<q[i].l)del(l++);
		const int md=q[i].md,len=r-l+1;
		for(int j=1;j<siz;++j)
		_2[j]=2LL*_2[j-1]%md;
		__2[1]=pow2(siz,md);
		for(int j=2;j<=siz;++j)
		__2[j]=__2[1]*1LL*__2[j-1]%md;
		int&ans=out[q[i].id];
		const int all=pw2(len,md);
		for(int j=1;j<=siz;++j)
		(ans+=1LL*(all-pw2(len-j,md)+md)*bvc[j]%md-md)+=ans>>31&md;
		for(int j:s){
			(ans+=1LL*(all-pw2(len-buc[j],md)+md)*lr[j]%md-md)+=ans>>31&md;
		}
	}
	for(int i=1;i<=m;++i)cout<<out[i]<<'\n';
	return 0;
}