染色图（数论分块模板）

数论分块模板（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;
}