结论
答案为
(
a
+
b
∑
k
=
1
a
k
)
×
∑
r
=
1
b
−
1
r
(a+b\sum \limits_{k=1} \limits^ak) \times \sum \limits_{r=1} \limits^{b-1} r
(a+bk=1∑ak)×r=1∑b−1r
还是建议自己手推一遍。过程见下。
题意
对于给定的 a a a 和 b b b,数 x x x 能成为“好数”的条件:
记 q = ⌊ x b ⌋ q=\left\lfloor\frac{x}{b}\right\rfloor q=⌊bx⌋, r = x m o d b r=x \bmod b r=xmodb,分别代表 x x x 除以 b b b 的商和余数。
- r ≠ 0 r\neq 0 r=0
- q = k r q=kr q=kr, k k k 为整数。
- 1 ≤ q r ≤ a 1 \leq \frac{q}{r} \leq a 1≤rq≤a
求所有好数的和。
推导
由条件的第二条可知, q = k r q=kr q=kr, k k k 为整数。
由带余除法的定义得:
x
=
q
b
+
r
=
k
r
b
+
r
=
(
k
b
+
1
)
r
x=qb+r=krb+r=(kb+1)r
x=qb+r=krb+r=(kb+1)r
将 q = k r q=kr q=kr 代入条件三得: 1 ≤ k ≤ a 1 \leq k \leq a 1≤k≤a,且根据余数的定义得 1 ≤ r 1 \leq r 1≤r < b < b <b。
有了范围,可以统计范围内的答案了。
记 s u m sum sum 为所有满足条件的 x x x 的和。则:
s u m = ∑ [ ( k b + 1 ) r ] sum=\sum\left[ (kb+1)r \right] sum=∑[(kb+1)r]
= ∑ k = 1 a ( k b + 1 ) × ∑ r = 1 b − 1 r =\sum\limits_{k=1}^{a} (kb+1) \times \sum\limits_{r=1}^{b-1} r =k=1∑a(kb+1)×r=1∑b−1r
= ( a + ∑ k = 1 a k b ) × ∑ r = 1 b − 1 r =(a+\sum\limits_{k=1}^{a} kb) \times \sum\limits_{r=1}^{b-1} r =(a+k=1∑akb)×r=1∑b−1r
= ( a + b ∑ k = 1 a k ) × ∑ r = 1 b − 1 r =(a+b\sum\limits_{k=1}^{a} k) \times \sum\limits_{r=1}^{b-1} r =(a+bk=1∑ak)×r=1∑b−1r
推到这样为止,剩下的就不难了。
代码
#include<bits/stdc++.h>
using namespace std;
const long long mod=1e9+7;
long long a,b;
signed main()
{
scanf("%lld%lld",&a,&b);
long long sigma1=a*(a+1)/2%mod;
long long sigma2=b*(b-1)/2%mod;
long long ans=(sigma1*b%mod+a)%mod*sigma2%mod;
printf("%lld\n",ans);
return 0;
}