BZOJ_2038

第一次写莫队了，之前偷懒一直没有写，听说还有一道历史观察啥的也可以做做，反正蛮水的。

但还是好好说一说。首先根号N分块，然后判断左端点在哪个块里面，按左端点的块为第一关键词，r为第二关键词排序。这样可以保证每一块里面r都是递增的。首先我们保证了r的复杂度是n根号哒，然后再是左端点，没跨过块的是根号，否则是暴力重建，还是根号哒。

所以总复杂度是n根号

 1 #include <cstdio>
 2 #include <iostream>
 3 #include <cstdlib>
 4 #include <algorithm>
 5 #include <cstring>
 6 #include <cmath>
 7 using namespace std;
 8 const int N = 50003;
 9 struct qus {
10     int l,r,b,key;
11 }query[N],ans[N];
12 int size,a[N],n,m,sta[N];
13 long long fenzi,fenmu,gcdd;
14 inline bool cmp(qus x,qus y) {
15     if(x.b < y.b) return true;
16     if(x.b > y.b) return false;
17     return x.r < y.r;
18 }
19 inline long long gcd(long long a,long long b) { return b==0?a:gcd(b,a%b); }
20 int main() {
21     scanf("%d%d",&n,&m);
22     size = sqrt(n);
23     for(int i = 1 ; i <= n ; ++i) scanf("%d",&a[i]);
24     for(int i = 1 ; i <= m ; ++i) {
25         scanf("%d%d",&query[i].l,&query[i].r);
26         query[i].b = query[i].l/size;
27         query[i].key = i;
28         if(query[i].l%size!=0) query[i].b++;
29     }
30     sort(query+1,query+1+m,cmp);
31     for(int i = 1 ; i <= m ; ++i) {
32         if(i==1 || query[i].b!=query[i-1].b) {
33             memset(sta,0,sizeof(sta));
34             fenzi = 0,fenmu = (query[i].r - query[i].l + 1) * (query[i].r - query[i].l);
35             for(int j = query[i].l ; j <= query[i].r ; ++j) sta[a[j]]++;
36             for(int j = 1 ; j <= n ; ++j) fenzi += sta[j]*(sta[j]-1);
37             gcdd = gcd(fenzi,fenmu);
38         } else {
39             fenmu = (query[i].r - query[i].l + 1) * (query[i].r - query[i].l);
40             for(int j = query[i-1].r+1; j <= query[i].r ; ++j) 
41             fenzi += sta[a[j]]*2,sta[a[j]]++;
42             for(int j = query[i-1].l ; j < query[i].l ; ++j) 
43             fenzi -= (sta[a[j]]-1)*2,sta[a[j]]--;
44             for(int j = query[i].l ; j < query[i-1].l ; ++j) 
45             fenzi += sta[a[j]]*2,sta[a[j]]++;
46             gcdd = gcd(fenzi,fenmu);
47         }
48         ans[query[i].key].l = fenzi/gcdd;
49         ans[query[i].key].r = fenmu/gcdd;
50     }
51     for(int i = 1 ; i <= m ; ++i) printf("%d/%d\n",ans[i].l,ans[i].r);
52 }        

View Code

这个代码还没有A啊，今天BZ挂了。贴上来玩耍

转载于:https://www.cnblogs.com/registerzxr/p/5130634.html