BZOJ 2956 模积和

题目链接：模积和

　　这种有模运算的题一般都要把取模运算给去掉，改成除法取下整的模式。即\(a\bmod b=a-\lfloor \frac{a}{b} \rfloor b\)

　　然后我们先把\(i=j\)的也统计进答案，最后再减去即可。接下来就是推式子时间(\(n \le m\))：

\begin{aligned}
&\sum_{i=1}^n\sum_{j=1}^m(n \bmod i)(m \bmod j) \\
=&\sum_{i=1}^n(n \bmod i)\sum_{j=1}^m(m \bmod j) \\
=&\sum_{i=1}^n(n-i\lfloor \frac{n}{i} \rfloor)\sum_{j=1}^m(m-j\lfloor \frac{m}{j} \rfloor) \\
=&(n^2-\sum_{i=1}^ni\lfloor \frac{n}{i} \rfloor)(m^2-\sum_{i=1}^mi\lfloor \frac{m}{i} \rfloor)
\end{aligned}

\begin{aligned}
&\sum_{i=1}^n(n \bmod i)(m \bmod i) \\
=&\sum_{i=1}^n(n-i\lfloor \frac{n}{i} \rfloor)(m-i\lfloor \frac{m}{i} \rfloor) \\
=&n^2m-m\sum_{i=1}^ni\lfloor \frac{n}{i} \rfloor-n\sum_{i=1}^ni\lfloor \frac{m}{i} \rfloor+\sum_{i=1}^n\lfloor \frac{n}{i} \rfloor\lfloor \frac{m}{i} \rfloor i^2
\end{aligned}

　　然后分块计算即可。上下两个式子一减就是答案。

　　下面贴代码：

#include<cstdio>
#define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)
#define mod 19940417
#define llg long long
#define min(a,b) (a<b?a:b)const int n6=7776;
int n,m,a[3];llg S(int x){return (1ll*x*(x+1)>>1)%mod;}
llg S2(int x){a[0]=x,a[1]=x+1,a[2]=2*x+1;for(int i=0;i<3;i++) if(a[i]%2==0){a[i]/=2;break;}for(int i=0;i<3;i++) if(a[i]%3==0){a[i]/=3;break;}return 1ll*a[0]*a[1]%mod*a[2]%mod;
}llg work(int x,int y){llg now=0;for(int i=1,nt;i<=x;i=nt+1)nt=min(y/(y/i),x),now+=1ll*(y/i)*(S(nt)-S(i-1));return now%mod;
}int main(){File("a");scanf("%d %d",&n,&m);if(n>m) n^=m^=n^=m;llg s1=work(n,n),ans;ans=(1ll*n*n-s1)%mod*(1ll*m*m%mod-work(m,m))%mod;ans+=s1*m+work(n,m)*n; ans%=mod;ans-=1ll*n*n%mod*m; ans%=mod;for(int i=1,nt;i<=n;i=nt+1){nt=min(n/(n/i),m/(m/i));ans-=(S2(nt)-S2(i-1))*(m/i)%mod*(n/i);ans%=mod;}printf("%lld",(ans+mod)%mod);return 0;
}

转载于:https://www.cnblogs.com/lcf-2000/p/6811864.html

BZOJ 2956 模积和相关推荐

BZOJ 2956 模积和（分块）
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=2956 [题目大意] 求∑∑((n%i)*(m%j))其中1<=i<=n,1 ...
BZOJ2956 (模积和)
题目:2956: 模积和题意:求∑∑((n mod i)*(m mod j))其中1<=i<=n,1<=j<=m,i≠j.然后mod 19940417 本题坑我太久啊,思路: ...
P2260 [清华集训2012]模积和，P2834 能力测验（二维除法分块）
P2260 [清华集训2012]模积和推导过程我们假定n<=mn <= mn<=m ∑i=1n∑j=1m(nmodi)(mmodj),i≠j\sum_{i = 1} ^{n} ...
数论分块专题复习（余数求和+模积和+Ice Rain+The Fool）
文章目录前提知识复习 T1:余数求和 title solution code T2:Ice Rain title solution code T3:The Fool title solution c ...
一句话题解（20170801~20170125）
8.1 bzoj 4720 noip2016 换教室 floyd预处理+期望(薛定谔的猫) bzoj 4318 OSU! 三次函数期望值从一次.二次推得 8.2 bzoj 1076 状压+期望DP ...
[颓废史]蒟蒻的刷题记录
QAQ蒟蒻一枚,其实我就是来提供水题库的. 以下记录从2016年开始. 1.1 1227: [SDOI2009]虔诚的墓主人树状数组+离散化 3132: 上帝造题的七分钟树状数组二维区间加减+查 ...
2019.4.summary
2019.4.1 BZOJ1061: [Noi2008]志愿者招募真心有点难QAQ https://www.byvoid.com/zhs/blog/noi-2008-employee 看void爷的 ...
数论与数学专题练习（一）(201802~201805)
数论与数学专题练习(一)(201802~201805) 手动博客搬家: 本文发表于20180224 00:33:28, 原地址https://blog.csdn.net/suncongbo/artic ...
银联高校极客挑战赛初赛第一场、第二场合集（3+3题）
第一场 A. 码队女朋友的王者之路测试用例测试1 4 5 2 2 00101 5 1 2 11110 6 2 3 000000 7 6 6 0000001 测试2 1 1 1 1 0 测试3 1 ...

BZOJ 2956 模积和

BZOJ 2956 模积和相关推荐

最新文章

热门文章