寒假训练1.18

普通莫队应用题
P2709

#include<bits/stdc++.h>
using namespace std;
const int N=5e4+40;
typedef long long ll;
int col[N];
ll num[N]; 
ll res,ans[N],b[N];
struct node{
	int l,r,i;
	bool operator<(const node& a)const {
		return (b[l]!=b[a.l])?(b[l]<b[a.l]):r<a.r;
	}
}q[N];
void add(int x,int k){
	res-=num[k]*num[k];
	num[k]++;
	res+=num[k]*num[k];
}
void del(int x,ll k){
	res-=num[k]*num[k];
	num[k]--;
	res+=num[k]*num[k];
}
int main(){
	int n,m,k;
	scanf("%d%d%d",&n,&m,&k);
	for(int i=1;i<=n;i++){
		scanf("%d",&col[i]);
	}
	int len=sqrt(n);
	for(int i=1;i<=n;i++){
		b[i]=i/len;
	}
	for(int i=0;i<m;i++){
		scanf("%d%d",&q[i].l,&q[i].r);
		q[i].i=i;
	}
	sort(q,q+m);
	int l=1,r=0;
	for(int i=0;i<m;i++){
		int x=q[i].l,y=q[i].r;
		while(x<l){
			--l;
			add(l,col[l]);
		}
		while(y>r){
			++r;
			add(r,col[r]);
		}
		while(x>l){
			del(l,col[l]);
			++l;
		}
		while(y<r){
			del(r,col[r]);
			--r;
		}
		ans[q[i].i]=res;
	}
	for(int i=0;i<m;i++){
		printf("%lld\n",ans[i]);
	}
	return 0;
}

想到了差分然后想到了方程就是没有想到bitset优化
结合我们队之前17年的ccpc杭州吧(记得好像是)有个点分治+bitset优化也是求的是相加起来是x的数

//差分加上莫队的存在
//bitset优化莫队的存在
//bitset的好处以及用法 
#include<bits/stdc++.h>
using namespace std;
const int N=1e5+50;
int num[N],ans[N],col[N];//一个是记录出正的个数恰好就等于负的个数 
int b[N];
struct node{
	int l,r,s,i,x;
	bool operator<(const node& a)const {
		return (b[l]!=b[a.l])?(b[l]<b[a.l]):r<a.r;
	}
}q[N]; 
bitset<N>aa,bb,cc,dd;
int main(){
	int n,m;
	scanf("%d%d",&n,&m);
	int len=sqrt(n);
	for(int i=1;i<=n;i++){
		scanf("%d",&col[i]); 
		b[i]=i/len;
	}
	for(int i=0;i<m;i++){
		scanf("%d%d%d%d",&q[i].s,&q[i].l,&q[i].r,&q[i].x);
		q[i].i=i;
	}
	sort(q,q+m);
	int s;
	int l=1,r=0;
	for(int i=0;i<m;i++){
		int x=q[i].l,y=q[i].r;
		while(x<l){
			--l;
			aa[col[l]]=1;//表示的是存在 
			dd[N-col[l]]=1;
			num[col[l]]++;
		}
		while(y>r){
			++r;
			num[col[r]]++;
			aa[col[r]]=1;
			dd[N-col[r]]=1;
		}
		while(x>l){
			num[col[l]]--;
			if(!num[col[l]]){
				aa[col[l]]=0;
				dd[N-col[l]]=0;
			}
			l++;
		}
		while(y<r){
			num[col[r]]--;
			if(!num[col[r]]){
				aa[col[r]]=0;
				dd[N-col[r]]=0;
			}
			r--;
		}
//		cout<<x<<" "<<y<<" :"<<endl;
//		for(int i=0;i<N;i++){
//			if(aa[i]){
//				cout<<i<<endl;
//			}
//		}
		s=q[i].s;
		int d=q[i].x; 
		if(s==1){
			bb=aa<<d;
			cc=bb&aa;
			if(cc.count()>0){
				ans[q[i].i]=1; 
			} 
		}else if(s==2){
			bb=dd>>(N-d);
			cc=bb&aa;
			if(cc.count()>0){
//				for(int i=0;i<N;i++){
//					if(aa[i]){
//						cout<<"aa"<<" "<<i<<endl;
//					}
//					if(dd[i]){
//						cout<<"dd "<<N-i<<" "<<i<<endl;
//					}
//					if(cc[i]){
//						cout<<x<<" "<<y<<" "<<d<<endl;
//						cout<<i<<endl;
//					}
//				}
				ans[q[i].i]=1;
			}
		}else {
			//这个的话就是直接进行枚举了
			for(int j=1;j*j<=d;j++){
				if(d%j)	continue;
				if(aa[j]&&aa[d/j]){
					ans[q[i].i]=1;
					break;
				}
			} 
		}
	}
	for(int i=0;i<m;i++){
		if(ans[i]){
			printf("hana\n");
		}else printf("bi\n");
	}
	return 0;
}

//树状数组+二次离散莫队

#include <bits/stdc++.h>
#define lowbit(x) (x&-x)
#define ll long long
using namespace std;
const int maxn=500005,maxa=100005,maxsa=325;
int n,m,a[maxn],s,ss,idx[maxn],Maxa,idxn,l=1,r,ccnt,idxx[maxa];
ll qz[maxn],ans[maxn],aans[maxn],sum1[maxa],sum2[maxsa],sum3[maxa],sum4[maxsa],qqz[maxn];
struct query{//询问 
	int l,r,id;
}p[maxn];
struct node{//二离 
	int l,r,op,id,x;
}que[maxn<<1];
int qr(){
	int x=0,f=0;
	char c=getchar();
	while(c>'9'||c<'0') f|=(c=='-'),c=getchar();
	while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
	return f?-x:x;
}
void add(int x){//O(sqrt(n))的修改 
	for(int i=x;idxx[i]==idxx[x];++i) ++sum1[i],sum3[i]+=x;
	for(int i=idxx[x];i<=idxx[Maxa];++i) ++sum2[i],sum4[i]+=x;
}
inline bool cmp(query x,query y){//奇偶排序 
	return idx[x.l]==idx[y.l]?(idx[x.l]&1?x.r<y.r:x.r>y.r):x.l<y.l;
}
inline bool cmp2(node x,node y){
	return x.x<y.x;
}
inline ll getsum(int x){
	return sum4[idxx[Maxa]]-sum4[idxx[x]-1]-sum3[x]+1ll*x*(sum2[idxx[x]-1]+(idxx[x-1]==idxx[x]?sum1[x-1]:0));//O(1) 查询 
}
int main() {
	n=qr(),m=qr();
	s=sqrt(n);
	for(int i=1;i<=n;++i){
		a[i]=qr();
		idx[i]=(i+s-1)/s;
		qqz[i]=qqz[i-1]+a[i];//最后要用 
		Maxa=max(Maxa,a[i]);//值域大小 
	}
	ss=sqrt(Maxa);
	for(int i=1;i<=Maxa;++i){//初始化值域块 
		idxx[i]=(i+ss-1)/ss;
	}
	for(int i=1;i<=n;++i){//求前缀贡献 
		qz[i]=qz[i-1]+getsum(a[i]);
		add(a[i]);
	}
	for(int i=1;i<=m;++i){
		p[i].l=qr(),p[i].r=qr(),p[i].id=i;
	}
	sort(p+1,p+m+1,cmp);//奇偶排序 
	for(int i=1;i<=m;++i){//二次离线转移部分，记得离线整体(打的 if) 
		if(l>p[i].l){
			que[++ccnt]=(node){p[i].l,l-1,1,i,r};
			ans[i]-=qz[l-1]-qz[p[i].l-1];
			l=p[i].l;
		}
		if(r<p[i].r){
			que[++ccnt]=(node){r+1,p[i].r,-1,i,l-1};
			ans[i]+=qz[p[i].r]-qz[r];
			r=p[i].r;
		}
		if(l<p[i].l){
			que[++ccnt]=(node){l,p[i].l-1,-1,i,r};
			ans[i]+=qz[p[i].l-1]-qz[l-1];
			l=p[i].l;
		}
		if(r>p[i].r){
			que[++ccnt]=(node){p[i].r+1,r,1,i,l-1};
			ans[i]-=qz[r]-qz[p[i].r];
			r=p[i].r;
		}
	}
	sort(que+1,que+ccnt+1,cmp2);//按那个端点排序 
	memset(sum1,0,sizeof(sum1));
	memset(sum2,0,sizeof(sum2));
	memset(sum3,0,sizeof(sum3));
	memset(sum4,0,sizeof(sum4));
	for(int i=1;i<=ccnt;++i){
		if(que[i].x==0) continue;//不用减没有贡献的 
		while(idxn<que[i].x) add(a[++idxn]);//转移 [1,l-1] 的 l-1 
		for(int j=que[i].l;j<=que[i].r;++j) ans[que[i].id]+=1ll*que[i].op*getsum(a[j]);//枚举，前缀和在前面。 
	}
	for(int i=2;i<=m;++i) ans[i]+=ans[i-1];//因为处理的只有转移的部分，所以要前缀和。 
	for(int i=1;i<=m;++i) aans[p[i].id]=ans[i]+qqz[p[i].r]-qqz[p[i].l-1];//先要全部加一遍，前面也可以处理。 
	for(int i=1;i<=m;++i) printf("%lld\n",aans[i]);
	return 0;
}                  

n根号mlogV复杂度写法

//这个是个模板题
//维护一个树状数组
//二次离散莫队 
#include<bits/stdc++.h>
using namespace std;
const int N=5e5+40;
typedef long long ll;
ll c[N],d[N];//需要的是两个
ll col[N],b[N],res,ans[N];
int n,m;
ll ma;
int lowbit(int x) {
	return x&(-x);
}
void add(int x,ll k,ll c[]) { //在某个位置之前加上k
	while(x<=ma) {
		c[x]+=k;
		x+=lowbit(x);
	}
}
ll get(int x,ll c[]) {
	ll res=0;
	while(x) {
		res+=c[x];
		x-=lowbit(x);
	}
	return res;
}

int read() {
	char ch=getchar();
	int x=0,cf=1;
	while(ch<'0'||ch>'9') {
		if(ch=='-') cf=-1;
		ch=getchar();
	}
	while(ch>='0'&&ch<='9') {
		x=(x<<3)+(x<<1)+(ch^48);
		ch=getchar();
	}
	return x*cf;
}


struct node {
	int l,r,i;
	bool operator<(const node& a)const {
		return (b[l]<b[a.l])?(b[l]<b[a.l]):r<a.r;
	}
} q[N];
int main() {
	n=read(),m=read();
	int len=sqrt(n);
	for(int i=1; i<=n; i++) {
		col[i]=read();
		ma=max(ma,col[i]);
		b[i]=i/len;
	}
	for(int i=0; i<m; i++) {
		scanf("%d%d",&q[i].l,&q[i].r);
		q[i].i=i;
	}
	sort(q,q+m);
	int l=1,r=0;
	for(int i=0; i<m; i++) {
		int x=q[i].l,y=q[i].r;
		ll sum,temp;
		//添加的
		while(x<l) {
			--l;
			sum=get(ma,d);//这个是总的答案
			temp=get(col[l],d);
			//第一部分是对于原来的变化还有一个新加的点
			res+=sum-temp+(get(col[l]-1,c)+1ll)*col[l];//还有一部分的话是现在
			add(col[l],col[l],d);
			add(col[l],1,c);//对于这个的话是数量上的差分
		}//直接开始进行一个添加的操作
		while(y>r) {
			++r;
			sum=get(ma,d);
			temp=get(col[r],d);
			res+=sum-temp+(get(col[r]-1,c)+1)*col[r];
			add(col[r],col[r],d);
			add(col[r],1,c);
		}
		while(x>l) {
			sum=get(ma,d),temp=get(col[l],d);
			res-=sum-temp+(get(col[l]-1,c)+1)*col[l];
			add(col[l],-col[l],d);
			add(col[l],-1,c);
			l++;
		}
		while(y<r) {
			sum=get(ma,d),temp=get(col[r],d);
			res-=sum-temp+(get(col[r]-1,c)+1)*col[r];
			add(col[r],-col[r],d);
			add(col[r],-1,c);
			--r;
		}
		ans[q[i].i]=res;
	}
	for(int i=0; i<m; i++) {
		printf("%lld\n",ans[i]);
	}
	return 0;
}

//还有一个题目的话是考虑出最少能够分成几个严格上升子序列

#include <cstdio>
#include <iostream>
#include <cmath>
#include <algorithm>
#include <vector>

using namespace std;

const int N=2e5+10;

int n,m,a[N],cnt[N],len,bl[N],ans[N],maxn,sum[N];//cnt[i]:数字i出现的个数；sum[i]:出现次数为i的数字个数
vector <int> v; 
struct query
{
	int l,r,id;
	bool operator <(const query x) const {return bl[l]^bl[x.l]? l<x.l:bl[l]&1? r<x.r:r>x.r;}
} q[N];

inline void add(int x)
{
	sum[++cnt[a[x]]]++;
	if(cnt[a[x]]>maxn) maxn=cnt[a[x]];
	return ;
}
inline void del(int x)
{
	if(sum[cnt[a[x]]] == 1 && maxn == cnt[a[x]]) maxn--;
	sum[cnt[a[x]]--]--;
	return ;
}

int main()
{
// 	freopen("work.in","r",stdin); freopen("work.out","w",stdout);
	scanf("%d%d",&n,&m);
	len=sqrt(n);
	for(int i=1;i<=n;i++) {scanf("%d",&a[i]); v.push_back(a[i]); bl[i]=(i-1)/len+1;}
	sort(v.begin(),v.end());
	v.erase(unique(v.begin(),v.end()),v.end());//离散化
	for(int i=1;i<=n;i++) a[i]=lower_bound(v.begin(),v.end(),a[i])-v.begin()+1;
	for(int i=1;i<=m;i++) {scanf("%d%d",&q[i].l,&q[i].r); q[i].id=i;}
	sort(q+1,q+1+m);
	int l=1,r=0;
	for(int i=1;i<=m;i++)
	{
		while(l > q[i].l) add(--l);
		while(r < q[i].r) add(++r);
		while(r > q[i].r) del(r--);
		while(l < q[i].l) del(l++);
		ans[q[i].id]=-maxn;
	}
	for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
// 	fclose(stdin); fclose(stdout);
	return 0;
}