2956: 模积和乘法逆元取模

最新推荐文章于 2024-09-16 15:57:29 发布

ws_fqk

最新推荐文章于 2024-09-16 15:57:29 发布

阅读量478

点赞数

CC 4.0 BY-SA版权

分类专栏： My Code

本文链接：https://blog.youkuaiyun.com/Phenix_2015/article/details/50672550

My Code 专栏收录该内容

349 篇文章

订阅专栏

本文解析了一道BZOJ题目，并提供了一份详细的代码实现。通过使用longlong类型避免了数据处理中的类型转换问题，针对不同的计算需求，设计了多个辅助函数进行模块化处理，有效解决了大数据场景下的计算难题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

建议直接开long long，否则各种强制转换类型可能会出问题。。我因为这调了好久，小数据拍的飞起大数据各种跪。。
公式见hzwer
bzoj题目描述坑。。其实还是挺简单的，分开计算就好了。。

#include<iostream>
#include<cstdio>
#define ll long long
#define p 19940417
#define ine2 9970209
#define ine6 3323403
using namespace std;
ll n,m;
ll calc(ll k,ll n)
{
    ll tmp=0;
    for (ll i=1,pos=0;i<=k;i=pos+1)
    {
        pos=min(n/(n/i),k);
        (tmp+=(n/i)%p*(((pos+1)*(pos)%p*ine2%p-(i-1)*i%p*ine2%p+p)%p)%p)%=p;
    }
    return (tmp+p)%p;
}
ll calc0(ll n,ll m)
{
    ll tmp=0;
    for (ll i=1,pos=0;i<=n;i=pos+1)
    {
        pos=min(n/(n/i),m/(m/i));
        (tmp+=(n/i)*(m/i)%p*((pos*(pos+1)%p*(pos*2+1)%p*ine6%p-(i-1)*i%p*(i*2-1)%p*ine6%p+p)%p)%p)%=p;
    }
    return (tmp+p)%p;
}
int main()
{
    scanf("%lld%lld",&n,&m);
    if (n>m) swap(n,m);
    ll t1=calc(n,n),t2=calc(m,m),t3=calc(n,m),t4=calc0(n,m);
    ll ans=((((ll)n*n%p-t1+p)%p)*(((ll)m*m%p-t2+p)%p)%p-((ll)n*n%p*m%p-(ll)t1*m%p-(ll)t3*n%p+t4%p+p*10)%p+p)%p;
    cout << ans << endl;
    return 0;
}