bzoj2956: 模积和

最新推荐文章于 2019-07-04 22:58:56 发布

*ACoder*

最新推荐文章于 2019-07-04 22:58:56 发布

阅读量353

点赞数

分类专栏： # 数学杂题

本文链接：https://blog.csdn.net/FSAHFGSADHSAKNDAS/article/details/71268796

版权

数学杂题专栏收录该内容

93 篇文章 0 订阅

订阅专栏

链接

　　http://www.lydsy.com/JudgeOnline/problem.php?id=2956

题解

　　这题难点在第一步，第一步想到后面就一路畅通了
　　先不考虑 $i\not =j$ 这个条件
　　

\sum i = 1 n \sum j = 1 m (n m o d i) (m m o d j) = \sum i = 1 n \sum j = 1 m (n - ⌊ n i ⌋) (m - ⌊ m j ⌋)

$\sum_{i=1}^n\sum_{j=1}^m(n\ mod\ i)(m\ mod\ j)\\ 　　=\sum_{i=1}^n\sum_{j=1}^m(n-\lfloor\frac{n}{i}\rfloor)(m-\lfloor\frac{m}{j}\rfloor)$
　　然后再减去

i=j $i=j$ 的情况（公式略）
　　画出柿子来搞一搞就好了。
　　 ~~博主真是懒~~

代码

//画柿子题 
#include <cstdio>
#include <algorithm>
#define ll long long
#define mod 19940417ll
using namespace std;
ll n, m, ans, tmp;
inline ll sum(ll a, ll b){return ((a+b)*(b-a+1)>>1)%mod;}
inline ll suan(ll x)
{
    if(x*(x+1)%3==0)return x*(x+1)/6%mod*(2*x+1)%mod;
    else return x*(x+1)/2%mod*(2*x+1)/3%mod;
}
inline ll sum2(ll a, ll b){return (suan(b)-suan(a-1))%mod;}
int main()
{
    ll i, last;
    scanf("%lld%lld",&n,&m);
    if(n>m)swap(n,m);
    for(i=1;i<=m;i=last+1)
    {
        last=m/(m/i);
        tmp=(tmp+m*(last-i+1)-m/i*sum(i,last))%mod;
    }
    for(i=1;i<=n;i=last+1)
    {
        last=n/(n/i);
        ans=(ans+(n*(last-i+1)-n/i*sum(i,last))%mod*tmp)%mod;
    }
    tmp=n*m%mod*n%mod;
    for(i=1;i<=n;i=last+1)
    {
        last=min(n/(n/i),m/(m/i));
        tmp=(tmp-n*(m/i)%mod*sum(i,last)-m*(n/i)%mod*sum(i,last))%mod;
        tmp=(tmp+(n/i)*(m/i)%mod*sum2(i,last))%mod;
    }
    ans=(ans-tmp)%mod;
    printf("%lld",(ans+mod)%mod);
    return 0;
}