【SDOI2014】【BZOJ3529】数表
Description
有一张N×m的数表,其第i行第j列(1 < =i < =礼。1 < =j < =m)的数值为
能同一时候整除i和j的全部自然数之和。给定a,计算数表中不大于a的数之和。
Input
输入包括多组数据。
输入的第一行一个整数Q表示測试点内的数据组数,接下来Q行,每行三个整数n,m,a(|a| < =10^9)描写叙述一组数据。
Output
对每组数据,输出一行一个整数。表示答案模2^31的值。
Sample Input
2
4 4 3
10 10 5
Sample Output
20
148
HINT
1 < =N.m < =10^5 。 1 < =Q < =2×10^4
Source
Round 1 Day 1
在反演上的处理和DzyLovesMath1是相似的属于同一类题目
问题在于询问有10^4个所以要BIT维护一下前缀和
推公式过程懒得放了…贴个Po姐课件的图吧
(从未见过如此丧病的反演..反演还要加数据结构T_T
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 100010
#define lowbit(x) (x&(-x))
#define LL long long
#define MAXINT (0x7fffffff)
#define GET (ch>='0'&&ch<='9')
using namespace std;
int T,maxn;
int c[MAXN];
bool not_prime[MAXN];
int prime[MAXN],mu[MAXN]={0,1},top;
int ans[MAXN];
struct Query
{
int n,m,a,id;
bool operator <(const Query& t)const {return a<t.a;}
}q[MAXN];
struct num
{
int a,b;
bool operator <(const num& x)const {return a==x.a?b<x.b:a<x.a;}
}f[MAXN];
void init()
{
for (int i=2;i<=maxn;i++)
{
if (!not_prime[i]) prime[++top]=i,mu[i]=-1;
for (int j=1;j<=top&&i*prime[j]<=maxn;j++)
{
not_prime[i*prime[j]]=1;mu[i*prime[j]]=-mu[i];
if (i%prime[j]==0) {mu[i*prime[j]]=0;break;}
}
}
for (int i=1;i<=maxn;i++)
{
for (int j=i;j<=maxn;j+=i) f[j].a+=i;
f[i].b=i;
}
}
void add(int x,int delta) {for (int i=x;i<=maxn;i+=lowbit(i)) c[i]+=delta;}
int query(int x)
{
int ret=0;
for (int i=x;i;i-=lowbit(i)) ret+=c[i];
return ret;
}
void in(int &x)
{
char ch=getchar();x=0;int flag=1;
while (!GET) flag=(ch=='-')?
-1:1,ch=getchar();
while (GET) x=x*10+ch-'0',ch=getchar();x*=flag;
}
void calc(int x)
{
int id=q[x].id,n=q[x].n,m=q[x].m,last=0,t=min(n,m);
for (int i=1;i<=t;i=last+1)
{
last=min(n/(n/i),m/(m/i));
ans[id]+=(n/i)*(m/i)*(query(last)-query(i-1));
}
}
int main()
{
in(T);int now=0;
for (int i=1;i<=T;i++) in(q[i].n),in(q[i].m),in(q[i].a),q[i].id=i,maxn=max(maxn,max(q[i].n,q[i].m));
init();sort(q+1,q+T+1);sort(f+1,f+maxn+1);
for (int i=1;i<=T;i++)
{
while (now+1<=maxn&&f[now+1].a<=q[i].a)
{
now++;
for (int j=f[now].b;j<=maxn;j+=f[now].b) add(j,f[now].a*mu[j/f[now].b]);
}
calc(i);
}
for (int i=1;i<=T;i++) printf("%d\n",ans[i]&MAXINT);
}
【SDOI2014】【BZOJ3529】数表的更多相关文章
- 【BZOJ3529】【SDOI2014】 数表
Time Limit: 10 Sec Memory Limit: 512 MB Description 有一张\(n×m\)的数表,其第i行第j列(\(,1 \le i \leq n,1 \le ...
- 【BZOJ3529】【SDOI2014】数表
Time Limit: 1000 ms Memory Limit: 512 MB Description 有一张n×m的数表,其第i行第j列 (1≤i≤n,1≤j≤m)的数值为能同时整除i和j的所有自 ...
- 【BZOJ3529】【SDOI2014】数表 (莫比乌斯反演+树状数组)
传送门 Description 有一张$n\times m$的数表,其第$i$行第$j$列 $(1≤i≤n,1≤j≤m)$ 的数值为能同时整除$i$和$j$的所有自然数之和.现在给定$a$,计算数表中 ...
- BZOJ 3529 【SDOI2014】 数表
题目链接:数表 我们一起来膜PoPoQQQ大爷的题解吧Orz 首先我们来考虑没有\(a\)的限制该怎么做.显然交换\(n\),\(m\)答案不变,所以后面默认\(n \le m\). 我们定义两个函数 ...
- 【SDOI2014】数表
题面 题解 不管$a$的限制 我们要求的东西是:($\sigma(x)$是$x$的约数个数和) $ \sum_{i=1}^n\sum_{j=1}^m\sigma(gcd(i,j)) $ 设$f(x)= ...
- [BZOJ3529]数表
假设$n\leq m$,我们先不考虑$\leq a$的限制 $\sum\limits_{i=1}^n\sum\limits_{j=1}^m\sigma((i,j))=\sum\limits_{T=1} ...
- 【JZOJ3623】【SDOI2014】数表(table) 树状数组+离线+莫比乌斯反演
题面 100 \[ Ans=\sum_{i=1}^n\sum_{j=1}^mg(gcd(i,j)) \] 其中, \[ g(d)=\sum_{i|d}i \] 我们注意到\(gcd(i,j)\)最多有 ...
- 【BZOJ3529】[Sdoi2014]数表 莫比乌斯反演+树状数组
[BZOJ3529][Sdoi2014]数表 Description 有一张N×m的数表,其第i行第j列(1 < =i < =礼,1 < =j < =m)的数值为能同时整除i和 ...
- [bzoj3529][Sdoi2014]数表_树状数组_莫比乌斯反演
数表 bzoj-3529 Sdoi-2014 题目大意:n*m的数表,第i行第j列的数是同时整除i和j的所有自然数之和.给定a,求数表中所有不超过a的和. 注释:$1\le n,m \le 10^5$ ...
- BZOJ3529 [Sdoi2014]数表 【莫比乌斯反演】
3529: [Sdoi2014]数表 Time Limit: 10 Sec Memory Limit: 512 MB Submit: 2151 Solved: 1080 [Submit][Status ...
随机推荐
- 性能测试工具——Mxdperfstat
Mxdperfstat是一款mxd性能检测工具,使用它来测试专题地图的性能非常不错! 获取工具 https://www.arcgis.com/home/item.html?id=a269d03aa1c ...
- PHP Warning: DOMDocument::loadHTML(): htmlParseEntityRef: expecting ';' in Entity,
$html = file_get_contents("http://www.somesite.com/"); $dom = new DOMDocument(); $dom-> ...
- ZH奶酪:PHP (爬虫)下载图片
原文地址:http://www.phpfensi.com/php/20140107/1128.html 通过图片地地址把图片保存到本址,这里我们直接通过readfile读取然后通过fopen保存即可, ...
- CSS 之 样式优先级机制
多重样式(Multiple Styles):如果外部样式.内部样式和内联样式同时应用于同一个元素,就是使多重样式的情况. 1.一般情况下,优先级从低到高,如下: 浏览器缺省(即浏览器默认) →(外部样 ...
- Ganglia监控Hadoop集群的安装部署
一. 安装环境 Ubuntu server 12.04 安装gmetad的机器:192.168.52.105 安装gmond的机器:192.168.52.31,192.168.52.32,192.16 ...
- Python取得系统进程列表
一.上代码 import psutil for proc in psutil.process_iter(): try: pinfo = proc.as_dict(attrs=['pid', 'name ...
- python 大小端数据转换
"6fe28c0ab6f1b372c1a6a246ae63f74f931e8365e15a089c68d6190000000000".decode('hex')[::-1].enc ...
- Influxdb数据存储
环境: CentOS6.5_x64 InfluxDB版本:1.1.0 InfluxDB存储引擎看起来很像一个LSM Tree,它包含预写日志和类似存储在LSM Tree中的SSTables只读数据. ...
- CentOS yum时出现"Could not retrieve mirrorlist"
问题描述: CentOS 6.x minimal(最小化) 安装, CentOS yum install net-tools 时出现"Could not retrieve mirrorlis ...
- LeetCode 225 Implement Stack using Queues(用队列来实现栈)(*)
翻译 用队列来实现栈的例如以下操作. push(x) -- 将元素x加入进栈 pop() -- 从栈顶移除元素 top() -- 返回栈顶元素 empty() -- 返回栈是否为空 注意: 你必须使用 ...