数论：BZOJ2956 (模积和)

题意：求∑∑((n mod i)*(m mod j))其中1<=i<=n,1<=j<=m,i≠j。然后mod 19940417

本题坑我太久啊，思路：

∑∑((n mod i) * (m mod j)) 1<=i<=n, 1<=j<=m, i≠j= ∑(n mod i) * ∑(m mod i) - ∑((n mod i) * (m mod i))= ∑(n-[n/i]*i) * ∑(m-[m/i]*i) - ∑(nm-([n/i]+[m/i])i+[n/i][m/i]*i*i)

[cpp]  view plain copy

1. #include <iostream>  
2. #include <algorithm>  
3.   
4. using namespace std;  
5.   
6. typedef long long LL;  
7.   
8. const LL MOD = 19940417;  
9.   
10. LL sum(LL n)  
11. {  
12.     return n*(n+1)%MOD*(2*n+1)%MOD*3323403%MOD;  
13. }  
14.   
15. LL Solve(LL m, LL n)  
16. {  
17.     LL ans=0,i,last;  
18.     for(i=1;i<=m;i=last+1)  
19.     {  
20.         last=min(m,n/(n/i));  
21.         ans += (n/i)*(i+last)%MOD*(last-i+1)%MOD*9970209%MOD;  
22.         ans %= MOD;  
23.     }  
24.     return ans;  
25. }  
26.   
27. int main()  
28. {  
29.     LL n, m;  
30.     cin>>n>>m;  
31.     if(n<m) swap(n, m);  
32.     LL ans=(n*n-Solve(n,n))%MOD*((m*m-Solve(m,m))%MOD);  
33.     ans+=-m*m%MOD*n%MOD+Solve(m, n)*m%MOD+Solve(m, m)*n%MOD;  
34.     ans%=MOD;  
35.     for(LL i=1,last;i<=m;i=last+1)  
36.     {  
37.         last=min(m,min(n/(n/i),m/(m/i)));  
38.         ans += -(n/i)*(m/i)%MOD*((sum(last)-sum(i-1))%MOD)%MOD;  
39.         ans %= MOD;  
40.     }  
41.     cout<<(ans%MOD+MOD)%MOD<<endl;  
42.     return 0;  
43. }