http://codeforces.com/problemset/problem/235/E (题目链接)

题意

  给出${a,b,c}$,求${\sum_{i=1}^a\sum_{j=1}^b\sum_{k=1}^cd(ijk)}$

extra

  有这样一个公式,就是约数个数和那道题的推广吧。$${\sum_{i=1}^a\sum_{j=1}^b\sum_{k=1}^cd(ijk)=\sum_{i=1}^a\sum_{j=1}^b\sum_{k=1}^c[gcd(i,j)=gcd(i,k)=gcd(j,k)=1]\lfloor\frac{a}{i}\rfloor\lfloor\frac{b}{j}\rfloor\lfloor\frac{c}{k}\rfloor}$$

  然并卵,这个式子我不会证也不会用。。还是直接推吧。

Solution

  莫比乌斯反演。

\begin{aligned}  \sum_{i=1}^a\sum_{j=1}^b\sum_{k=1}^cd(ijk)=\sum_{i=1}^{ab}f(i)\sum_{j=1}^cd(ij)   \end{aligned}

  其中${f(n)=\sum_{i=1}^a\sum_{j=1}^b[ab=n]}$。

\begin{aligned}     & \sum_{i=1}^{ab}f(i)\sum_{j=1}^cd(ij)  \\  =&\sum_{i=1}^{ab}f(i)\sum_{j=1}^c\sum_{u|i}\sum_{v|j}[gcd(u,v)=1]  \\  =&\sum_{u=1}^{ab}\sum_{v=1}^c[gcd(u,v)=1]\sum_{i=1}^{\lfloor{ab/u}\rfloor}f(iu)\lfloor\frac{c}{v}\rfloor     \end{aligned}

  我们令${S(n)=\sum_{i=1}^{\lfloor{ab/n}\rfloor}f(in)}$,并把变量${u,v}$换成${i,j}$,因为${u,v}$看起来太丑了→_→。

\begin{aligned}     & \sum_{i=1}^{ab}\sum_{j=1}^c[gcd(i,j)=1]S(i)\lfloor\frac{c}{j}\rfloor  \\  =&\sum_{i=1}^{ab}\sum_{j=1}^c\sum_{t|i,t|j}μ(t)S(i)\lfloor\frac{c}{j}\rfloor  \\  =&\sum_{t=1}^cμ(t)\sum_{i=1}^{\lfloor{ab/t}\rfloor}S(it)\sum_{j=1}^{\lfloor{c/t}\rfloor}\lfloor\frac{c}{jt}\rfloor     \end{aligned}

  看到这个式子是不是感到了满满的套路,我们故技重施,令${P(n)=\sum_{i=1}^{\lfloor{ab/n}\rfloor}S(in)}$,${Q(n)=\sum_{i=1}^n\lfloor\frac{n}{i}\rfloor}$

\begin{aligned}     \sum_{t=1}^cμ(t)P(t)Q(\lfloor\frac{c}{t}\rfloor)     \end{aligned}

  这样,我们${O(ab)}$的求出${f}$,之后按顺序${O(ab\log ab)}$的求出${S}$和${P}$,然后${O(c^2)}$的求出${Q}$,当然如果想分段也可以分段。最后只需要从${1}$枚举到${c}$就可以得出答案了。

细节

  时限卡得有点紧,你需要常数优化→_→

代码

#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<queue>
#define LL long long
#define inf 2147483640
#define MOD (1ll<<30)
#define Pi acos(-1.0)
#define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
using namespace std; const int maxn=4000010,maxm=2010;
LL s[maxn],f[maxn],P[maxm],Q[maxm];
int a,b,c,p[maxm],vis[maxm],mu[maxm]; void inc(LL &a,LL b) {
a+=b;if (a>MOD) a-=MOD;
}
int main() {
scanf("%d%d%d",&a,&b,&c);
mu[1]=1;
for (int i=2;i<=c;i++) {
if (!vis[i]) p[++p[0]]=i,mu[i]=-1;
for (int j=1;j<=p[0] && i*p[j]<=c;j++) {
vis[i*p[j]]=1;
if (i%p[j]==0) {mu[i*p[j]]=0;break;}
else mu[i*p[j]]=-mu[i];
}
}
for (int i=1;i<=a;i++)
for (int j=1;j<=b;j++) f[i*j]++;
for (int i=1;i<=a*b;i++)
for (int j=1;j<=a*b/i;j++) inc(s[i],f[i*j]);
for (int i=1;i<=c;i++)
for (int j=1;j<=a*b/i;j++) inc(P[i],s[i*j]);
for (int i=1;i<=c;i++)
for (int j=1;j<=i;j++) inc(Q[i],i/j);
LL ans=0;
for (int i=1;i<=c;i++)
ans=(ans+mu[i]*P[i]*Q[c/i]%MOD)%MOD;
printf("%lld",(ans+MOD)%MOD);
return 0;
}

【codeforces 235E】 Number Challenge的更多相关文章

  1. 【Codeforces 466C】Number of Ways

    [链接] 我是链接,点我呀:) [题意] 让你把数组分成3个连续的部分 每个部分的和要一样 问你有多少种分法 [题解] 先处理出来num[i] 表示i..n这里面有多少个j 满足aft[j] = af ...

  2. 【CodeForces 660D】Number of Parallelograms(n个点所能组成的最多平行四边形数量)

    You are given n points on a plane. All the points are distinct and no three of them lie on the same ...

  3. 【codeforces 415D】Mashmokh and ACM(普通dp)

    [codeforces 415D]Mashmokh and ACM 题意:美丽数列定义:对于数列中的每一个i都满足:arr[i+1]%arr[i]==0 输入n,k(1<=n,k<=200 ...

  4. 【codeforces 514A】Chewbaсca and Number

    [题目链接]:http://codeforces.com/contest/514/problem/A [题意] 允许你把每个数字翻转 ->x变成9-x 然后问你能够变成的最小的数字是什么; 不能 ...

  5. 【codeforces 805D】Minimum number of steps

    [题目链接]:http://codeforces.com/contest/805/problem/D [题意] 给你一个字符串; 里面只包括a和b; 让你把里面的"ab"子串全都去 ...

  6. 【codeforces 546D】Soldier and Number Game

    time limit per test3 seconds memory limit per test256 megabytes inputstandard input outputstandard o ...

  7. 【codeforces 768C】Jon Snow and his Favourite Number

    [题目链接]:http://codeforces.com/contest/768/problem/C [题意] 给你n个数字; 让你每次把这n个数字排序; 然后对奇数位的数字进行异或操作,然后对新生成 ...

  8. 【26.67%】【codeforces 596C】Wilbur and Points

    time limit per test2 seconds memory limit per test256 megabytes inputstandard input outputstandard o ...

  9. 【codeforces 766B】Mahmoud and a Triangle

    time limit per test2 seconds memory limit per test256 megabytes inputstandard input outputstandard o ...

随机推荐

  1. git解决代码提交冲突

    树冲突文件名修改造成的冲突,称为树冲突.比如,A同事把文件改名为A.C,B同事把同一个文件改名为B.C,那么B同事将这两个commit合并时,会产生冲突.如果最终确定用B同事的文件名,那么解决办法如下 ...

  2. 如何判断Map中的key或value是什么类型

    在上班写工具类时,遇到了一个问题,将xml文件的节点都放入map容器中时,map的value也是一个map,导致取map的value时,需要判断这个value的数据类型,用到了一下说的这些知识: 对于 ...

  3. 笨办法学Python - 习题4: Variables and Names

    1.习题 4: 变量(variable)和命名 学习目标:了解Python中变量的定义,学习定义简明易记的变量名 变量:变量是存储内存中的值,就是每定义一个变量就会在内存中开辟一个空间.基于变量的类型 ...

  4. logout命令详解

    基础命令学习目录首页 logout指令让用户退出系统,其功能和login指令相互对应.语法 logout

  5. prototype原型(待完善)

    模式:prototype  解决向量的深浅克隆 #pragma once #ifndef _PROTOTYPE_H_ #define _PROTOTYPE_H_ class Prototype{ pu ...

  6. 每天学一点easyui①

    引入js和css文件 <script type="text/javascript" src="js/jquery-easyui-1.4.3/jquery.min.j ...

  7. java 框架 面试

    Java—SSH(MVC)1. 谈谈你mvc的理解MVC是Model—View—Controler的简称.即模型—视图—控制器.MVC是一种设计模式,它强制性的把应用程序的输入.处理和输出分开.MVC ...

  8. php之 常用的 流程管理

    1.流程管理的用法是什么样的? 2.怎么发起想要的流程? 3.审批的人要是怎么审批通过? 4.流程审核是不是要挨个走过? 一.要有数据库的内容的 肯定会有表的,首先就是用户表了,然后就是流程表,用户编 ...

  9. AttributeError: module ‘tensorflow.python.ops.nn’ has no attribute ‘leaky_relu’

    #AttributeError: module 'tensorflow.python.ops.nn' has no attribute 'leaky_relu' 的原因主要是版本的问题 解决方法是更新 ...

  10. ListViewAnimations使用时报错NoClassDefFoundError: com.nineoldandroids.animation.Animator

    见 https://github.com/nhaarman/ListViewAnimations/issues/294 解决: Add this to your dependencies in you ...