【BZOJ3529】数表(莫比乌斯反演,树状数组)

题解

首先不管\(A\)的范围的限制

要求的东西是

\[\sum_{i=1}^n\sum_{j=1}^m\sigma(gcd(i,j))
\]

其中\(\sigma(x)\)表示\(x\)的约数之和

约数之和是一个积性函数,可以线性筛

具体的做法请参考皮皮亮的Blog

根据常见的套路

把\(gcd\)给提出来

\[\sum_{d=1}^n\sigma(d)\sum_{i=1}^n\sum_{j=1}^m[gcd(i,j)=d]
\]

后面那个东西不用说了吧。。

我原来已经推过很多遍啦

可以参考原来推过的式子

所以我直接写啦

\[\sum_{d=1}^n\sigma(d)\sum_{i=1}^{n/d}\mu(i)[\frac{n}{id}][\frac{m}{id}]
\]

还是一样的,设\(T=id\)

然后提出来

\[\sum_{T=1}^n[\frac{n}{T}][\frac{m}{T}]\sum_{d|T}\sigma(d)\mu(\frac{T}{d})
\]

如果没有\(A\)的限制

这就可以直接筛后面的那个狄利克雷卷积

然后\(O(\sqrt n)\)的回答每一组询问

现在有了\(A\)的限制,貌似筛出来也没有什么用了

考虑到\(n,m<=10^5\)应该可以乱搞了

所以干脆就先筛出\(\sigma\)和\(\mu\)

然后,这个限制还是很不好搞???

范围这么小

那我就暴力呀

离线询问后按照\(A\)排序

然后把\(\sigma\)排序

依次插到相对应的位置

因为要求前缀,就搞一个树状数组来算前缀和就行啦

最后是取膜,为啥要用\(long\ long\)

\(int\)自然溢出就行啦

暴力出奇迹

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define MAX 100000
inline int read()
{
int x=0,t=1;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=-1,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return x*t;
}
bool zs[MAX+1];
int pri[MAX],tot;
int mu[MAX+1],sig[MAX+1],sumd[MAX+1],powd[MAX+1];
int T,N,ans[MAX];
struct Ask{int n,m,A,id;}q[MAX];
int p[MAX+1],c[MAX+1];
bool operator<(Ask a,Ask b){return a.A<b.A;}
void pre(int N)
{
zs[1]=true;sig[1]=mu[1]=p[1]=1;
for(int i=2;i<=N;++i)
{
p[i]=i;
if(!zs[i])pri[++tot]=i,mu[i]=-1,sig[i]=i+1,sumd[i]=i+1,powd[i]=i;
for(int j=1;j<=tot&&i*pri[j]<=N;++j)
{
zs[i*pri[j]]=true;
if(i%pri[j])
{
mu[i*pri[j]]=-mu[i];
sig[i*pri[j]]=sig[i]*sig[pri[j]];
sumd[i*pri[j]]=1+pri[j];
powd[i*pri[j]]=pri[j];
}
else
{
powd[i*pri[j]]=powd[i]*pri[j];
sumd[i*pri[j]]=sumd[i]+powd[i*pri[j]];
sig[i*pri[j]]=sig[i]/sumd[i]*sumd[i*pri[j]];
break;
}
}
}
}
bool cmp(int a,int b){return sig[a]<sig[b];}
inline int lowbit(int x){return x&(-x);}
void Add(int x,int w){while(x<=N)c[x]+=w,x+=lowbit(x);}
int Query(int x){int ret=0;while(x)ret+=c[x],x-=lowbit(x);return ret;}
int Solve(int x,int y)
{
int now=0,lst=0,ret=0;
for(int i=1,j;i<=x;i=j+1)
{
j=min(x/(x/i),y/(y/i));
now=Query(j);
ret+=(x/i)*(y/i)*(now-lst);
lst=now;
}
return ret;
}
int main()
{
T=read();
for(int i=1;i<=T;++i)
{
q[i].id=i,q[i].n=read(),q[i].m=read(),q[i].A=read();
if(q[i].n>q[i].m)swap(q[i].n,q[i].m);
N=max(N,q[i].n);
}
pre(N);
sort(&q[1],&q[T+1]);sort(&p[1],&p[N+1],cmp);
for(int i=1,j=1;i<=T;++i)
{
for(;j<=N&&sig[p[j]]<=q[i].A;++j)
for(int k=p[j];k<=N;k+=p[j])
if(mu[k/p[j]])Add(k,sig[p[j]]*mu[k/p[j]]);
ans[q[i].id]=Solve(q[i].n,q[i].m);
}
for(int i=1;i<=T;++i)
{
if(ans[i]<0)ans[i]+=2147483647,ans[i]++;
printf("%d\n",ans[i]);
}
return 0;
}

【BZOJ3529】数表(莫比乌斯反演,树状数组)的更多相关文章

  1. 【BZOJ3529】[Sdoi2014]数表 莫比乌斯反演+树状数组

    [BZOJ3529][Sdoi2014]数表 Description 有一张N×m的数表,其第i行第j列(1 < =i < =礼,1 < =j < =m)的数值为能同时整除i和 ...

  2. BZOJ3529: [Sdoi2014]数表(莫比乌斯反演 树状数组)

    题意 题目链接 Sol 首先不考虑\(a\)的限制 我们要求的是 \[\sum_{i = 1}^n \sum_{j = 1}^m \sigma(gcd(i, j))\] 用常规的套路可以化到这个形式 ...

  3. BZOJ 3259 [Sdoi2014]数表 (莫比乌斯反演 + 树状数组)

    3529: [Sdoi2014]数表 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 2321  Solved: 1187[Submit][Status ...

  4. BZOJ 3529 [Sdoi2014]数表 (莫比乌斯反演+树状数组+离线)

    题目大意:有一张$n*m$的数表,第$i$行第$j$列的数是同时能整除$i,j$的所有数之和,求数表内所有不大于A的数之和 先是看错题了...接着看对题了发现不会做了...刚了大半个下午无果 看了Po ...

  5. BZOJ 3529: [Sdoi2014]数表 [莫比乌斯反演 树状数组]

    3529: [Sdoi2014]数表 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 1399  Solved: 694[Submit][Status] ...

  6. BZOJ 3529 [Sdoi2014]数表 ——莫比乌斯反演 树状数组

    $ans=\sum_{i=1}^n\sum_{j=1}^n\sigma(gcd(i,j))$ 枚举gcd为d的所有数得到 $ans=\sum_{d<=n}\sigma(d)*g(d)$ $g(d ...

  7. luogu3312 [SDOI2014]数表 (莫比乌斯反演+树状数组)

    link \(\sum_{i=1}^n\sum_{j=1}^m[s(\gcd(i,j))\le a]s(\gcd(i,j))\) \(=\sum_{p=1}^ns(p)[s(p)\le a]\sum_ ...

  8. bzoj 3529 数表 莫比乌斯反演+树状数组

    题目大意: 有一张N×m的数表,其第i行第j列(1 < =i < =礼,1 < =j < =m)的数值为能同时整除i和j的所有自然数之和.给定a,计算数表中不大于a的数之和. ...

  9. 【BZOJ3529】【莫比乌斯反演 + 树状数组】[Sdoi2014]数表

    Description 有一张N×m的数表,其第i行第j列(1 < =i < =礼,1 < =j < =m)的数值为 能同时整除i和j的所有自然数之和.给定a,计算数表中不大于 ...

  10. 【BZOJ3529】【SDOI2014】数表 (莫比乌斯反演+树状数组)

    传送门 Description 有一张$n\times m$的数表,其第$i$行第$j$列 $(1≤i≤n,1≤j≤m)$ 的数值为能同时整除$i$和$j$的所有自然数之和.现在给定$a$,计算数表中 ...

随机推荐

  1. 网络配置:linux学习第一篇

    1.      先使用dhclient获取ip 再使用命令ip addr查看获取到的ip 2.      设置静态IP 编辑网卡配置文件,路径: 3.      重启网络服务 命令:systemctl ...

  2. 【Tools】Pycharm2017 windows安装与修改中文界面教程

    [windows] 1.到官网下载Pycharm最新版 https://www.jetbrains.com/pycharm/download/#section=windows 2.安装激活 Pycha ...

  3. 【Unity3D技术文档翻译】第1.8篇 AssetBundles 问题及解决方法

    上一章:[Unity3D技术文档翻译]第1.7篇 AssetBundles 补丁更新 本章原文所在章节:[Unity Manual]→[Working in Unity]→[Advanced Deve ...

  4. 【项目记录】-路灯监测 gmap.net

    需求 2016年5月,客户要求在地图上显示路灯及数据,分析数据生成报表,以便查看分析路灯情况. 选型 国外项目就不考虑国内的地图了,开始想使用google的web地图,考虑到地图上标记物过多影响性能及 ...

  5. python进阶学习笔记(二)

    1.模块和包的概念 python的解决方案是把同名的模块放到不同的包中 1.1,导入模块 要使用一个模块,我们必须首先导入该模块.Python使用import语句导入一个模块.例如,导入系统自带的模块 ...

  6. java-redis字符类数据操作示例(一)

    对于大部分程序猿来讲,学习新知识重在编码实践,于我也是这样.现在初识redis,一直看文章难免感觉是浮光掠影,印象不深.所以间隙中,将自己的测试代码整理成博客,旨在加深记忆并提醒自己对待编程要用心沉下 ...

  7. dubbox系列【二】——dubbox admin环境搭建

    1.环境 JDK7+ zookeeper 3.3.6 tomcat 7+ 2.搭建步骤 (1)安装zookeeper 单机模式:参考这里. 伪集群模式:参考这里. ①单机或伪集群模式均可. ②笔者使用 ...

  8. IntentService源码

    原文地址IntentService源码分析 @Override public void onCreate() { super.onCreate(); HandlerThread thread = ne ...

  9. Oracle总结【视图、索引、事务、用户权限、批量操作】

    前言 在Oracle总结的第一篇中,我们已经总结了一些常用的SQL相关的知识点了...那么本篇主要总结关于Oralce视图.序列.事务的一些内容... 在数据库中,我们可以把各种的SQL语句分为四大类 ...

  10. css入门二-常用样式

    css入门二-常用样式总结 基本标签样式 背景色background-color 高度height; 宽度width; 边框对齐以及详细设定举例 width/*宽度*/: 80%; height/*高 ...