BZOJ 3781 小B的询问 莫队
2017-08-05 17:07
399 查看
3781: 小B的询问
Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 937 Solved: 632
[Submit][Status][Discuss]
Description
小B有一个序列,包含N个1~K之间的整数。他一共有M个询问,每个询问给定一个区间[L..R],求Sigma(c(i)^2)的值,其中i的值从1到K,其中c(i)表示数字i在[L..R]中的重复次数。小B请你帮助他回答询问。Input
第一行,三个整数N、M、K。第二行,N个整数,表示小B的序列。
接下来的M行,每行两个整数L、R。
Output
M行,每行一个整数,其中第i行的整数表示第i个询问的答案。Sample Input
6 4 31 3 2 1 1 3
1 4
2 6
3 5
5 6
Sample Output
69
5
2
HINT
对于全部的数据,1<=N、M、K<=50000题解:
今天做水题了吗?做了。
第一道莫队的题23333巨水无比。
莫队算法就是先将所有的询问都记录下来,按照区间排序,然后每次根据询问的区间不同,在左端点和右端点加加减减。
这道题就是先减去要修改的端点处的值的平方,然后修改和,再加上修改后的和的平方。
(上午瞎搞了一下fread读入,无耻地copy了大佬的代码)
/************************************************************** Problem: 3781 User: LaLaLa112138 Language: C++ Result: Accepted Time:1272 ms Memory:24444 kb ****************************************************************/ #include<cstdio> #include<cstring> #include<cmath> #include<ctime> #include<algorithm> #define ll long long #define N 100010 using namespace std; namespace readIn{ const int maxn = 20*1024*1024; char buf[maxn]; char *s,*t; inline char Getchar(){ if(s==t){ t=(s=buf)+fread(buf,1,maxn,stdin); if(s==t) return EOF; } return *s++; } inline int read(){ int x=0,f=1;char ch=Getchar(); while(ch<'0'||ch>'9') {if(ch=='-') f=-1;ch=Getchar();} while(ch>='0'&&ch<='9') {x=x*10+ch-'0';ch=Getchar();} return x*f; } } int n,m,k; int a ,c ; int blo; ll ans ; struct node{ int l,r,id; int bl; }q ; bool cmp(const node &a,const node &b){ return (a.bl==b.bl&&a.r<b.r)||a.bl<b.bl; } ll cnt=0; void modify(int x,int d){ cnt-=c[a[x]]*c[a[x]]; c[a[x]]+=d; cnt+=c[a[x]]*c[a[x]]; } int main(){ //freopen("test.txt","r",stdin); n=readIn::read(),m=readIn::read(),k=readIn::read(); blo=sqrt(n); for(register int i=1;i<=n;i++) a[i]=readIn::read(); for(register int i=1;i<=m;i++) q[i].l=readIn::read(),q[i].r=readIn::read(),q[i].id=i,q[i].bl=(q[i].l-1)/blo+1; sort(q+1,q+m+1,cmp); int l=1,r=0; for(register int i=1;i<=m;i++){ while(l>q[i].l) {l--;modify(l,1);} while(r<q[i].r) {r++;modify(r,1);} while(l<q[i].l) {modify(l,-1);l++;} while(r>q[i].r) {modify(r,-1);r--;} ans[q[i].id]=cnt; } for(register int i=1;i<=m;i++){ printf("%lld\n",ans[i]); } return 0; }
相关文章推荐
- [bzoj] 3781 小B的询问 || 莫队
- BZOJ[3781]小B的询问 莫队
- [bzoj3781]小B的询问 智障莫队
- bzoj3781 小B的询问(莫队)
- bzoj3781: 小B的询问(莫队)
- 【bzoj3781】小B的询问 莫队算法
- 【BZOj 3781】小B的询问 莫队
- bzoj 3781: 小B的询问 莫队算法+分块
- BZOJ 3781 小B的询问 序列莫队算法
- BZOJ 3781: 小B的询问 莫队算法
- bzoj3781 小B的询问 莫队
- 【模板】BZOJ 3781: 小B的询问 莫队算法
- BZOJ_3781_小B的询问_莫队
- 【BZOJ3781】小B的询问 莫队
- bzoj 3781 小B的询问(莫队算法)
- [BZOJ3781]小B的询问(莫队)
- [bzoj3781]小B的询问_莫队
- 【BZOJ】3781 小B的询问 莫队
- BZOJ 3781: 小B的询问 [莫队]
- 洛谷P2709 BZOJ 3781 小B的询问 (莫队)