現在我們只要求出g(n)就能解決此題了，ans(n)列舉因子nlogn篩一下就好了
公式秒推，分塊求g(n):o(n*sqrtn*logn)秒t，嘗試分塊打表200*5000求g(n)：還是t了
嘗試oeis才發現g(n)有遞推公式g(n)=g(n-1)+d(n-1)+1，d(n)是n的因子個數，腦洞不夠啊
如果是向下取整求和也有公式=sum(d(i))
g++ 405ms
#include#define ll long long
using namespace std;
const int mod=1e9+7;
ll mu[1000005];bool isp[1000005];int pri[1000005];
ll ans[1000005];ll g[1000005];ll d[1000005];
int p=0;
void init(ll n)
else mu[i*pri[j]]=-mu[i];
} }}void init2(int n)
int main()
莫比烏斯反演 二 莫比烏斯反演定理
首先設兩個任意函式f x 和f x 定義運算 f x sum f d 這時就可以用f x 表示f x f 1 f 1 f 2 f 1 f 2 f 3 f 3 f 1 f 4 f 4 f 2 f 1 f 5 f 5 f 1 f 6 f 6 f 3 f 2 f 1 這時可以試著用f x 表示f x f ...
莫比烏斯反演
首先 莫比烏斯函式有個性質 d n d 1 n 1 0 n 1 證明 n 1時，不做多餘說明。n 1 根據唯一分解定理，可以分解n ki 1pai i 對於那些含平方因子也就是存在ai 不為1的數，它的函式值為0，對答案沒有任何貢獻。所以我們來看看那些是互異素數乘積的數，每乙個成為它約數的數是什麼樣...
莫比烏斯反演
定理 f n 和f n 是定義在非負整數集合上的兩個函式，並且滿足條件f n d nf d 那麼我們得到結論f n d n d f n d 在上面的公式中有乙個函式 d 它的定義如下 1 若d 1，那麼 d 1 2 若d p1 p2 p k 均為互異素數，那麼 d 1 k 3 其它情況下 d 0 對...