数表(bzoj 3529)

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

/*
    在PoPoQQQ的ppt上看的,感觉各种化简很神。
*/
#include<cstdio>
#include<iostream>
#include<algorithm>
#define N 1000010
using namespace std;
int prime[N],num,mul[N],t[N],mark[N],ans[N],mx;
struct node{int n,m,a,id;}q[N];
pair<int,int> F[N];
bool operator<(node a,node b){
    return a.a<b.a;
}
void add(int x,int val){
    for(int i=x;i<=mx;i+=i&-i) t[i]+=val;
}
int query(int x){
    int tmp=0;
    for(int i=x;i;i-=i&-i) tmp+=t[i];
    return tmp;
}
void get_prime(){
    mul[1]=1;
    for(int i=2;i<=mx;i++){
        if(!mark[i]) prime[++num]=i,mul[i]=-1;
        for(int j=1;j<=num&&i*prime[j]<=mx;j++){
            mark[i*prime[j]]=1;
            mul[i*prime[j]]=-mul[i];
            if(i%prime[j]==0){mul[i*prime[j]]=0;break;}
        }
    }
    for(int i=1;i<=mx;i++)
        for(int j=i;j<=mx;j+=i)
            F[j].first+=i;
    for(int i=1;i<=mx;i++) F[i].second=i;
}
void solve(int x){
    int id=q[x].id,n=q[x].n,m=q[x].m;
    for(int i=1,j;i<=q[x].n;i=j+1){
        j=min(n/(n/i),m/(m/i));
        ans[id]+=(n/i)*(m/i)*(query(j)-query(i-1));
    }
}
int main(){
    int Q;scanf("%d",&Q);
    for(int i=1;i<=Q;i++){
        scanf("%d%d%d",&q[i].n,&q[i].m,&q[i].a);
        if(q[i].n>q[i].m) swap(q[i].n,q[i].m);
        mx=max(mx,q[i].n);q[i].id=i;
    }
    get_prime();
    sort(q+1,q+Q+1);
    sort(F+1,F+mx+1);
    int now=0;
    for(int i=1;i<=Q;i++){
        while(now+1<=mx&&F[now+1].first<=q[i].a){
            now++;
            for(int j=F[now].second;j<=mx;j+=F[now].second)
                add(j,F[now].first*mul[j/F[now].second]);
        }
        solve(i);
    }
    for(int i=1;i<=Q;i++)
        printf("%d\n",ans[i]&0x7fffffff);
    return 0;
}
时间: 2024-12-28 21:14:44

数表(bzoj 3529)的相关文章

BZOJ 3529: [Sdoi2014]数表

Description 有一张N×m的数表,其第i行第j列(1 < =i < =n,1 < =j < =m)的数值为能同时整除i和j的所有自然数之和.给定a,计算数表中不大于a的数之和. HINT 1 < =N.m < =10^5  , 1 < =Q < =2×10^4 Solution 先假设a的限制不存在 正面把答案强行写出来是这样的 其中F(i)为i的约数和,可以线性筛处理 然后慢慢化简其中的各个部分 令g(i)为1<=x<=n,1<

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

3529: [Sdoi2014]数表 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 1399  Solved: 694[Submit][Status][Discuss] Description 有一张N×m的数表,其第i行第j列(1 < =i < =礼,1 < =j < =m)的数值为能同时整除i和j的所有自然数之和.给定a,计算数表中不大于a的数之和. Input 输入包含多组数据.    输入的第一行一个整数Q表示测试点内的数据

bzoj 3529 数表

题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=3529 题目大意:令F(i)为i的约数和,多次询问对于1<=x<=n,1<=y<=m,F(gcd(x,y))<=a的所有数对(x,y),求ΣF(gcd(x,y))%(2^31) n,m<=10^5,a<=10^9 首先如果不考虑a的限制 令g(i)为1<=x<=n,1<=y<=m,gcd(x,y)=i的数的个数 那么显然有 利用线性筛

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

题目大意: 有一张N×m的数表,其第i行第j列(1 < =i < =礼,1 < =j < =m)的数值为能同时整除i和j的所有自然数之和.给定a,计算数表中不大于a的数之和. http://wenku.baidu.com/link?url=1zHluup-GXHdByoQXhMRwRu22Uu15y4DztIr1_hKVxjHJmuLQF4_01UQhLEOR7RJIpsGyfD_5fXrx9DE7sY6JeukaNUY83In097GjUOmZ7K ppt课件中讲的很仔细了 1

【bzoj 3529】【sdoi 2014】数表

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 Sa

BZOJ 3529 数表(莫比乌斯反演)

http://www.lydsy.com/JudgeOnline/problem.php?id=3529 思路:令F(i)为i的约数和, 1<=x<=n,1<=y<=m G(i)为i|gcd(x,y)的个数 g(i)为i=gcd(x,y)的个数 G(i)=floor(n/i)*floor(m/i) g(i)=Σu(d/i)*G(d) (i|d) F(i)可以用线性筛O(n)处理 ans=(1<=i<=min(n,m)) Σg(i)*F(i)=Σ F(i)*ΣG(d)*u

BZOJ 3529([Sdoi2014]数表-莫比乌斯反演)

有一张N×m的数表,其第i行第j列(1 < =i < =n,1 < =j < =m)的数值为能同时整除i和j的所有自然数之和.给定a,计算数表中不大于a的数之和模2^31的值模2^31的值. 1<=n.m<=10 5 ,Q<=2×10 4  组询问 记k的约数和f(k) 求∑ n i=1 ∑ m j=1 gcd(i,j)[f(gcd(i,j))≤a] #include<bits/stdc++.h> using namespace std; #defin

●BZOJ 3529 [Sdoi2014]数表

题链: http://www.lydsy.com/JudgeOnline/problem.php?id=3529 题解: 莫比乌斯反演. 按题目的意思,令$f(i)$表示i的所有约数的和,就是要求: $ANS=\sum f(gcd(i,j)),满足1 \leq i \leq n,1 \leq j \leq m,且 f(gcd(i,j))\leq a$ 首先 $f(i)$ 应该还是比较好推的,利用其为积性函数的特点,可以在线性筛时完成计算. 令$g[k]$表示$gcd(i,j)=k$的$(i,j)

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

题目大意:有一张$n*m$的数表,第$i$行第$j$列的数是同时能整除$i,j$的所有数之和,求数表内所有不大于A的数之和 先是看错题了...接着看对题了发现不会做了...刚了大半个下午无果 看了Po姐的题解(Orzzz)才搞懂这道题,搞清楚了莫比乌斯反演的两种经典的卷积形式的不同之处 令$\sigma(i)$表示i的约数个数和 如果去掉A这个限制,则题目是让我们求$\sum_{i=1}^{n}\sum_{j=1}^{m}\sigma(gcd(i,j))$ 考虑如何正确转化式子,让我们能够把不大