数论分块/整除分块

P1403 [AHOI2005] 约数研究 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

思路:板子.

第一次写，太死板了，写了个n*sqrt(n)的，TLE了
实际上题目问的问题可以转换为在[1,n]中，每个约数出现了多少次
每个约数出现的次数为n/x
这样直接套板子就是o(sqrt(n))的.
 1  2  3  4  5  6,7  8,9,10  [11,21]---约数
21 10  7  5  4   3     2        1   ---约数个数
不能太死板了
P1403 [AHOI2005] 约数研究
https://www.luogu.com.cn/problem/P1403
void solve(){               o(sqrt(n))
    int n; cin>>n;
    int ans=0;
    int r=0;
    for(int l=1;l<=n;l=r+1){
        r=n/(n/l);
        ans+=(r-l+1)*(n/l);
    }
    cout<<ans;
}

P2424 约数和 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

思路:用1到Y的贡献减去1到X-1的贡献，即是X到Y的贡献.

要清楚 L是什么 和 X/L是什么
ps:本题题目为"约数和"
P2424 "约数和"
https://www.luogu.com.cn/problem/P2424
void solve(){
    int x,y; cin>>x>>y;
    int r=0,ans=0;
    for(int l=1;l<=x-1;l=r+1){          [1,x-1]的贡献
        r=(x-1)/((x-1)/l);
        ans-=(r-l+1)*(l+r)/2*((x-1)/l);   等差数列*个数
    }
    r=0;
    for(int l=1;l<=y;l=r+1){          [1,y]的贡献
        r=y/(y/l);
        ans+=(r-l+1)*(l+r)/2*(y/l);   等差数列*个数
    }
    cout<<ans;
}

P2261 [CQOI2007] 余数求和 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

思路:

G(n,k)=∑(k%i) (i=[1,n])
化简k%i可得:ans=n*k-∑(k/i)*i   (i=[1,n])
那么∑(k/i)*i的贡献可以用数论分块求得

G(n,k)=∑(k%i) (i=[1,n])
化简k%i可得:ans=n*k-∑(k/i)*i   (i=[1,n])
那么∑(k/i)*i的贡献可以用数论分块求得
P2261 [CQOI2007] 余数求和
https://www.luogu.com.cn/problem/P2261
void solve(){
    int n,k; cin>>n>>k;
    int r=0,ans=n*k;
    for(int l=1;l<=n;l=r+1){
        if(l>k) break;   细节1:不能除以0
        r=min(n,k/(k/l));  细节2:r不能超过n,因为题目求的是l=[1,n]
        ans-=(r-l+1)*(l+r)/2*(k/l);  等差数列*(k/l)
    }
    cout<<ans;
}

K-取模_2022年中国高校计算机大赛-团队程序设计天梯赛（GPLT）上海理工大学校内选拔赛（重现赛）@IR101 (nowcoder.com)

思路:这题题意跟上一题一样，但是数据范围更大，到了1e12.就算一直取模也会爆long long，得用__int128才行.

用__int128才能ac.因为算贡献的时候已经超过long long类型的范围了.
∑(n%i) (i=[1,n])
化简得:ans=n*n-∑(n/i)*i   (i=[1,n])
取模
https://ac.nowcoder.com/acm/contest/74679/K
void solve(){
    const int mod=998244353;
    int m; cin>>m;
    __int128 n=m;
    __int128 r=0,res=n*n;
    for(int l=1;l<=n;l=r+1){
        r=n/(n/l);
        res-=(r-l+1)*(l+r)/2*(n/l)%mod,res+=mod,res%=mod;
    }
    int ans=res%mod;
    cout<<ans;
}