数论
文章平均质量分 82
_lifehappy_
这个作者很懒,什么都没留下…
展开
-
G. GCD Festival(莫比乌斯、欧拉函数)
G. GCD Festival∑i=1n∑j=1ngcd(ai,aj)gcd(i,j)∑d=1nd∑i=1nd∑j=1ndgcd(aid,ajd)[gcd(i,j)=1]∑d=1nd∑k=1ndμ(k)∑i=1nkd∑j=1nkdgcd(aikd,ajkd)T=kd∑T=1n∑i=1nT∑j=1nTgcd(aiT,ajT)∑d∣Tdμ(Td)∑T=1nϕ(T)∑i=1nT∑j=1nTgcd(aiT,ajT)\sum_{i = 1} ^{n} \sum_{j = 1} ^{n} \gcd(a原创 2021-10-07 20:19:43 · 561 阅读 · 0 评论 -
K. Easy Sigma(类欧几里得)
K. Easy Sigma∑i=1n(−1)⌊i×k⌋,(n≤109,k≤104)\sum_{i = 1} ^{n} (-1) ^{\lfloor i \times \sqrt k \rfloor}, (n \le 10 ^ 9, k \le 10 ^ 4)\\i=1∑n(−1)⌊i×k⌋,(n≤109,k≤104)考虑(−1)x=1−2×(xmod 2)=1−2(x−2×x2)=1−2x+4×⌊x2⌋(-1) ^{x} = 1 - 2 \times (x \mod 2) = 1 - 2(x原创 2021-10-07 20:18:40 · 317 阅读 · 0 评论 -
类欧几里得(模板题推导)
类欧几里得设三个函数f(a,b,c,n)=∑i=0na×i+bc,g(a,b,c,n)=∑i=0ni×a×i+bc,h(a,b,c,n)=∑i=0n(a×i+bc)2f(a, b, c, n) = \sum\limits_{i = 0} ^{n} \frac{a \times i + b}{c}, g(a, b, c, n) = \sum\limits_{i = 0} ^{n} i \times \frac{a \times i + b}{c}, h(a, b, c, n) = \sum\limits_{原创 2021-10-02 19:59:46 · 306 阅读 · 0 评论 -
#138. 类欧几里得算法
#138. 类欧几里得算法以下除法均为向下取整, 定义f(a,b,c,n,k1,k2)=∑x=0nxk1(a×x+bc)k2f(a, b, c, n, k_1, k_2) = \sum\limits_{x = 0} ^{n} x ^{k_1} \left(\frac{a \times x + b}{c}\right) ^ {k_2}f(a,b,c,n,k1,k2)=x=0∑nxk1(ca×x+b)k2。∑x=0nxk1(a×x+bc)k2,(1≤n,a,c≤109,0≤b≤109,0≤k1+原创 2021-10-02 18:00:41 · 281 阅读 · 0 评论 -
I. Rise of Shadows(类欧几里得)
I. Rise of Shadows一天有HHH个小时,MMM分钟,问,有多少个整数分钟,满足时针与分针的角度≤α\le \alpha≤α,α=2πAHM\alpha = \frac{2 \pi A}{HM}α=HM2πA。∑i=0H−1∑j=0M−1[∣2π(i×M+j)HM−2πjM ∣≤2πAHM]∑i=0H−1∑j=0M−1[∣i×M+j−H×j∣≤A]H×M−∑i=0H−1∑j=0M−1[i×M+j−H×j>A]−∑i=0H−1∑j=0M−1[i×M+j−H×j<−A]原创 2021-10-02 16:38:20 · 864 阅读 · 1 评论 -
C - Maximize GCD(简单数论)
C - Maximize GCD给定长度为n,(2≤3×105)n, (2 \le 3 \times 10 ^ 5)n,(2≤3×105)的数组a,(1≤ai≤3×105)a, (1 \le a_i \le 3 \times 10 ^ 5)a,(1≤ai≤3×105),一个数字K,(1≤K≤1018)K, (1 \le K \le 10 ^{18})K,(1≤K≤1018),我们可以对数组aaa进行最多kkk次操作,每次操作选定一个i,(1≤i≤n)i, (1 \le i \le n)i,(1≤i≤n原创 2021-09-23 19:59:59 · 321 阅读 · 0 评论 -
E. Mocha and Stars(莫比乌斯反演、简单dp)
E. Mocha and Stars∑a1=l1r1∑a2=l2r2⋯∑an=lnrn[a1+a2+⋯+an≤m][gcd(a1,a2,…,an)=1]\sum_{a_1 = l_1} ^{r_1} \sum_{a_2 = l_2} ^{r_2} \dots \sum_{a_n = l_n} ^{r_n} [a_1 + a_2 + \dots + a_n \le m][\gcd(a_1, a_2, \dots, a_n) = 1]\\a1=l1∑r1a2=l2∑r2⋯an=ln∑原创 2021-08-16 20:28:37 · 570 阅读 · 0 评论 -
Convolution(2021牛客暑期多校训练营4)
Convolution定义a⊕b=a×bgcd(a,b)2a \oplus b = \frac{a \times b}{\gcd(a, b) ^ 2}a⊕b=gcd(a,b)2a×b,bm=∑i=1n∑j=1nai×jc[i⊕j=m]b_m = \sum\limits_{i = 1} ^{n} \sum\limits_{j = 1} ^{n}a_i \times j ^ c [i \oplus j = m]bm=i=1∑nj=1∑nai×jc[i⊕j=m],我们要求b1⊕b2⊕⋯⊕bnb_1原创 2021-08-08 13:50:32 · 286 阅读 · 0 评论 -
P3700 [CQOI2017]小Q的表格(反演、分块)
P3700 [CQOI2017]小Q的表格给定一个大小为n×nn \times nn×n的表格,初始时i,ji, ji,j位置上填的是f(i,j)=i×jf(i, j) = i \times jf(i,j)=i×j,有mmm个操作,每次操作给定a,b,x,ka, b, x, ka,b,x,k,把格子a,ba, ba,b上的值改成xxx,求∑i=1k∑j=1kf(i,j)\sum\limits_{i = 1} ^{k} \sum\limits_{j = 1} ^{k} f(i, j)i=1∑kj=1∑k原创 2021-08-07 19:24:00 · 228 阅读 · 0 评论 -
HDU6956-Pass!(2021杭电多校一)(BSGS)
Pass!f(1)=0,f(2)=n−1,f(t)=(n−2)×f(t−1)+(t−1)×f(t−2)f(1) = 0, f(2) = n - 1, f(t) = (n - 2) \times f(t - 1) + (t - 1) \times f(t - 2)f(1)=0,f(2)=n−1,f(t)=(n−2)×f(t−1)+(t−1)×f(t−2),考虑对通项两边同时加一个x×f(t−1)x \times f(t - 1)x×f(t−1)。可以得到f(t)+x×f(t−1)=(n−1+x)×(f(t原创 2021-07-21 23:37:58 · 625 阅读 · 0 评论 -
J. Product of GCDs(莫比乌斯反演)(2021牛客暑期多校训练营2)
Product of GCDs∏d=1nd∑[gcd(s1d,s2d,…,skd)=1]∏d=1nd∑i=1ndμ(i)Cf[id]k\prod_{d = 1} ^{n} d ^{\sum[\gcd(\frac{s_1}{d}, \frac{s_2}{d}, \dots, \frac{s_k}{d}) = 1]}\\\prod_{d = 1} ^{n} d ^{\sum\limits_{i = 1} ^{\frac{n}{d}} \mu(i) C_{f[id]} ^{k}}\\d=1∏nd∑[g原创 2021-07-21 20:19:46 · 317 阅读 · 0 评论 -
P5175 数列(矩阵快速幂)
P5175 数列anb=(x×an−1+y×an−2)2x2×an−12+y2×an−22+2×x×y×an−1an−2x2×an−12+y2×an−22+2×x×y×an−2(x×an−2+y×an−3)x2×an−12+y2×an−22+2×x×y×(x×an−22+y×an−2×an−3)a_n ^ b = \left(x \times a_{n - 1} + y \times a_{n - 2}\right) ^ 2\\x ^ 2 \times a_{n - 1} ^ 2 + y ^ 2 \t原创 2021-07-13 20:04:05 · 250 阅读 · 0 评论 -
E. Surprise me!(莫比乌斯反演 + 虚树 DP)
E. Surprise me!∑i=1n∑j=1nϕ(ai×aj)d(i,j)设pai=i∑i=1n∑j=1nϕ(i×j)d(pi,pj)∑i=1n∑j=1nϕ(i)ϕ(j)ϕ(gcd(i,j))×gcd(i,j)×d(pi,pj)∑d=1ndϕ(d)∑i=1nd∑j=1ndϕ(id)ϕ(jd)×d(pid,pjd)[gcd(i,j)=1]∑d=1ndϕ(d)∑k=1ndμ(k)∑i=1nkd∑j=1nkdϕ(ikd)ϕ(jkd)×d(pikd,pjkd)T=kd∑T=1n(∑i=1nT∑j=1nT原创 2021-07-08 21:19:04 · 308 阅读 · 3 评论 -
数列递推(牛客练习赛83)(数学、分块)
数列递推给定f(0)f(0)f(0),定义fn=∑i=1nf(nmod i)f_n = \sum\limits_{i = 1} ^{n} f_{(n \mod i)}fn=i=1∑nf(nmodi),求f1,f2,f3,…,fn−1,fnf_1, f_2, f_3, \dots, f_{n - 1}, f_{n}f1,f2,f3,…,fn−1,fn。∑i=1nf(nmod i)∑i=1nf(n−nii)\sum_{i = 1} ^{n} f_{(n \mod i)}\\\sum_原创 2021-05-24 15:10:51 · 283 阅读 · 0 评论 -
A. Slackline Adventure(思维 + 莫比乌斯)(2018-2019 ACM-ICPC Brazil Subregional Programming Contest)
A. Slackline Adventure考虑枚举每个矩形的长跟宽,再统计这个矩形能在坐标轴上出现几次(同行相邻同列相邻的单独算),然后有如下式子:2×∑i=1n−1∑j=1m−1(n−i)(m−j)[gcd(i,j)=1][L2≤i2+j2≤R2]2 \times \sum_{i = 1} ^{n - 1} \sum_{j = 1} ^{m - 1}(n - i)(m - j)[\gcd(i, j) = 1][L ^ 2 \leq i ^ 2 + j ^ 2 \leq R ^ 2]\\2×i=原创 2021-05-01 23:45:02 · 349 阅读 · 0 评论 -
Mr. Panda and Kakin(拓展欧几里得 + O(1)快速乘)
Mr. Panda and Kakin给定n,cn, cn,c,要我们找到nnn是两个相邻质数的乘积,要我们找到xxx,满足x230+3≡c(modn)x ^{2 ^{30} + 3} \equiv c \pmod nx230+3≡c(modn),1010≤n≤1018,0<c<n10 ^{10} \leq n \leq 10 ^ {18}, 0 < c < n1010≤n≤1018,0<c<n,考虑得到230+32 ^{30} + 3230+3模ϕ(n)\phi(n原创 2021-04-28 19:20:48 · 336 阅读 · 0 评论 -
小 Q 与函数求和 1(牛客练习赛 81 E)
小 Q 与函数求和 1∑i=1n∑j=1nϕ(ijgcd(i,j)K)∑i=1n∑j=1ngcd(i,j)Kϕ(ij)∑i=1n∑j=1ngcd(i,j)Kϕ(i)ϕ(j)gcd(i,j)ϕ(gcd(i,j))∑i=1n∑j=1ngcd(i,j)K+1ϕ(i)ϕ(j)ϕ(gcd(i,j))∑d=1ndK+1inv(ϕ(d))∑i=1nd∑j=1ndϕ(id)ϕ(jd)[gcd(i,j)=1]∑d=1ndK+1inv(phi(d))∑k=1ndϕ(k)∑i=1nkdϕ(ikd)∑j=1nkd原创 2021-04-27 19:39:59 · 275 阅读 · 0 评论 -
D. Steps to One(概率DP,莫比乌斯反演)
D. Steps to One设f[i]f[i]f[i]为gcd\gcdgcd为iii,还需要多少个数,那么有f[i]=1+∑j=1mf[gcd(i,j)]mf[i] = 1 + \frac{\sum\limits_{j = 1} ^{m} f[\gcd(i, j)]}{m}f[i]=1+mj=1∑mf[gcd(i,j)],f[1]=0f[1] = 0f[1]=0,考虑化简∑j=1mf[gcd(i,j)]\sum\limits_{j = 1} ^{m} f[\gcd(i, j)]j=1∑mf原创 2021-04-14 17:28:10 · 306 阅读 · 0 评论 -
ZOJ The Sum of Unitary Totient(min_25 筛)
The Sum of Unitary Totient积性函数,满足质数点是多项式,直接 min_25 了,由于单次求解,所以使用递归的 min_25 会较快。#include <bits/stdc++.h>using namespace std;const int N = 1e5 + 10;int prime[N], a[N], id1[N], id2[N], cnt, n, T, m;long long sum1[N], sum2[N], g1[N], g2[N], s[N原创 2021-04-13 20:39:27 · 233 阅读 · 0 评论 -
P6271 [湖北省队互测2014]一个人的数论(莫比乌斯反演 + 伯努利数)
P6271 [湖北省队互测2014]一个人的数论∑i=1nim[gcd(i,n)=1]∑d∣nμ(d)dm∑i=1ndim由伯努利数可知∑i=0nim=1m+1∑i=0mCm+1iBi(n+1)m−i+1设fi=1m+1Bm−i+1Cm+1i,则有∑d∣nμ(d)dm(∑i=1m+1fi(nd)i+(nd)m)∑i=1m+1fini∑d∣nμ(d)dm−i+(nm∑d∣nμ(d))考虑后项∑d∣nμ(d)dm−i,迪利克雷卷积乘积,为积性函数F(n)=∑d∣nμ(d)dm−i,F(1)=1,F(pk)=原创 2021-04-10 11:35:19 · 235 阅读 · 0 评论 -
min_25 推导及例题总结
min_25 筛一个亚线性筛,复杂度大概是O(n34logn)O(\frac{n ^{\frac{3}{4}}}{ \log n})O(lognn43)。使用min_25min\_25min_25求前缀和,有两个基本特征:① 积性函数,② 满足质数点为多项式。算法思路给定n≤1011n \leq 10 ^ {11}n≤1011,设f(n)f(n)f(n)为一个积性函数,质数点为多项式,求∑i=1nf(i)\sum\limits_{i = 1} ^{n} f(i)i=1∑nf(i)。先原创 2021-03-23 19:54:30 · 494 阅读 · 0 评论 -
#3328. PYXFIB(单位根反演)
#3328. PYXFIB∑i=0⌊nk⌋Cni×k×Fi×k∑i=0nCni×Fi×[i≡0(modk)]i≡0(modk),单位根反演有1k∑j=0k−1wkij1k∑i=0nCni×Fi∑j=0k−1wkij\sum_{i = 0} ^{\lfloor \frac{n}{k} \rfloor} C_{n} ^{i \times k} \times F_{i \times k}\\\sum_{i = 0} ^{n} C_n ^{i} \times F_i \times [i \equiv 0 \p原创 2021-03-23 14:57:51 · 253 阅读 · 0 评论 -
P5591 小猪佩奇学数学(单位根反演)
P5591 小猪佩奇学数学∑i=0n(in)×pi×⌊ik⌋⌊ik⌋=i−i%kk1k∑i=0n(in)×pi×(i−i%k)1k∑i=0n(in)×pi×i−1k∑i=0n(in)×pi(imod k)\sum_{i = 0} ^{n} (_i ^ n) \times p ^ i \times \lfloor \frac{i}{k} \rfloor\\\lfloor \frac{i}{k} \rfloor = \frac{i - i \% k}{k}\\\frac{1}{k} \sum_{i =原创 2021-03-21 13:36:35 · 270 阅读 · 0 评论 -
单位根反演小记
单位根反演一个等式:[n∣a]=1n∑k=0n−1wnak[n \mid a] = \frac{1}{n} \sum\limits_{k = 0} ^{n - 1}w_n ^{ak}[n∣a]=n1k=0∑n−1wnak证明:wnaw_n ^ awna是nnn次单位根的aaa次方,所以这里是一个公比为wnaw_n ^ awna的等比数列,考虑公比不为111,等比数列求和1nwn0(wnan−1)wna−1\frac{1}{n} \frac{w_n ^ 0(w_n ^ {an} - 1)}{原创 2021-03-21 13:34:04 · 276 阅读 · 0 评论 -
1190 最小公倍数之和 V2
1190 最小公倍数之和 V2∑i=ablcm(i,b)∑i=abi×bgcd(i,b)b∑d∣b∑i=⌈ad⌉bdi[gcd(i,bd)=1]b∑d∣b∑k∣bdμ(k)k∑i=⌈⌈ad⌉k⌉abkib∑T∣n∑i=⌈aT⌉bTi∑k∣Tμ(k)kb∑T∣n(bT+⌈aT⌉)×(bT−⌈aT⌉+1)2∑k∣Tμ(k)k设f(n)=∑d∣nμ(d)d,f(1)=1,f(p)=1−p,f(pk)=1−p,且为积性函数\sum_{i = a} ^{b} lcm(i, b)\\\sum_{i = a} ^原创 2021-03-13 20:44:44 · 267 阅读 · 0 评论 -
HDU 6428 Problem C. Calculate(积性函数)
Problem C. Calculateϕ=ϕ∗ϵ=ϕ∗μ∗Iϕ(n)=∑d∣n(ϕ∗μ)(d)设g(n)=∑d∣n(ϕ∗μ)(d)∑i=1A∑j=1B∑k=1Cϕ(gcd(i,j2,k3))∑i=1A∑j=1B∑k=1C∑d∣i,d∣j2,d∣k3(ϕ∗μ)(d)∑d=1A(ϕ∗μ)(d)∑i=1A∑j=1B∑k=1C[d∣i,d∣j2,d∣k3]\phi = \phi * \epsilon = \phi * \mu * I\\\phi(n) = \sum_{d \mid n} (\phi*\mu)原创 2021-03-06 15:01:08 · 266 阅读 · 0 评论 -
#2693. jzptab
jzptab∑i=1n∑j=1mlcm(i,j)∑i=1n∑j=1mijgcd(i,j)∑d=1nd∑i=1nd∑j=1mdij[gcd(i,j)=1]∑d=1nd∑k=1ndk2μ(k)∑i=1nkdi∑j=1mkdjT=kd,f(n)=∑i=1ni∑T=1nf(nT)f(mT)(T∑k∣Tμ(k)k)设g(n)=n∑d∣nμ(d)d先令g(n)=g(n)ng(1)=1,g(p)=μ(1)+μ(p)p=1−p,g(pk,k≥2)=1−p同时是积性函数,可以O(n)求得,最后再乘上数组下标即可\su原创 2021-03-02 20:36:46 · 233 阅读 · 0 评论 -
2020 ICPC 济南 F. Gcd Product
Gcd ProductCm=∑i=1mAgcd(i,m)Bgcd(k+1−i,m)∑d1∣mAd1∑d2∣mBd2∑i=1m([gcd(id1,md1)=1][d1∣i])([gcd(m+1−id2,md2)=1][d2∣m+1−i])∑d1∣mAd1∑d2∣mBd2∑k1∣md1μ(k1)∑k2∣md2μ(k2)∑i=1m([d1∣i][k1∣id1])([d2∣m+1−i][k2∣m+1−id2])T1=d1×k1,T2=d2×k2∑T1∣m∑d1∣T1Ad1μ(T1d1)∑T2∣m∑d2∣T原创 2021-02-21 21:32:32 · 451 阅读 · 1 评论 -
NUMTRYE - Number Theory (Easy)
NUMTRYE - Number Theory (Easy)Hard 版本就是用 pollard_rho 分解质因子。f(n)=∏(pi2ei+1+1)f(n) = \prod(p_i ^{2e_i + 1} + 1)f(n)=∏(pi2ei+1+1),g(n)=∑i=1nngcd(n,i)g(n) = \sum\limits_{i = 1} ^{n} \frac{n}{\gcd(n, i)}g(n)=i=1∑ngcd(n,i)n,pip_ipi是nnn的质因子,eie_iei是pip_i原创 2021-02-20 15:27:44 · 204 阅读 · 0 评论 -
美味果冻(牛客练习赛53B)
美味果冻∑i=1n∑j=1ii×⌊ij⌋j∑i=1n∑j=inj×⌊ji⌋i\sum_{i = 1} ^{n} \sum_{j = 1} ^{i} i \times \lfloor \frac{i}{j} \rfloor ^ j\\\sum_{i = 1} ^{n} \sum_{j = i} ^{n} j \times \lfloor \frac{j}{i} \rfloor ^ i\\i=1∑nj=1∑ii×⌊ji⌋ji=1∑nj=i∑nj×⌊ij⌋i接下来只需要从小到大枚举iii,因原创 2021-01-22 22:14:49 · 284 阅读 · 2 评论 -
二次剩余模板
二次剩余模板x2≡n(modp)x ^ 2 \equiv n \pmod px2≡n(modp)/* Author : lifehappy*/#pragma GCC optimize(2)#pragma GCC optimize(3)#include <bits/stdc++.h>#define mp make_pair#define pb push_back#define endl '\n'#define mid (l + r >> 1)#define原创 2020-08-27 16:13:24 · 235 阅读 · 0 评论 -
拓展欧几里德算法模板
ll ex_gcd(ll a, ll b, ll & x, ll & y) { if(!b) { x = 1, y = 0; return a; } ll gcd = ex_gcd(b, a % b, x, y); ll temp = x; x = y; y = temp - a / b * y; return gcd;}int main() { // freopen("in.txt", "原创 2020-07-24 17:57:37 · 253 阅读 · 0 评论 -
pollard_rho 模板
ll gcd(ll a, ll b) { return b ? gcd(b, a % b) : a;}ll quick_mult(ll a, ll b, ll mod) { ll ans = 0; while(b) { if(b & 1) ans = (ans + a) % mod; a = (a + a) % mod; b >>= 1; } return ans;}ll quick_po原创 2020-07-24 16:05:48 · 328 阅读 · 0 评论 -
miller_rabin 模板
miller_rabin 模板ll quick_mult(ll a, ll b, ll mod) { ll ans = 0; while(b) { if(b & 1) ans = (ans + a) % mod; a = (a + a) % mod; b >>= 1; } return ans;}ll quick_pow(ll a, ll n, ll mod) { ll ans = 1;原创 2020-07-22 18:39:56 · 342 阅读 · 0 评论 -
欧拉筛模板
欧拉筛模板const int N = 1e5 + 10;vector<int> prime;bool st[N];void init() { st[0] = st[1] = 1; for(int i = 2; i < N; i++) { for(int j = 0; j < prime.size() && i * prime[j] < N; j++) { st[i * prime[j]] = 1原创 2020-07-22 18:08:17 · 298 阅读 · 0 评论 -
1575 Gcd and Lcm
1575 Gcd and Lcm∑i=1n∑j=1i∑k=1ilcm(gcd(i,j),gcd(i,k))设f(n)=∑i=1n∑j=1nlcm(gcd(i,n),gcd(j,n))f(p)=3p2−3p+1f(pk)=(2k+1)(p2k−p2k−1)+pk−1\sum_{i = 1} ^{n} \sum_{j = 1} ^{i} \sum_{k = 1} ^{i} lcm(gcd(i, j), gcd(i, k))\\设f(n) = \sum_{i = 1} ^{n} \sum_{j = 1} ^原创 2020-10-27 22:24:50 · 341 阅读 · 0 评论 -
Min_25筛有关求解次小质因子
#188. 【UR #13】Sanrd题意化简就是求次小质因子,这一步我们可以在Min_25筛的ans计算中得到,S(n, j)表示的是最小质因子大于等于primejprime_jprimej的加上质数的答案贡献,要满足次小质因子,一定有除去这个数之后只剩下质数了,显然当我们计算到S(n, j)时,我们可以计算primej−1prime_{j - 1}primej−1对答案的贡献,接着我们还需要加上质数的次方项对答案的贡献,最后我们就可以得到答案了。/* Author : lifehapp原创 2020-09-23 21:12:05 · 246 阅读 · 0 评论 -
「LibreOJ Round #11」Misaka Network 与求和(杜教筛 + Min_25)
#572. 「LibreOJ Round #11」Misaka Network 与求和推式子∑i=1n∑j=1nf(gcd(i,j))k∑d=1nf(d)k∑i=1nd∑j=1nd[gcd(i,j)=1]∑d=1nf(d)k∑K=1ndμ(k)(nKd)2t=Kd∑t=1n(nt)2∑d∣tf(d)kμ(td)我们记f(x)k=F(x)上面式子后半部分是一个迪利克雷卷积形式:F∗μ所以我们卷上一个I,有F∗μ∗I=F∗ϵ=F得到后半部分的前缀和S(n)=∑i=1nF(i)−∑i=2nS(ni)\sum原创 2020-09-23 21:13:29 · 429 阅读 · 0 评论 -
Easy Math(ACM-ICPC 2018 徐州赛区网络预赛)(递归 + 杜教筛)
Easy Math推式子∑i=1mμ(in)∑i=1mμ(i×nd×d),d是n的一个质因子i,d互质项有(−∑i=1mμ(i×nd)),由于减去了多余的非互质项,所以加上,−∑i=1mμ(i×nd)+∑i=1mdμ(i×d×nd)−∑i=1mμ(i×nd)+∑i=1mdμ(i×n)\sum_{i = 1} ^{m} \mu(in)\\\sum_{i = 1} ^{m} \mu(i \times \frac{n}{d} \times d),d是n的一个质因子\\i, d互质项有(-\sum_{i原创 2020-09-12 11:03:57 · 276 阅读 · 0 评论 -
B. Product(2019ICPC西安邀请赛)(杜教筛)
Product推式子原式=m∑i=1n∑j=1n∑k=1ngcd(i,j)[k∣gcd(i,j)]∑i=1n∑j=1n∑k=1ngcd(i,j)[k∣gcd(i,j)]∑k=1n∑i=1nk∑j=1nkgcd(ik,jk)∑k=1nk∑i=1nk∑j=1nkgcd(i,j)∑k=1nk∑d=1nkd∑t=1nkdμ(t)(nkdt)2∑k=1nk∑d=1nkϕ(d)(nkd)2在这个点直接输出答案tle了几个点,所以还得化简。t=kd∑t=1n(nt)2∑d∣tϕ(d)td后面一项写成迪利克雷卷积形式f原创 2020-09-11 23:06:32 · 388 阅读 · 0 评论