[SDOI2014]数表
题目描述
有一张N*m的数表,其第i行第j列(1 < =i < =礼,1 < =j < =m)的数值为能同时整除i和j的所有自然数之和。给定a,计算数表中不大于a的数之和。
输入输出格式
输入格式:
输入包含多组数据。 输入的第一行一个整数Q表示测试点内的数据组数,接下来Q行,每行三个整数n,m,a(|a| < =10^9)描述一组数据。
输出格式:
对每组数据,输出一行一个整数,表示答案模2^31的值。
输入输出样例
2
4 4 3
10 10 5
20
148
说明
1 < =N.m < =10^5 , 1 < =Q < =2*10^4
题解:
令F[i]为i的因数和,g[i]为gcd=i的数量
则ans=∑∑F[gcd(i,j)]
=∑dF[d]g[d]
g[d]=∑i∑j[gcd(i,j)==d] (i<=n,j<=m)
=∑∑[gcd(i,j)==1] (i<=n/d,j<=m/d)
=∑∑∑pμ(p) (p|i,p|j)
=∑pμ(p)∑∑1 (p<=n/d,i<=n/dp,j<=m/dp)
=∑pμ(p)[n/dp][m/dp]
ans=∑dF[d]∑pμ(p)[n/dp][m/dp]
令k=dp
g[d]=∑k[n/k][m/k]μ(k/d) (k<=n,d|k)
ans=∑k[n/k][m/k]∑dF[d]μ(k/d) (d|k)
令f[k]=∑dF[d]μ(k/d) (d|k)
处理处f的前缀和就行了,对于每个F[i],把它的f[x*i]全部加上F[i]*mu[x]
至于小于a的条件,离线按a排序,将F按从小到大顺序加入维护
分块就不说了
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
struct data
{
int id,n,m,a;
}a[];
int Mod=(<<);
int f[],mu[],ms[],prime[],tot;
int c[],p[],now=,ans[];
bool vis[];
bool cmp2(data a,data b)
{
return a.a<b.a;
}
bool cmp(int a,int b)
{
return f[a]<f[b];
}
void mobius()
{int i,j;
mu[]=;
ms[]=;f[]=;
for (i=;i<=;i++)
{
if (vis[i]==)
{
mu[i]=-;
f[i]=i+;
ms[i]=i;
prime[++tot]=i;
}
for (j=;j<=tot,prime[j]*i<=;j++)
{
vis[i*prime[j]]=;
if (i%prime[j])
{
ms[i*prime[j]]=prime[j];
f[prime[j]*i]=f[i]*f[prime[j]];
mu[i*prime[j]]=-mu[i];
}
else
{
ms[prime[j]*i]=ms[i]*prime[j];
if (ms[i]==i) f[prime[j]*i]=(ms[prime[j]*i]*prime[j]-)/(prime[j]-);
else f[prime[j]*i]=f[i/ms[i]]*f[prime[j]*ms[i]];
mu[i*prime[j]]=;
break;
}
}
}
}
void add(int x,int d)
{
while (x<=)
{
c[x]+=d;
x+=(x&(-x));
}
}
int query(int x)
{
int s=;
while (x)
{
s+=c[x];
x-=(x&(-x));
}
return s;
}
void solve(int x)
{int j,pos,i,lasts,ss;
for (;f[p[now]]<=a[x].a;now++)
{
for (j=;p[now]*j<=;j++)
add(p[now]*j,f[p[now]]*mu[j]);
}
pos=;
lasts=;
for (i=;i<=a[x].n;i=pos+)
{
pos=min(a[x].n/(a[x].n/i),a[x].m/(a[x].m/i));
ss=query(pos);
ans[a[x].id]=(ans[a[x].id]+(a[x].n/i)*(a[x].m/i)*(ss-lasts));
//cout<<ans[a[x].id]<<' ';
lasts=ss;
}
while (ans[a[x].id]<) ans[a[x].id]+=Mod;
//cout<<endl;
}
int main()
{int T,i,j;
cin>>T;
for (i=;i<=T;i++)
{
scanf("%d%d%d",&a[i].n,&a[i].m,&a[i].a);
if (a[i].n>a[i].m) swap(a[i].n,a[i].m);
a[i].id=i;
}
mobius();
for (i=;i<=;i++) p[i]=i;
sort(p+,p+,cmp);
sort(a+,a+T+,cmp2);
for (i=;i<=T;i++)
solve(i);
for (i=;i<=T;i++)
printf("%d\n",ans[i]);
}
[SDOI2014]数表的更多相关文章
- BZOJ 3529: [Sdoi2014]数表 [莫比乌斯反演 树状数组]
3529: [Sdoi2014]数表 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 1399 Solved: 694[Submit][Status] ...
- 【BZOJ 3529】 [Sdoi2014]数表 (莫比乌斯+分块+离线+树状数组)
3529: [Sdoi2014]数表 Description 有一张N×m的数表,其第i行第j列(1 < =i < =礼,1 < =j < =m)的数值为能同时整除i和j的所有 ...
- BZOJ 3259 [Sdoi2014]数表 (莫比乌斯反演 + 树状数组)
3529: [Sdoi2014]数表 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 2321 Solved: 1187[Submit][Status ...
- 洛咕3312 [SDOI2014]数表
洛咕3312 [SDOI2014]数表 终于独立写出一道题了...真tm开心(还是先写完题解在写的) 先无视a的限制,设\(f[i]\)表示i的约数之和 不妨设\(n<m\) \(Ans=\su ...
- 【BZOJ3529】[Sdoi2014]数表 莫比乌斯反演+树状数组
[BZOJ3529][Sdoi2014]数表 Description 有一张N×m的数表,其第i行第j列(1 < =i < =礼,1 < =j < =m)的数值为能同时整除i和 ...
- 洛谷 P3312 [SDOI2014]数表 解题报告
P3312 [SDOI2014]数表 题目描述 有一张\(N*M\)的数表,其第\(i\)行第\(j\)列(\(1\le i \le n\),\(1 \le j \le m\))的数值为能同时整除\( ...
- BZOJ3529 [Sdoi2014]数表 【莫比乌斯反演】
3529: [Sdoi2014]数表 Time Limit: 10 Sec Memory Limit: 512 MB Submit: 2151 Solved: 1080 [Submit][Status ...
- BZOJ[Sdoi2014]数表 莫比乌斯反演
[Sdoi2014]数表 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 2383 Solved: 1229[Submit][Status][Disc ...
- 洛谷P3312 - [SDOI2014]数表
Portal Solution 共\(T(T\leq2\times10^4)\)组测试数据.给出\(n,m(n,m\leq10^5),a(a\leq10^9)\),求\[ \sum_{i=1}^n\s ...
- bzoj [SDOI2014]数表 莫比乌斯反演 BIT
bzoj [SDOI2014]数表 莫比乌斯反演 BIT 链接 bzoj luogu loj 思路 \[ \sum\limits_{i=1}^{n}\sum\limits_{j=1}^{m}a*[f[ ...
随机推荐
- 高级软件工程2017第5次作业—— 团队项目:需求改进&系统设计
Deadline:2017-10-23(周一) 21:00pm 注:以下内容参考 集大作业 1.评分规则: 按时交 - 有分,检查的项目包括后文的四个方面 需求&原型改进 - 20分 系统设计 ...
- 2017-2018-1 Java演绎法 小组成员贡献量汇总
[第一周]贡献量(31) [说明] 完成情况 是指 每次是否全部完成分配的任务,如果全部完成贡献量记为1,否则记为0,与贡献量(时间量)相加计算贡献比例,由于前十周有具体的任务分配,Alpha阶段(第 ...
- Node入门教程(4)第三章:第一个 Nodejs 程序
第一个 Nodejs 程序 本教程仅适合您已经有一定的JS编程的基础或者是后端语言开发的基础.如果您是零基础,建议您先学一下老马的前端免费视频教程 第一步:创建项目文件夹 首先创建 demos 文件夹 ...
- centos 6.5安装并配置mysql
折腾了半天终于把mysql安装并配置好了,以下是安装步骤和遇到问题的解决方式 1.查看机器上是否已经安装了mysql或其相关项 # yum list installed | grep mysql如果安 ...
- 04_Linux目录文件操作命令1(mv ls cd...)_我的Linux之路
上一节已经给大家讲了Linux的目录结构,相信大家已经对Linux的整个目录结构有所了解 现实中,服务器(包含Linux,Unix,windows server)一般都摆放在机房里,因为一个机房摆放了 ...
- python 基础 set 集合类型补充
为啥今天又重提这个数据类型呢?平时用的少,等要用起来的时候才发现,自己对这块啥都不知道了,so,今天就把这块再梳理一下咯. 一.set集合,是一个无序且不重复的元素集合.这一点是非常重要的. 二.集合 ...
- wordpress | WP Mail SMTP使用QQ邮箱发布失败的解决办法
在使用contact form 7插件时遇到邮件发送失败的问题,经过检查发现是因为服务器不支持mail()函数,判断是否支持mail()函数可以参考http://www.diyzhan.com/201 ...
- 使用cxf创建webservice 出现timeOut的问题,设置spring超时时间
<?xml version="1.0" encoding="UTF-8"?> <beans xmlns="http://www.sp ...
- linux下的Shell编程(8)自定义函数
Shell Script中也可以使用自定义的函数,其语法形式如下: functionname() { - }
- SpringCloud的部署模型
http://www.th7.cn/Program/java/201608/919853.shtml