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/dgcd里面同除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 50050using nam...

传送门然后就是整除分块+杜教筛#include<bits/stdc++.h>#define N 5000050#define Mod 1000000007#define LL long longusing 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 1000000007using namespace std;int prim[N], tot, is...

传送门   用杜教筛筛就可以了(g 设为Id^2) #include<bits/stdc++.h>#define N 10000050#define LL long longusing 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 ...