#include using namespace std; typedef signed long long ll; #undef _P #define _P(...) (void)printf(__VA_ARGS__) #define FOR(x,to) for(x=0;x<(to);x++) #define FORR(x,arr) for(auto& x:arr) #define ITR(x,c) for(__typeof(c.begin()) x=c.begin();x!=c.end();x++) #define ALL(a) (a.begin()),(a.end()) #define ZERO(a) memset(a,0,sizeof(a)) #define MINUS(a) memset(a,0xff,sizeof(a)) //------------------------------------------------------- int P,Q; int L,R; ll A[101010],S[101010]; int nex[101010]; ll NS[101010]; void solve() { int i,j,k,l,r,x,y; string s; cin>>P>>Q; for(i=1;i<=100001;i++) { A[i]=P%i; S[i]=S[i-1]+A[i]; nex[i]=P/(i+1)+1; } for(i=1;i<=100000;i++) { L=nex[i]; R=nex[i-1]+1; NS[i]=1LL*P*(R-L+1)-1LL*i*(R+L)*(R-L+1)/2; } while(Q--) { cin>>L>>R; ll ret=0; if(R<=100000) { ret=S[R]-S[L-1]; } else { if(L<=100000) { ret+=S[100000]-S[L-1]; L=100001; } x=P/L; y=P/R; while(x!=y) { if(L==nex[x]) { ret+=NS[x]; } else { int ne=nex[x-1]+1; ret+=1LL*P*(ne-L+1)-1LL*x*(ne+L)*(ne-L+1)/2; } x--; L=nex[x]; } if(L<=R) { int x=P/R; ret+=1LL*P*(R-L+1)-1LL*x*(R+L)*(R-L+1)/2; } } cout<