Crash的数字表格 BZOJ 2154 / jzptab BZOJ 2693

Posted 2020-09-13 草根柴鸡

jzptab

【问题描述】

求：

 

多组询问

【输入格式】

一个正整数T表示数据组数

接下来T行 每行两个正整数 表示N、M

【输出格式】

T行 每行一个整数 表示第i组数据的结果

【样例输入】

1
4 5

【样例输出】

122

【数据范围】

T <= 10000
N, M<=10000000

---

题解：

即后面那个部分为 H[T]，H[T]是积性函数，求详细证明的话将T和d展开为质因数次幂相乘的形式，考虑线性筛中枚举的质数与被筛数的性质即可

 1 #include<cmath>
 2 #include<cstdio>
 3 #include<cstdlib>
 4 #include<cstring>
 5 #include<iostream>
 6 #include<algorithm>
 7 using namespace std;
 8 typedef long long ll;
 9 const int maxn = 1e7 + 1;
10 const int mod = 1e8 + 9;
11 int cnt;
12 int h[maxn];
13 int pri[maxn];
14 int sum[maxn];
15 bool vis[maxn];
16 inline void Scan(int &x)
17 {
18     char c;
19     bool o = false;
20     while(!isdigit(c = getchar())) o = (c != \'-\') ? o : true;
21     x = c - \'0\';
22     while(isdigit(c = getchar())) x = x * 10 + c - \'0\';
23     if(o) x = -x;
24 }
25 inline void Sieve()
26 {
27     h[1] = 1;
28     for(int i = 2; i <= maxn; ++i)
29     {
30         if(!vis[i]) pri[++cnt] = i, h[i] = ((-(long long) i * i % mod) + mod + i) % mod;
31         for(int j = 1; j <= cnt; ++j)
32         {
33             int s = pri[j];
34             long long k = (long long) i * s;
35             if(k > maxn) break;
36             vis[k] = true;
37             if(!(i % s))
38             {
39                 h[k] = (long long) s * h[i] % mod;
40                 break;
41             }
42             else h[k] = (long long) h[s] * h[i] % mod;
43         }
44     }
45     for(int i = 1; i <= maxn; ++i) sum[i] = (sum[i - 1] + h[i]) % mod;
46 }
47 inline int Sum(int n, int m)
48 {
49     return ((long long) n * (n + 1) >> 1) % mod * (((long long) m * (m + 1) >> 1) % mod) % mod;
50 }
51 inline int Mobius(int n, int m)
52 {
53     int res = 0, last = 0;
54     if(n > m) swap(n, m);
55     for(int i = 1; i <= n; i = last + 1)
56     {
57         last = min(n / (n / i), m / (m / i));
58         res = (res + (long long) Sum(n / i, m / i) * ((sum[last] - sum[i - 1] + mod) % mod) % mod) % mod;
59     }
60     return res;
61 }
62 int main()
63 {
64     Sieve();
65     int n;
66     Scan(n);
67     int a, b;
68     while(n--)
69     {
70         Scan(a), Scan(b);
71         printf("%d\\n", Mobius(a, b));
72     }
73 }