数表

数表

时间限制: 1 Sec 内存限制: 512 MB
题目描述
有一张n×m的数表，其第i行第j列（1≤i≤n，1≤j≤m）的数值为能同时整除i和j的所有自然数之和。给定a，计算数表中不大于a的数之和。
输入
输入包含多组数据。
输入的第一行一个整数Q表示测试点内的数据组数，接下来Q行，每行三个整数n，m，a(|a|≤109)描述一组数据。

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

样例输入
2
4 4 3
10 10 5

样例输出
20
148

来源
sdoi2014R1D1

题解

第一次写数论题题解，写的比较差表怪我。
题目要求一个式子：sigma(f(gcd(i,j)))，i<=n,j<=m,f(gcd(i,j))<=a；
显然有：g(i)=miu(d/i)*[n/d][m/d]，其中，d|i，[]为下取整，g(i)表示n,m内gcd=i的数对的组数。
所以，ans=sigma(f(i)*g(i))，f(i)表示i的因数和，且f(i)<=a。
进而，ans=sigma([n/d][m/d])sigma(f[i]*miu(d/i))(i|d)。
由于本题有a的限制，并且有多组询问。所以考虑将操作离线，按照a的大小从小到大跑。
用一个树状数组来维护f[i]*miu[d/i]的和。
对于公式前半部分，通过枚举商的经典做法操作即可。

代码

#include<iostream> 
#include<cstdio> 
#include<cstdlib> 
#include<cstring> 
#include<string> 
#include<cmath> 
#include<algorithm> 
#define N 100010
using namespace std;
int Q,maxn,cnt,ans[N],prime[N],u[N],check[N],c[N];
struct info{
  int x,val;
  bool operator<(const info &p)
  const{return val<p.val;}
}f[N];
struct node{
  int n,m,a,pos;
  bool operator<(const node &p)
  const{return a<p.a;}
}t[N];

class bit
{
  public:
  void modify(int x,int val)
  {
    for(;x<=maxn;x+=x&-x)c[x]+=val;
  }
  int qry(int x)
  {
    int res=0;
    for(;x;x-=x&-x)res+=c[x];
    return res;
  }
}T;

int main()
{
  int n,m,a,pos;
  scanf("%d",&Q);
  for(int i=1;i<=Q;i++)
  {
    scanf("%d%d%d",&n,&m,&a);
    if(n>m)swap(n,m);maxn=max(maxn,n);
    t[i]=(node){n,m,a,i};
  }
  u[1]=1;
  for(int i=2;i<=maxn;i++)
  {
    if(!check[i])prime[++cnt]=i,u[i]=-1;
    for(int j=1;j<=cnt&&i*prime[j]<=maxn;j++)
    {
      check[i*prime[j]]=1;
      if(!(i%prime[j])){u[i*prime[j]]=0;break;}
      u[i*prime[j]]=-u[i];
    }
  }
  for(int i=1;i<=maxn;f[i].x=i,i++)
    for(int j=i;j<=maxn;j+=i)f[j].val+=i;
  sort(f+1,f+maxn+1);sort(t+1,t+Q+1);
  for(int i=1,st=1;i<=Q;i++)
  {
    n=t[i].n;m=t[i].m;a=t[i].a;pos=t[i].pos;
    for(;f[st].val<=a&&st<=maxn;st++)
      for(int j=f[st].x;j<=maxn;j+=f[st].x)
        T.modify(j,u[j/f[st].x]*f[st].val);
    for(int l=1,r;l<=n;l=r+1)
    {
      r=min(n/(n/l),m/(m/l));
      ans[pos]+=(n/l)*(m/l)*(T.qry(r)-T.qry(l-1));
    }
  }
  for(int i=1;i<=Q;i++)printf("%d\n",ans[i]&0x7fffffff);
  return 0;
}