莫比乌斯反演
v5zsq
新的一天,新的不会
展开
-
HDU 1695 GCD(莫比乌斯反演)
Description 求1<=x<=b,1<=y<=d且gcd(x,y)=k的对数,(x,y)和(y,x)视为一组 Input 第一行一整数T表示用例组数,每组用例五个整数a,b,c,d,k (T<=3000,a=c=1,b,d<=100000) Output 对于每组用例,输出1<=x<=b,1<=y<=d且gcd(x,y)=k的对数 Sample Input 2 1 3 1原创 2016-08-03 13:56:42 · 534 阅读 · 0 评论 -
CodeForces 235 E.Number Challenge(莫比乌斯反演+数论)
Description求∑i=1a∑j=1b∑k=1cd(ijk),a,b,c≤2000\sum\limits_{i=1}^a\sum\limits_{j=1}^b\sum\limits_{k=1}^cd(ijk),a,b,c\le 2000Input三个整数a,b,c(1≤a,b,c≤2000)a,b,c(1\le a,b,c\le 2000)Output输出结果模1073741824107374原创 2018-01-02 13:27:35 · 712 阅读 · 0 评论 -
CodeForces 803 F.Coprime Subsequences(莫比乌斯反演+组合数学)
Description给出一长度为nn的序列aia_i,问该序列的子序列中有多少子序列的gcdgcd是11Input第一行输入一整数nn表示序列长度,之后输入nn个整数aia_i表示该序列(1≤n,ai≤105)(1\le n,a_i\le 10^5)Output输出aa序列的gcdgcd为11的子序列,结果模109+710^9+7Sample Input3 1 2原创 2018-01-14 19:31:17 · 379 阅读 · 0 评论 -
CodeForces 839 D.Winter is here(莫比乌斯反演+组合数学)
Description给出nn个正整数a1,...,ana_1,...,a_n,从中选出kk个数,若其gcdgcd大于11,则贡献为k⋅gcdk\cdot gcd,求贡献和Input第一行一整数nn,之后输入nn个正整数ai(1≤n≤2⋅105,1≤ai≤106)a_i(1\le n\le 2\cdot 10^5,1\le a_i\le 10^6)Output输出贡献和,结果模原创 2018-01-14 20:29:04 · 543 阅读 · 0 评论 -
HDU 6159 Graph Of Zhuper(dp+CDQ分治+FFT+生成函数+莫比乌斯反演)
DescriptionGG为一个nn个点带标号无向图,无重边无自环,假设其有kk个连通分支,第ii个连通分支点数为sizeisize_i,定义GG的权值Zhu(G)=gcd(size1,size2,...,sizek)Zhu(G)=gcd(size_1,size_2,...,size_k),求∑Zhu(G)\sum Zhu(G)Input第一行一整数TT表示用例组数,每组用例输入两个整原创 2018-01-22 20:03:35 · 715 阅读 · 0 评论 -
CodeForces 439 E.Devu and Birthday Celebration(莫比乌斯反演+组合数学)
Description问满足gcd(a1,a2,...,af)=1,∑i=1fai=ngcd(a1,a2,...,af)=1,∑i=1fai=ngcd(a_1,a_2,...,a_f)=1,\sum\limits_{i=1}^fa_i=n的序列aaa的个数Input第一行一整数TTT表示用例组数,每组用例输入两个整数n,fn,fn,f表示一组查询(1≤T≤105,1≤f≤n≤105)(...原创 2018-01-08 21:57:25 · 444 阅读 · 0 评论 -
HDU 5942 Just a Math Problem(莫比乌斯反演+容斥原理)
Description已知f(k)f(k)f(k)为kkk的素因子个数,g(k)=2f(k)g(k)=2f(k)g(k)=2^{f(k)},给一正整数nnn,求∑i=1ng(i)∑i=1ng(i)\sum\limits_{i=1}^ng(i)Input第一行一整数TTT表示用例组数,每组用例输入一整数n(1≤T≤50,1≤n≤1012)n(1≤T≤50,1≤n≤1012)n(1\le ...原创 2018-04-08 11:26:26 · 316 阅读 · 0 评论 -
HDU 6340 Problem I. Delightful Formulas(莫比乌斯反演+数论+FFT)
Description给出n,Kn,Kn,K,令si=∑j≤ijKsi=∑j≤ijKs_i=\sum\limits_{j\le i}j^K,求∑1≤i≤nsi[gcd(i,n)=1]∑1≤i≤nsi[gcd(i,n)=1]\sum\limits_{1\le i\le n}s_i[gcd(i,n)=1]Input第一行一整数TTT表示用例组数,每组用例首先输入一整数KKK,之后输入一整数...原创 2018-09-06 08:16:54 · 368 阅读 · 0 评论 -
HDU 6363 bookshelf(数论+莫比乌斯反演)
Description将nnn本书等概率随机放在kkk层书架上,如果第iii层书架上有cnticnticnt_i本书,那么该层书架的稳定值为stablei=f[cnti]stablei=f[cnti]stable_i=f[cnt_i],其中f[0]=0,f[1]=1,f[i]=f[i−1]+f[i−2]f[0]=0,f[1]=1,f[i]=f[i−1]+f[i−2]f[0]=0,f[1]=1,...原创 2018-09-09 21:41:15 · 263 阅读 · 0 评论 -
HDU 6390 GuGuFishtion(数论+莫比乌斯反演)
Description求∑a=1m∑b=1nφ(ab)φ(a)φ(b)∑a=1m∑b=1nφ(ab)φ(a)φ(b)\sum\limits_{a=1}^m\sum\limits_{b=1}^n\frac{\varphi(ab)}{\varphi(a)\varphi(b)}Input第一行一整数TTT表示用例组数,每组用例输入三个整数n,m,pn,m,pn,m,p,其中ppp是素数(...原创 2018-09-09 21:45:18 · 316 阅读 · 0 评论 -
SPOJ 20173 DIVCNT2 - Counting Divisors (square)(数论+莫比乌斯反演)
Description给出一整数nnn,求S(n)=∑i=1nσ0(i2)S(n)=∑i=1nσ0(i2)S(n)=\sum\limits_{i=1}^n\sigma_0(i^2),其中σ0(x)σ0(x)\sigma_0(x)为xxx的因子数Input第一行一整数TTT表示用例组数,每组用例输入一整数n(1≤T≤104,1≤n≤1012)n(1≤T≤104,1≤n≤1012)n(1\...原创 2018-09-11 08:57:22 · 344 阅读 · 0 评论 -
HDU 6428 Problem C. Calculate(数论+莫比乌斯反演)
Description求∑i=1A∑j=1B∑k=1Cφ(gcd(i,j2,k3))(mod 230)\sum\limits_{i=1}^{A}\sum\limits_{j=1}^B\sum\limits_{k=1}^C\varphi(gcd(i,j^2,k^3))(mod\ 2^{30})i=1∑Aj=1∑Bk=1∑Cφ(gcd(i,j2,k3))(mod 230)...原创 2018-09-19 16:41:07 · 612 阅读 · 0 评论 -
Newcoder 148 E.Rikka with Equation(数论+莫比乌斯反演)
Description对于一个长度为nnn的正整数序列AAA和一个正整数mmm,定义f(A,m)f(A,m)f(A,m)为满足同余方程∑i=1nAixi≡0(mod m)\sum\limits_{i=1}^n A_ix_i\equiv 0(mod\ m)i=1∑nAixi≡0(mod m)的向量xi∈[0,m]x_i\in [0,m]xi∈[0,m]的个数。现...原创 2018-09-25 20:18:52 · 288 阅读 · 0 评论 -
BZOJ 3561 DZY Loves Math VI(莫比乌斯反演)
Description给定正整数n,mn,m,求∑i=1n∑j=1mlcm(i,j)gcd(i,j)\sum\limits_{i=1}^n\sum\limits_{j=1}^m lcm(i,j)^{gcd(i,j)}Input一行两个整数n,mn,m,1≤n,m≤5e51\le n,m\le 5e5Output一个整数,为答案模109+710^9+7后的值Sample Input5 4Sample原创 2017-12-22 15:55:26 · 376 阅读 · 0 评论 -
HDU 6088 Rikka with Rock-paper-scissors(莫比乌斯反演+组合数学+FFT)
Description两个人玩儿n局石头剪刀布,每局每个人出石头剪刀布的概率相同,玩nn局后假设第一个人赢aa局,第二个人赢bb局,那么得分就是gcd(a,b)gcd(a,b),设得分的期望为ss,求s×22ns\times 2^{2n},结果模modmodInput第一行一整数TT表示用例组数,每组用例输入两个整数nn和modmod表示游戏局数和模数(1≤T≤20,1≤n≤105,108≤mod≤原创 2017-09-05 22:32:50 · 873 阅读 · 0 评论 -
BZOJ 3529 数表(莫比乌斯反演+BIT)
Description 有一张n×m的数表,其第i行第j列(1<=i<=n,1<=j<=m)的数值为能同时整除i和j的所有自然数之和,给定a,计算数表中不大于a的数之和 Input 输入包含多组数据,输入的第一行一个整数Q表示测试点内的数据组数,接下来Q行,每行三个整数n,m,a(|a| < =10^9)描述一组数据 (1<=n,m<=10^5,1<=Q<=2×10^4) Output原创 2016-08-03 10:41:04 · 469 阅读 · 0 评论 -
BZOJ 2440 完全平方数(莫比乌斯反演+容斥原理+二分)
Description 小X自幼就很喜欢数。但奇怪的是,他十分讨厌完全平方数。他觉得这些数看起来很令人难受。由此,他也讨厌所有是完全平方数的正整数倍的数。然而这丝毫不影响他对其他数的热爱。这天是小X的生日,小W想送一个数给他作为生日礼物。当然他不能送一个小X讨厌的数。他列出了所有小X不讨厌的数,然后选取了第 K个数送给了小X。小X很开心地收下了。 然而现在小W却记不起送给小X的是哪个数了。你能帮他原创 2016-08-03 11:22:16 · 702 阅读 · 0 评论 -
SPOJ 4491 PGCD - Primes in GCD Table(莫比乌斯反演)
Description 给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对 Input 第一行一个整数T 表述数据组数接下来T行,每行两个正整数,表示N, M (T<=10000,N,M<=10000000) Output T行,每行一个整数表示第i组数据的结果 Sample Input 2 10 10 100 100原创 2016-08-03 12:36:58 · 504 阅读 · 0 评论 -
ZOJ 3435 Ideal Puzzle Bobble(莫比乌斯反演)
Description 给出3个数a,b,c, 定义一个立方体,这个立方体有a*b*c个点(一个端点在(1,1,1)而且处于第一象限),每个点的坐标都是整数(x,y,z),求经过坐标(1,1,1)和另外任意一个点(x1,y1,z1)的不同的直线有多少条 Input 多组用例(不超过200组),每组用例占一行包括三个整数a,b,c (2<=a,b,c<=1000000) Output 对于原创 2016-08-03 12:45:29 · 531 阅读 · 0 评论 -
SPOJ 7001 VLATTICE - Visible Lattice Points(莫比乌斯反演)
Description 求经过坐标(0,0,0)和另外任意一个点(x1,y1,z1)的不同的直线有多少条(0<=x1,y1,z1<=n) Input 第一行为一整数T表示用例组数,每组用例占一行为一整数n(T<=50,1<=n<=1000000) Output 对于每组用例,输出经过坐标(0,0,0)和另外任意一个点(x1,y1,z1)的不同的直线有多少条 Sample Inpu原创 2016-08-03 12:49:40 · 429 阅读 · 0 评论 -
BZOJ 2693 jzptab(莫比乌斯反演)
Description 求 Input 一个正整数T表示数据组数,接下来T行每行两个正整数表示N和M (T<=10000,N,M<=10000000) Output T行,每行一个整数表示第i组数据的结果 Sample Input 1 4 5 Sample Output 122 Solution Code#include<cstdio>#include<iostrea原创 2016-08-03 09:48:51 · 663 阅读 · 0 评论 -
BZOJ 2154 Crash的数字表格(莫比乌斯反演)
Description 今天的数学课上,Crash小朋友学习了最小公倍数(Least Common Multiple)。对于两个正整数a和b,LCM(a, b)表示能同时被a和b整除的最小正整数。例如,LCM(6, 8) = 24。回到家后,Crash还在想着课上学的东西,为了研究最小公倍数,他画了一张N*M的表格。每个格子里写了一个数字,其中第i行第j列的那个格子里写着数为LCM(i, j)。一原创 2016-08-03 13:15:06 · 651 阅读 · 0 评论 -
BZOJ 2301 Problem b(莫比乌斯反演)
Description 对于给出的n个询问,每次求有多少个数对(x,y),满足a≤x≤b,c≤y≤d且gcd(x,y)=k,gcd(x,y)函数为x和y的最大公约数 Input 第一行一个整数n,接下来n行每行五个整数,分别表示a,b,c,d,k Output 共n行,每行一个整数表示满足要求的数对(x,y)的个数 Sample Input 2 2 5 1 5 1 1 5 1 5原创 2016-08-03 13:23:09 · 422 阅读 · 0 评论 -
BZOJ 2820 YY的GCD(莫比乌斯反演)
Description 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对,kAc这种傻×必然不会了,于是向你来请教……多组输入 Input 第一行一个整数T 表述数据组数接下来T行,每行两个正整数,表示N, M (T<=10000,N,M<=10000000) Output T行,每行一个整数表示第i组原创 2016-08-03 11:06:01 · 487 阅读 · 0 评论 -
HDU 6053 TrickGCD(莫比乌斯反演)
Description 一个长度为n的序列b[1]~b[n]序满足gcd(b[l],…,b[r])>=2,1<=l<=r<=n,1<=b[i]<=a[i],给出a[1]~a[n],问满足条件的b序列个数 Input 第一行一整数T表示用例组数,每组用例首先输入一整数n表示序列长度,之后输入n个整数a[i] (1<=T<=10,1<=n,a[i]<=1e5) Output 对于每组用例,输出原创 2017-07-30 14:10:56 · 601 阅读 · 0 评论 -
HDU 6120 All Kill(数论+FFT+莫比乌斯反演)
Description 给一个非负序列x1,2,...,nx_{1,2,...,n},值小于32677,令yi,j=xi∗xj mod 32677y_{i,j}=x_{i}*x_{j} \ mod\ 32677,问有多少六元组(a,b,c,d,e,f)(a,b,c,d,e,f)满足gcd(ya,b,yc,d)=gcd(yc,d,ye,f)=gcd(ye,f,ya,b)=1gcd(y_{a,b},y原创 2017-08-16 16:06:35 · 598 阅读 · 0 评论 -
HDU 6134 Battlestation Operational(积性函数+莫比乌斯反演)
Description求f(n)=∑i=1n∑j=1i⌈ij⌉[(i,j)=1]f(n)=\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{i}\lceil \frac{i}{j}\rceil [(i,j)=1]Input多组用例,每组用例输入一整数nn,以文件尾结束(1≤n≤106)(1\leq n\leq10^{6})Output对于每组用例,输出f(n)f(n),原创 2017-08-17 18:41:38 · 1886 阅读 · 3 评论 -
HDU 6439 Congruence equation(莫比乌斯反演)
Description给出一个长度为kkk的序列A1,...,AkA_1,...,A_kA1,...,Ak,定义f(m)f(m)f(m)为满足以下条件的CCC序列的数量:1.若Ai=−1A_i=-1Ai=−1则CiC_iCi可以取[0,m)[0,m)[0,m)中任意整数,否则Ci=Ai%mC_i=A_i\%mCi=Ai%m2.同余方程∑i=1kCixi≡1(mod m)...原创 2018-09-27 08:18:46 · 618 阅读 · 0 评论