hdu 6076 huntian oy 杜教筛

打表观察得到,gcd(i,j)==1时,gcd(i^a−j^a,i^b−j^b)的值为i - j。所以,你发现这个题跟ab就没关系了...

变成去求∑∑(i-j)[gcd(i,j) == 1]了。有一个显然的结论,gcd(i,j) == gcd(i-j,i)。

设k为i-j,则变成

∑(i 1->n)∑k(1->i-1)[gcd(i,k) == 1]。

又因为i和i本身,一定不互质,所以变成

∑(i 1->n)∑k(1->i)[gcd(i,k) == 1]。然后有个结论,1…N中与N互质的数的和,ans=N*phi(N)/2。

然后就变成了求∑phi(i)*i,求前缀和想到了杜教筛。卷积一下函数g(x) = x,就行了。

 1 #include <cstdio>
 2 #include <map>
 3 #include <cmath>
 4 using namespace std;
 5 typedef long long ll;
 6 const int MAXN = 1000100,mo = 1e9 + 7,inv2 = (mo + 1) / 2,inv6 = (mo + 1) / 6;
 7
 8 int T,n,a,b,maxn,phi[MAXN],sum[MAXN],pri[MAXN];
 9 bool vis[MAXN];
10 map<int,int> f;
11 int solve(int n)
12 {
13     if (n <= maxn)
14         return sum[n];
15     if (f.count(n))
16         return f[n];
17     int ans = (ll)n * (n + 1) % mo * (n * 2 + 1) % mo * inv6 % mo;
18     for (int l = 2,r;l <= n;l = r + 1)
19     {
20         r = n / (n / l);
21         int tp = (ll)(l + r) * (r - l + 1) % mo * inv2 % mo;
22         ans -=  (ll)tp * solve(n / l) % mo;
23         if (ans < 0)
24             ans += mo;
25     }
26     return f[n] = ans;
27 }
28 void init()
29 {
30     maxn = 1000000;
31     phi[1] = 1;
32     int tot = 0;
33     for (int i = 2;i <= maxn;i++)
34     {
35         if (vis[i] == false)
36         {
37             pri[++tot] = i;
38             phi[i] = i - 1;
39         }
40         for (int j = 1;j <= tot && i * pri[j] <= maxn;j++)
41         {
42             vis[i * pri[j]] = true;
43             if (i % pri[j] != 0)
44                 phi[i * pri[j]] = phi[i] * (pri[j] - 1);
45             else
46             {
47                 phi[i * pri[j]] = phi[i] * pri[j];
48                 break;
49             }
50         }
51     }
52     for (int i = 1;i <= maxn;i++)
53         sum[i] = (sum[i - 1] + (ll)phi[i] * i % mo) % mo;
54 }
55
56 int main()
57 {
58     init();
59     for (scanf("%d",&T);T != 0;T--)
60     {
61         scanf("%d%d%d",&n,&a,&b);
62         printf("%d\n",(ll)(solve(n) - 1 + mo) % mo * inv2 % mo);
63     }
64     return 0;
65 }

原文地址:https://www.cnblogs.com/iat14/p/11407855.html

时间: 2024-09-29 06:59:20

hdu 6076 huntian oy 杜教筛的相关文章

CCPC 2019 网络赛 HDU huntian oy (杜教筛)

1005 huntian oy (HDU 6706) 题意: 令,有T次询问,求 f(n, a, b). 其中 T = 10^4,1 <= n,a,b <= 1e9,保证每次 a,b互质. 思路: 首先我们需要知道 公式: gcd(a^n - b^n, a^m - b^m) = a^(gcd(m,n)) - b^(gcd(m,n)) 由a,b互质,原式即为 f(n, a, b) = ∑∑ (i-j)*[(i,j)=1] = ∑ (i*∑ [(i, j)=1] ) - ∑∑ j*[(i, j)=

HDU6706 CCPC 2019网络赛 huntian oy 推式子+杜教筛

CCPC 2019 网络赛 HDU 6706 huntian oy 标签 奇奇怪怪的数论结论 杜教筛 前言 我的csdn和博客园是同步的,欢迎来访danzh-博客园~ 简明题意 给定n,a,b,求: \[\sum_{i=1}^n\sum_{j=1}^igcd(i^a-j^a,i^b-j^b)[gcd(i,j)=1]\%(10^9+7)\] 思路 首先有一个结论: \[gcd(i^a-j^a,i^b-j^b)=i^{gcd(a,b)}-j^{gcd(a,b)}\] 上面的结论对于i,j互质是成立的

杜教筛 学习总结

看了看唐老师的blog,照猫画虎的做了几道题目,感觉对杜教筛有些感觉了 但是稍微有一点难度的题目还是做不出来,放假的时候争取都A掉(挖坑ing) 这篇文章以后等我A掉那些题目之后再UPD上去就好啦 由于懒得去写怎么用编辑器写公式,所以公式就准备直接copy唐老师的啦 首先积性函数和完全积性函数什么的就不再多说了 列举常见的积性函数: 1.约数个数函数和约数个数和函数 2.欧拉函数phi 3.莫比乌斯函数mu 4.元函数e 其中e(n)=[n==1] 5.恒等函数I 其中I(n)=1 6.单位函数

【bzoj 4176】 Lucas的数论 莫比乌斯反演(杜教筛)

Description 去年的Lucas非常喜欢数论题,但是一年以后的Lucas却不那么喜欢了. 在整理以前的试题时,发现了这样一道题目“求Sigma(f(i)),其中1<=i<=N”,其中 表示i的约数个数.他现在长大了,题目也变难了. 求如下表达式的值: 一行一个整数ans,表示答案模1000000007的值. Sample Input 2 Sample Output 8 HINT 对于100%的数据n <= 10^9. 题解: 解锁新技能:杜教筛. 再复习一下: 若$F(n)=\s

【51nod-1239&amp;1244】欧拉函数之和&amp;莫比乌斯函数之和 杜教筛

题目链接: 1239:http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1239 1244:http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1244 杜教筛裸题,不过现在我也只会筛这俩前缀和... $$s(n)=\sum _{i=1}^{n}f(i)$$ 那么就有: $$\sum_{i=1}^{n}f(i)\lfloor \frac{n}{i} \

【BZOJ4916】神犇和蒟蒻 杜教筛

[BZOJ4916]神犇和蒟蒻 Description 很久很久以前,有一只神犇叫yzy; 很久很久之后,有一只蒟蒻叫lty; Input 请你读入一个整数N;1<=N<=1E9,A.B模1E9+7; Output 请你输出一个整数A=\sum_{i=1}^N{\mu (i^2)}; 请你输出一个整数B=\sum_{i=1}^N{\varphi (i^2)}; Sample Input 1 Sample Output 1 1 题解:哎?上面的那个东西好像一直是1?(废话),然后 设j=i/d,

【bzoj4176】Lucas的数论 莫比乌斯反演+杜教筛

题目描述 去年的Lucas非常喜欢数论题,但是一年以后的Lucas却不那么喜欢了. 在整理以前的试题时,发现了这样一道题目“求Sigma(f(i)),其中1<=i<=N”,其中 表示i的约数个数.他现在长大了,题目也变难了. 求如下表达式的值: 其中f(ij)表示ij的约数个数. 他发现答案有点大,只需要输出模1000000007的值. 输入 第一行一个整数n. 输出 一行一个整数ans,表示答案模1000000007的值. 样例输入 2 样例输出 8 题解 莫比乌斯反演+杜教筛 首先有个神奇

杜教筛进阶+洲阁筛讲解+SPOJ divcnt3

Part 1:杜教筛进阶在了解了杜教筛基本应用,如$\sum_{i=1}^n\varphi(i)$的求法后,我们看一些杜教筛较难的应用.求$\sum_{i=1}^n\varphi(i)*i$考虑把它与$id$函数狄利克雷卷积后的前缀和.$$\sum_{i=1}^n\sum_{d|i}\varphi(d)*d*\frac id=\sum_{i=1}^ni^2$$枚举$T=\frac id$,原式化为$$\sum_{T=1}^nT\sum_{d=1}^{\lfloor\frac nT\rfloor}

【51Nod 1237】最大公约数之和 V3 莫比乌斯反演+杜教筛

题意 求$\sum_{i=1}^{n}\sum_{j=1}^{n}(i,j)$ 枚举约数 $$ \begin{align} ans &=\sum_{d=1}^{n}\sum_{i=1}^{n}\sum_{j=1}^{n}[(i,j)=d] \ &=\sum_{d=1}^{n}\sum_{i=1}^{\lfloor \frac{n}{d} \rfloor}\sum_{j=1}^{\lfloor \frac{n}{d} \rfloor}[(i,j)=1] \ \end{align} $$ 利用