模积和(bzoj 2956)

求∑∑((n mod i)*(m mod j))其中1<=i<=n,1<=j<=m,i≠j。

对于100%的数据n,m<=10^9。

 

 

      

对于而言,他的个数很小,最大是63244,因此可以枚举分块确定区域 [ l , n/(n/l) ] 

 

c++ code:

#include <bits/stdc++.h>
 
using namespace std;
const int mod=19940417;
const int ni=3323403;
typedef long long ll;
 
ll sum(ll l,ll r)
{
    return (r-l+1)*(l+r)/2%mod;
}
ll cal(ll n)
{
    ll l,r,ans=n*n%mod;
    l=1;
    while(l<=n)
    {
        r=n/(n/l);
        ans=(ans-(n/l)*sum(l,r)%mod+mod)%mod;
        l=r+1;
    }
    return ans;
}
 
ll pre(ll l,ll r)
{
    return (r*(r+1)%mod*(2*r+1)%mod*ni%mod-l*(l+1)%mod*(2*l+1)%mod*ni%mod+mod)%mod;
}
 
int main()
{
    ll m,n;
    scanf("%lld%lld",&m,&n);
    if(m<n) swap(m,n);
    ll ans=(cal(n)*cal(m)%mod-n*n%mod*m%mod+mod)%mod;
    ll l=1,r;
    while(l<=n)
    {
        ll x,y;
        r=min(n/(n/l),m/(m/l));
        x=n/l;y=m/l;
        ans=(ans+sum(l,r)*(n*y%mod+m*x%mod)%mod-x*y%mod*pre(l-1,r)%mod+mod)%mod;
        l=r+1;
    }
    printf("%lld\n",ans);
    return 0;
}

 

转载于:https://www.cnblogs.com/lemon-jade/p/9028615.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值