染色圖（數論分塊模板）

數論分塊模板（2sqrt(n）)適用於|n / i|向下取餘
因爲n/i一共有2sqrt(n)個值所以時間複雜度爲根號n
適用範圍：只要出現的式子有n/i則可以進行分塊

	//nl代表總的左邊界，nr是總的右邊界
	for(long long l=nl,r;l<=nr;l=r+1){
            r=min(n/(n/l),nr);
            long long t=n/l;
            //l代表當前塊的左邊界，r代表右邊界，t代表爲分塊的第幾塊 
        }

前根n樸素，後根n分塊
附數據便於理解 1-1000的分塊情況
1 1
2 2
3 3
4 4
5 5
6 6
7 7
8 8
9 9
10 10
11 11
12 12
13 13
14 14
15 15
16 16
17 17
18 18
19 19
20 20
21 21
22 22
23 23
24 24
25 25
26 26
27 27
28 28
29 29
30 30
31 31
32 32
33 33
34 34
35 35
36 37
38 38
39 40
41 41
42 43
44 45
46 47
48 50
51 52
53 55
56 58
59 62
63 66
67 71
72 76
77 83
84 90
91 100
101 111
112 125
126 142
143 166
167 200
201 250
251 333
334 500
501 1000

染色圖

#include<iostream>
#include<cstdio>
using namespace std;
const long long MOD=998244353;
long long pow_mod(long long a,long long n,long long m){
    long long ans=1;
    while(n){
        if(n&1){
            ans=(ans*a)%m;
        }
        a=(a*a)%m;
        n>>=1;
    }
    return ans;
}
int main(void){
 int T;
 long long n,nl,nr;
 cin>>T;
 long long Inv=pow_mod(2,MOD-2,MOD);
 while(T--){
  scanf("%lld%lld%lld",&n,&nl,&nr);
  //-2nt+nn-n+t(t+1)i
  //t*(t+1)i-2nt
  long long ans=((n*n-n)%MOD)*((nr-nl+1))%MOD;
  for(long long l=nl,r;l<=nr;l=r+1){
            r=min(n/(n/l),nr);
            long long t=n/l;
            ans=(ans+t*(t+1)%MOD*(r+l)%MOD*(r-l+1)%MOD*Inv)%MOD;
            ans=(ans-2*n*t%MOD*(r-l+1)%MOD+MOD)%MOD;
        }
        ans=ans*Inv%MOD;
        printf("%lld\n",ans);
 }
 return 0;
}