题目描述

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

输入输出格式

输入格式:

输入包含多组数据。 输入的第一行一个整数Q表示测试点内的数据组数,接下来Q行,每行三个整数n,m,a(|a| < =10^9)描述一组数据。

输出格式:

对每组数据,输出一行一个整数,表示答案模2^31的值。

输入输出样例

输入样例#1:

2
4 4 3
10 10 5
输出样例#1:

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]=∑ij[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]数表的更多相关文章

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

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

  2. 【BZOJ 3529】 [Sdoi2014]数表 (莫比乌斯+分块+离线+树状数组)

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

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

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

  4. 洛咕3312 [SDOI2014]数表

    洛咕3312 [SDOI2014]数表 终于独立写出一道题了...真tm开心(还是先写完题解在写的) 先无视a的限制,设\(f[i]\)表示i的约数之和 不妨设\(n<m\) \(Ans=\su ...

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

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

  6. 洛谷 P3312 [SDOI2014]数表 解题报告

    P3312 [SDOI2014]数表 题目描述 有一张\(N*M\)的数表,其第\(i\)行第\(j\)列(\(1\le i \le n\),\(1 \le j \le m\))的数值为能同时整除\( ...

  7. BZOJ3529 [Sdoi2014]数表 【莫比乌斯反演】

    3529: [Sdoi2014]数表 Time Limit: 10 Sec Memory Limit: 512 MB Submit: 2151 Solved: 1080 [Submit][Status ...

  8. BZOJ[Sdoi2014]数表 莫比乌斯反演

    [Sdoi2014]数表 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 2383  Solved: 1229[Submit][Status][Disc ...

  9. 洛谷P3312 - [SDOI2014]数表

    Portal Solution 共\(T(T\leq2\times10^4)\)组测试数据.给出\(n,m(n,m\leq10^5),a(a\leq10^9)\),求\[ \sum_{i=1}^n\s ...

  10. bzoj [SDOI2014]数表 莫比乌斯反演 BIT

    bzoj [SDOI2014]数表 莫比乌斯反演 BIT 链接 bzoj luogu loj 思路 \[ \sum\limits_{i=1}^{n}\sum\limits_{j=1}^{m}a*[f[ ...

随机推荐

  1. Alpha第八天

    Alpha第八天 听说 031502543 周龙荣(队长) 031502615 李家鹏 031502632 伍晨薇 031502637 张柽 031502639 郑秦 1.前言 任务分配是VV.ZQ. ...

  2. c语言第1次作业

    一.PTA实验作业 题目1:7-3 温度转换 本题要求编写程序,计算华氏温度150°F对应的摄氏温度.计算公式:C=5×(F−32)/9,式中:C表示摄氏温度,F表示华氏温度,输出数据要求为整型. 1 ...

  3. 本地通知UILocalNotification

    1.增加一个本地推送 //设置20秒之后  ]; //chuagjian一个本地推送 UILocalNotification *noti = [[[UILocalNotification alloc] ...

  4. HTTP协议以及HTTP2.0/1.1/1.0区别

    HTTP协议以及HTTP2.0/1.1/1.0区别 一.简介 摘自百度百科: 超文本传输协议(HTTP,HyperText Transfer Protocol)是互联网上应用最为广泛的一种网络协议.所 ...

  5. Python内置函数(62)——exec

    英文文档: exec(object[, globals[, locals]]) This function supports dynamic execution of Python code. obj ...

  6. emqtt 试用(六)系统主题

    $SYS-系统主题 EMQ 消息服务器周期性发布自身运行状态.MQTT 协议统计.客户端上下线状态到 $SYS/ 开头系统主题. $SYS 主题路径以 "$SYS/brokers/{node ...

  7. Scala入门(1)Linux下Scala(2.12.1)安装

    Scala入门(1)Linux下Scala(2.12.1)安装 一.文件准备 1.1 文件名称 scala-2.12.1.tgz 1.2 下载地址 http://www.scala-lang.org/ ...

  8. Jetty入门(1-1)Jetty入门教程

    一.Jetty是什么? 1.Jetty 是一个Java语言编写的,开源的Servlet容器和应用服务器. Jetty 极度轻量级.高便携性.功能强大.灵活和扩展性好,而且支持各种技术如SPDY.Web ...

  9. 阿里云API网关(1)服务网关的产品概述

    网关指南: https://help.aliyun.com/document_detail/29487.html?spm=5176.doc48835.6.550.23Oqbl 网关控制台: https ...

  10. Python学习之dict和set

    #coding=utf-8 # dict dict= {'bob': 40, 'andy': 30} print dict['bob'] # 通过dict提供的get方法,如果key不存在,可以返回N ...