FSYo 的博客

传送门以下整点皆在第一象限讨论令 f(x)=∑i≥1∑j≥1[i2+j2=x2]f(x)=\sum_{i\ge 1}\sum_{j\ge 1}[i^2+j^2=x^2]f(x)=∑i≥1​∑j≥1​[i2+j2=x2]，要求的就是 ∑i=1Nf(i)\sum_{i=1}^Nf(i)∑i=1N​f(i)考虑把整点 (x,y)(x,y)(x,y) 用复平面的高斯整数 x+yix+yix+yi 表示，而 (x+yi)(x−yi)=x2+y2(x+yi)(x-yi)=x^2+y^2(x+yi)(x−yi

今天被打爆了。。。题解单独写题意：求∑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...

T≤100,n≤1e6,∑k≤1e6T\le 100,n\le 1e6,\sum k\le 1e6T≤100,n≤1e6,∑k≤1e6简要题解：首先可以暴力 dpdpdp，fi,jf_{i,j}fi,j​ 表示到了第 iii 维，乘积为 jjj 的所有方案数的和fi,j=∑l∣jfi−1,l(jl)bif_{i,j}=\sum_{l|j}f_{i-1,l}(\frac{j}{l})^{b...

今天做了一道 powerful numberpowerful\ numberpowerful number 筛的题，暴力可以用 min25min25min25，然后我板子大挂 3 个大样例全过…题意：f(x)=1(x=1)f(x)=1(x=1)f(x)=1(x=1)f(x)=2Ω(x)(x>1)f(x)=2^{\Omega(x)}(x>1)f(x)=2Ω(x)...

传送门要求求积性函数 f(x) 的前缀和，f(prim) 是一个关于prim的简单多项式，f(prim^k) 可以快速计算求法：Min25 筛分为两部分，第一部分处理素数的幂在的前缀和即 Min25 的核心思想就是考虑小于根号n的质数可以出去大于根号n的合数的贡献因此我们可以利用这个性质避免计算大于根号n的数的幂设为已经筛了 j 个倍数，<= i ...

传送门考虑每一个 Min(x) 时的贡献，发现就是 Min25 筛时没有被考虑的数的和 / Min(x)然后Min25 筛时顺便筛出答案即可#include<bits/stdc++.h>#define N 1000050using namespace std;typedef long long ll;ll l, r, prim[N]; int tot; boo...

传送门, 然后就是Min25的板子题，要特殊考虑 2 的影响#include<bits/stdc++.h>#define N 400050using namespace std;typedef long long ll;const int Mod = 1e9 + 7, inv2 = (Mod+1)/2;ll add(ll a, ll b){ return (a + ...

传送门首先 g 可以整除分块，然后杜教筛处理 g考虑如何用 Min25 筛出 f，一个 d 有 f(d) ^ k 的贡献，当且仅当有一个 >= 它的质数在它前面所有的方案数就是 >= 它的质数个数于是 Min25 先处理出质数个数，然后加上合数贡献时算一下 f(d)具体一点就是#include<bits/stdc++.h>#define ...

9.29 重学 min25fi,jf_{i,j}fi,j​ 表示已经筛去前 j 个素数的倍数后 ≤i\le i≤i 的数的 k 次幂之和fi,j=∑x=2ixk[Min(x)>pj∣x∈P]f_{i,j}=\sum_{x=2}^i x^k[Min(x)>p_j|x\in P]fi,j​=∑x=2i​xk[Min(x)>pj​∣x∈P]转移：fi,j=fi,j−1−pjk(f...