Time Limit: 10 Sec  Memory Limit: 128 MB

Description

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

Input

第一行两个数n,m。

Output

一个整数表示答案mod 19940417的值

Sample Input

3 4

Sample Output

1

样例说明
  
答案为(3 mod 1)*(4 mod 2)+(3 mod 1) * (4 mod 3)+(3 mod 1) * (4 mod 4) + (3 mod 2) * (4 mod 1) +

(3 mod 2) * (4 mod 3) + (3 mod 2) * (4 mod 4) + (3 mod 3) * (4 mod 1) + (3 mod 3) * (4 mod 2) +

(3 mod 3) * (4 mod 4) = 1

数据规模和约定
  对于100%的数据n,m<=10^9

这里是链接:【BZOJ】2956:模积和

这里是题解:

首先,暴力枚举将会很凄惨:O(nm)。早就 GG ( Time Limit Exceeded )了。

所以从公式入手:原公式是:

展开为:

观察式子:[n mod i](同理 m mod j)

根据mod的定义可以将上式写成这个样子:

[n/i]向下取整,就是C++中的整型(n/i),然后再乘以 i 就相当于下图灰色区域,

再用n减掉就能得到mod后的值。【下图模拟mod的转化

所以式子就可以简化为:

注意:因为题目中i!=j,所以当i、j 相同就直接减掉】

然而这样还是O(n^2)的复杂度。

所以继续化简:将第一个式子的Σ移动,使时间复杂度变为O(n)

这里是最终式子:

能够移动Σ的证明:设n-[n/i]*i为Xi,m-[m/j]*j为Yj。

(下图有点错误:n、m的值应该是不同的,但是n、m不同的证明也是这样的。)

这个是原式子展开的ans:

  这个是化简后展开的ans:

显然它们的ans值是相等的。那么,第一步化简式子已经完成了。

虽然移动Σ已经将复杂度降低到O(n),但很不幸的是依然过不了

考虑如何优化

低于O(n)的复杂度一般就三种:O(1)、O(logn)、O(√n)

注意最终式子,都有一个式子like this:[n/i](其中n为一个定值,i是从1到n的一个变量。)

但是这里有个很美妙的事情就是:

假如n==1000时,i在91~100之间,n/i的值都是为10的;

假如n==100000000时,i在9090910~10000000之间,n/i的值都是为10的。

这里因为在某一个区间中的值都是相等的。所以我们可以很愉快地利用分块的思想

那么怎么算这个块的大小呢?

假设有一个块里面的[n/i]的值都是为k,那么其区间就是:[(n/(k+1))+1,n/k].

推导:因为[n/i]是向下取整的,所以k*i<=n,n/k为定值,所以i<=n/k,但i一定也有下界,所以i>n/(k+1),

即i>=n/(k+1)+1.

注意:分块求值套用公式的时候需要除法,并不能先取模,然而不先取模会爆long long.

但好在除的数是固定的6,所以就直接在求平方和的时候,MOD开大6倍,最后再模回去就行了。

(其实反过来也就是网上普遍流传的3323403,是[mod/6]的值)

这里是代码:

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define LL long long
#define mod 19940417
#define MOD 119642502
using namespace std;
LL n,m,tmp1,tmp2,ans,ine; LL sum(LL x){
return x*(x+)/%mod;
} LL SUM(LL x){
return x*(x+)%MOD*(*x+)%MOD/;
} LL calc(LL x){
LL ans=x*x%mod;
for(LL i=;i<=x;i=tmp1+){
tmp1=x/(x/i);//(x/i)求k的值,n/k为块的上界
ans=(ans-(x/i)*(sum(tmp1)-sum(i-)+mod)%mod)%mod;
}
return ans%mod;
} int main(){
scanf("%lld %lld",&n,&m);
if(n>m) swap(n,m);
ans=calc(n)*calc(m)%mod;
//处理i,j相同情况
for(LL i=;i<=n;i=tmp2+){
tmp2=min(n/(n/i),m/(m/i));
ans=(ans-n*m%mod*(tmp2-i+)%mod
+m*(n/i)%mod*(sum(tmp2)-sum(i-)+mod)%mod
+n*(m/i)%mod*(sum(tmp2)-sum(i-)+mod)%mod
-(n/i)*(m/i)%mod*(SUM(tmp2)-SUM(i-)+mod)%mod+*mod)%mod;
}
printf("%lld",ans%mod);
return ;
}

【BZOJ】2956

梦想总是要有的,万一实现了呢?

【BZOJ】2956:模积和的更多相关文章

  1. BZOJ 2956 模积和 (数学推导+数论分块)

    手动博客搬家: 本文发表于20170223 16:47:26, 原地址https://blog.csdn.net/suncongbo/article/details/79354835 题目链接: ht ...

  2. BZOJ 2956 模积和

    题目链接:http://61.187.179.132/JudgeOnline/problem.php?id=2956 题意:给出n和m.计算: 思路: i64 n,m; i64 cal(i64 m,i ...

  3. [Bzoj 2956] 模积和 (整除分块)

    整除分块 一般形式:\(\sum_{i = 1}^n \lfloor \frac{n}{i} \rfloor * f(i)\). 需要一种高效求得函数 \(f(i)\) 的前缀和的方法,比如等差等比数 ...

  4. BZOJ 2956 模积和(分块)

    [题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=2956 [题目大意] 求∑∑((n%i)*(m%j))其中1<=i<=n,1 ...

  5. bzoj 2956: 模积和 ——数论

    Description 求∑∑((n mod i)*(m mod j))其中1<=i<=n,1<=j<=m,i≠j. Input 第一行两个数n,m. Output 一个整数表 ...

  6. 【BZOJ】2956: 模积和

    题意 求\(\sum_{i=1}^{n} \sum_{j=1}^{m} (n \ mod \ i)(m \ mod \ j)[i \neq j] \ mod \ 19940417\), \((n, m ...

  7. 「BZOJ 2956」模积和

    「BZOJ 2956」模积和 令 \(l=\min(n,m)\).这个 \(i\neq j\) 非常不优雅,所以我们考虑分开计算,即: \[\begin{aligned} &\sum_{i=1 ...

  8. BZOJ_2956_模积和_数学

    BZOJ_2956_模积和_数学 Description 求∑∑((n mod i)*(m mod j))其中1<=i<=n,1<=j<=m,i≠j. Input 第一行两个数 ...

  9. P2260 [清华集训2012]模积和

    P2260 [清华集训2012]模积和 整除分块+逆元 详细题解移步P2260题解板块 式子可以拆开分别求解,具体见题解 这里主要讲的是整除分块(数论分块)和mod不为素数时如何求逆元 整除分块:求Σ ...

随机推荐

  1. Maven学习(六)-----Maven仓库的详细介绍

    Maven仓库的详细介绍 在Maven中,任何一个依赖.插件或者项目构建的输出,都可以称之为构件.Maven在某个统一的位置存储所有项目的共享的构件,这个统一的位置,我们就称之为仓库.(仓库就是存放依 ...

  2. 一个IP可以登几个拼多多后台 拼多多如何推广营销

    想要在拼多多上开双店?担心一根网线会引起IP冲突?那么一根网线可以登几个拼多多后台呢?有数据显示,挂双店是没有问题的,甚至可以多开.那么解决了一根网线的事情,要怎么对自己的店铺进行营销推广呢?下面是小 ...

  3. C++自学第一课:函数

    此贴并非教学,主要是自学笔记,所述内容只是些许个人学习心得的记录和备查积累,难以保证观点正确,也不一定能坚持完成. 如不幸到访,可能耽误您的时间,也难及时回复,贴主先此致歉.如偶有所得,相逢有缘,幸甚 ...

  4. Selenium2+python自动化-xpath定位语法

    前言    在上一篇简单的介绍了用工具查看目标元素的xpath地址,工具查看比较死板,不够灵活,有时候直接复制粘贴会定位不到.这个时候就需要自己手动的去写xpath了,这一篇详细讲解xpath的一些语 ...

  5. MySQL双主复制

    原文发表于cu:2017-06-12 本文简单介绍MySQL双主复制原理及1个简单是双主复制验证. 一.MySQL双主复制原理 1. 双主复制原理 master-master复制的两台服务器,既是ma ...

  6. gets函数的完美替代

    众所周知 在C语言中scanf用来读取一行字符串时遇到空格或回车会停止 而若要读入一行带空格的字符串时 有些人会用gets来代替 然而,gets的最大问题在于:会读取超过数组长度上限个字符,而超出长度 ...

  7. apply新用法,最大值查找

    要找到数组中的最大或最小值,可以像下面这样使用apply() var values=[1,2,3,4,5,6,7,8]; var max = Math.max.apply(Math,values); ...

  8. Linux学习—导航文件系统

    与windows相同,Linux操作系统也是以被称之为分层目录结构的方式来组织文件的.这意味着文件是以树形结构的目录中进行组织的,该树形结构目录可能包含文件和其他目录.文件系统的第一个目录叫做根目录, ...

  9. DoItYourself!

    在杨老师的勉励下,我准备开始“自己”写程序.速度很慢,不过在写的过程中对于用到的几个函数更加熟悉.也尝试多学一点,学透一点.遇到不会的函数,语法不清楚的,还是会百度,不过会自己再敲一遍.重复下去. 下 ...

  10. “Hello World!”团队第六周的第二次会议

    今天是我们团队“Hello World!”团队第六周召开的第二次会议.博客内容: 一.会议时间 二.会议地点 三.会议成员 四.会议内容 五.todo list 六.会议照片 七.燃尽图 八.代码 一 ...