bzoj2956 模积和

最新推荐文章于 2018-04-20 13:26:00 发布

Euryale_

最新推荐文章于 2018-04-20 13:26:00 发布

阅读量483

点赞数

分类专栏：数论文章标签：数论分块

本文链接：https://blog.csdn.net/u010336344/article/details/65935910

版权

数论专栏收录该内容

25 篇文章 0 订阅

订阅专栏

题目链接：bzoj2956
题目大意：
求 $\sum_{i=1}^{n}\sum_{j=1}^{m}(n\mod i)\times (m\mod j),i≠j$

题解：
数论+分块
看清题啊 $i≠j$ ！
所以式子就变成：

\sum i = 1 n \sum j = 1 m (n mod i) \times (m mod j) - \sum i = 1 m i n (n, m) (n mod i) \times (m mod j)

$\sum_{i=1}^{n}\sum_{j=1}^{m}(n\mod i)\times (m\mod j)-\sum_{i=1}^{min(n,m)}(n\mod i)\times (m\mod j)$
像 bzoj1257那样处理
于是就变成了：

\sum i = 1 n \sum j = 1 m (n - i \times ⌊ n i ⌋) \times (m - j \times ⌊ m j ⌋) - \sum i = 1 m i n (n, m) (n - i \times ⌊ n i ⌋) \times (m - j \times ⌊ m j ⌋)

$\sum_{i=1}^{n}\sum_{j=1}^{m}(n-i\times \lfloor\frac{n}{i}\rfloor)\times (m-j\times \lfloor\frac{m}{j}\rfloor)-\sum_{i=1}^{min(n,m)}(n-i\times \lfloor\frac{n}{i}\rfloor)\times (m-j\times \lfloor\frac{m}{j}\rfloor)$
傻傻的我把两边都拆开直接搞了。
于是对于

∑ni=1∑mj=1i×j×⌊ni⌋×⌊mj⌋ $\sum_{i=1}^n\sum_{j=1}^{m}i\times j\times\lfloor\frac{n}{i}\rfloor\times\lfloor \frac{m}{j}\rfloor$ 表示打了个

O(n√m−−√) $O(\sqrt{n}\sqrt{m})$ 的…诶好蠢啊ww
应该把式子化成：

\sum i = 1 n (n - i \times ⌊ n i ⌋) \sum j = 1 m (m - j \times ⌊ m j ⌋) - \sum i = 1 m i n (n, m) (n - i \times ⌊ n i ⌋) \times (m - j \times ⌊ m j ⌋))

$\sum_{i=1}^{n}(n-i\times \lfloor\frac{n}{i}\rfloor)\sum_{j=1}^{m} (m-j\times \lfloor\frac{m}{j}\rfloor)-\sum_{i=1}^{min(n,m)}(n-i\times \lfloor\frac{n}{i}\rfloor)\times (m-j\times \lfloor\frac{m}{j}\rfloor))$
这样前面的就可以把

n $n$ 、

m $m$ 提出来之后像bzoj1257那样分块搞，后半式的话拆开之后也是类似的，用分块来做(像 bzoj2301那样)。

总结：
有时可以考虑内外层求和对调啊
还有分块的套路要熟悉。要有意识看到 $\lfloor\frac{n}{i}\rfloor$ 、 $\lfloor\frac{n}{i}\rfloor\times\lfloor\frac{m}{j}\rfloor$ 之后用分块啊

话说，这题mod得我真恶心，因为mod得不够多WA了两三次！

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;

const LL ny=3323403;
const LL mod=19940417;
LL mymin(LL x,LL y){return (x<y)?x:y;}
LL get(LL n,LL lim)
{
    LL ret=0,r;
    for (LL i=1;i<=lim;i=r+1)
    {
        LL x=n/i;r=mymin(lim,n/x);
        ret=(ret+x%mod*((r-i+1)*(r+i)/2%mod))%mod;
    }return ret;
}
int main()
{
    //freopen("a.in","r",stdin);
    //freopen("a.out","w",stdout);
    LL n,m,ans;
    scanf("%lld%lld",&n,&m);
    ans=(((n*n)%mod-get(n,n)+mod)%mod)*(((m*m)%mod-get(m,m)+mod)%mod)%mod;

    LL lim=mymin(n,m),x,y,r,lx,ly;
    ans=(ans-(lim*n%mod*m%mod)+mod)%mod;
    for (LL i=1;i<=lim;i=r+1)
    {
        x=n/i;y=m/i;
        r=mymin(n/x,m/y);

        ans=(ans+m*x%mod*((r-i+1)*(r+i)/2%mod))%mod;
        ans=(ans+n*y%mod*((r-i+1)*(r+i)/2%mod))%mod;

        lx=((r*(r+1)%mod*(2*r+1)%mod)*ny)%mod;
        ly=(((i-1)*i%mod*(2*i-1)%mod)*ny)%mod;
        ans=(ans-(x*y%mod*((lx-ly+mod)%mod))%mod)%mod;

        ans=(ans+mod)%mod;
    }

    while (ans<0) ans+=mod;
    printf("%lld\n",ans%mod);
    return 0;
}