FSYo 的博客

今天被打爆了。。。题解单独写 题意：求 ∑k=1n∑i=1k∑j=1kgcd(i,j,k)\sum_{k=1}^n\sum_{i=1}^k\sum_{j=1}^kgcd(i,j,k)k=1∑n​i=1∑k​j=1∑k​gcd(i,j,k) 简单反演 Ans=∑k∑d∣kd∑i∑j[gcd(i,j,k)=d]=∑k∑d∣kd∑ik/d∑jk/d[gcd(i,j,k/d)=1]=∑k∑d∣kd...

传送门 考虑什么样的 xy\frac{x}{y}yx​ 可以成为纯循环小数 设其循环节为 LLL，那么有 xy∗kL−xy\frac{x}{y}*k^L-\frac{x}{y}yx​∗kL−yx​ 为整数 每一对贡献在 gcd(x,y)=1gcd(x,y)=1gcd(x,y)=1 的时候统计，于是上面这个条件可以转换为 ∃L,s.t,kL−1≡0(mod y)\exists L,s.t,...

题意：给定 n,m,dn,m,dn,m,d，求 ∑i=1n∑j=1m(gcd(i,j)B)\sum_{i=1}^n\sum_{j=1}^m\binom{gcd(i,j)}{B}i=1∑n​j=1∑m​(Bgcd(i,j)​) 直接反演 ∑d=1n(dB)∑l=1ndμ(l)⌊ndl⌋⌊mdl⌋\sum_{d=1}^n\binom{d}{B}\sum_{l=1}^{\frac{n}{d}}\mu (...

传送门 先从一道简单题入手 (LCM sum) 求  令k=i/d gcd里面同除d 我们发现枚举到n/d 与 d是等价的 这道题就可以转化为 杜教筛筛第一坨, 然后整除分块 #include<bits/stdc++.h> #define N 1000050 #define LL long long #define Mod 1000...

传送门 反演一波 枚举k,l 发现n/t 总共有根号n个取值, 每算一次前缀和是根号n的 可以证明时间复杂度是n^3/4, mu用杜教筛 #include<bits/stdc++.h> #define N 2000050 #define Mod 1000000007 #define LL long long #define M 50050 using nam...

传送门 然后就是整除分块+杜教筛 #include<bits/stdc++.h> #define N 5000050 #define Mod 1000000007 #define LL long long using namespace std; int n,k,L,R; int prim[N], isp[N], tot, mu[N]; map<int,int&...

传送门 记住这个性质 那么原式就是   考虑求  然后整除分块一波, 求 的前缀和就可以了, 杜教筛解决 #include<bits/stdc++.h> #define N 5000050 #define LL long long #define Mod 1000000007 using namespace std; int prim[N], tot, is...

传送门       用杜教筛筛就可以了(g 设为Id^2)  #include<bits/stdc++.h> #define N 10000050 #define LL long long using namespace std; LL P,n; int prim[N],isp[N],tot; LL phi[N]; map<LL,LL> Phi; LL in...

传送门 杜教筛解决这类问题:  求            其中   我们寻找一个函数 g, 看能不能求出   看看我们要什么 (g为积性函数 -> g(1) = 1)  于是要求的变成 我们选出g 使第一个很好求, 然后第二个可以整除分块递归下去 本题的解决方案   当 f 为 mu 时, 考虑 令g = I, 那么   当 f 为 phi ...