[SDOI2014]数表

• 题目大意：有t组数据每组数据给出3个数，n,m,a,求
• 
  
  
  

  其中（n,m<=10^5,a<=10^9,t<=2*10^4）

• 解题思路：首先，我们忽略a的限制，也就是想办法求

• 
  

• 现在，我们可以换一下枚举顺序，先枚举gcd，设在n，m固定的情况下，gcd为i的格点的个数为g（i）原式就变成了

• 其中
  

• 
  

• 这时上一波反演就变成了

• 所以我们现在要求的经过转化变成了

• 
  

• 交换一下枚举顺序

• 设函数，那么，我们要求的就是

• 发现和的值都是最多只有级别个，所以我们可以求出f函数的前缀和，然后对于每一种不同的值单独计算在求和，再回到题目中a的限制，当有a的限制的时候可以离线按询问的a值排序，将函数的值也排序，然后每次增加a值的时候就对f函数的前缀和做贡献，可以注意到，若新增贡献的值为那么会受到影响的f函数值，j都i的倍数或i本身，对于一个1-n的约数数量和是nlogn级别的
  

• 具体证明：考虑1是所有n个数的约数，2是个数的约数，3是个数的约数，对于数字i，i是个数的约数，那么

• 1-n的约数数量和相当于求然而
  

• 

  然后当式子大于n时需要前1项就可以了，当式子大于2n时需要前3项，式子大于3n时要前7项，式子大于4n时要前15项……式子大于kn时要前-1项，而整个式子一共才n项，所以式子的值不会超过nlogn。

  我们不能暴力每次O（n）修改前缀和，应该用树状数组维护，这样这道题就解决了。主要的时间消耗在，求值和修改前缀和上，求值的时间复杂度为修改前缀和的时间复杂度是，时间上勉强卡过。
•                                                                       推荐番：《你的名字。》

最后：上代码

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
struct data
{long long a,n,m,i;};
data a[20005];
struct dat
{long long v,i;};
dat F[100005];
long long tre[100005],tiao,u[100005],zhi[100005],pp,line[100005];
long long h[1000],t[1000],hh[1000],tt[1000],jie[1000],jiee[1000],q,qq,dd,dang,shuchu[20005];
long long h1[1000],t1[1000],hh1[1000],tt1[1000],jie1[1000],jiee1[1000],gg,ans;
void jia(long long wei,long long k)
{
	while (wei<=100000)
	{
		tre[wei]+=k;
		wei+=wei & (-wei);
    }
}
long long qiu (long long l,long long r)
{
	long long p1=0,p2=0;
	l--;
	while (l>0)
	{
		p1+=tre[l];
		l-=l & (-l);
	}
	while (r>0)
	{
		p2+=tre[r];
		r-=r & (-r);
	}
	return p2-p1;
}
bool cmp(data a,data b)
{return a.a<b.a;}
bool tmp(dat a,dat b)
{return a.v<b.v;}
void zuo(long long n,long long m,long long weiwei)
{
	q=0;qq=0;
	for (long long i=1;i*i<=n;i++)	
   {
   	qq++;
   	 jie[i]=i;h[i]=n/(i+1)+1;t[i]=n/i;
   }
   for (long long i=1;i*i<n;i++)
   {
   	if (n/i!=jie[qq]){q++;jiee[i]=n/i;hh[i]=i;tt[i]=i;}
   }
   for (long long i=q;i>=1;i--) {qq++;jie[qq]=jiee[i];h[qq]=hh[i];t[qq]=tt[i];}
   dd=qq;
   q=0;qq=0;
   for (long long i=1;i*i<=m;i++)	
   {
   	qq++;
   	 jie1[i]=i;h1[i]=m/(i+1)+1;t1[i]=m/i;
   }
   for (long long i=1;i*i<m;i++)
   {
   	if (m/i!=jie1[qq]){q++;jiee1[i]=m/i;hh1[i]=i;tt1[i]=i;}
   }
   for (long long i=q;i>=1;i--) {qq++;jie1[qq]=jiee1[i];h1[qq]=hh1[i];t1[qq]=tt1[i];}
   gg=0;ans=0;q=dd;
   while (gg<min(n,m))
   {
   	 if (t[q]<t1[qq]){ans+=qiu(gg+1,t[q])*jie[q]*jie1[qq];gg=t[q];q--;}
   	 else {ans+=qiu(gg+1,t1[qq])*jie[q]*jie1[qq];gg=t1[qq];qq--;}
   }
   shuchu[weiwei]=ans%(1 << 31);
}
int main()
{
   cin>>tiao;
   for (long long i=1;i<=tiao;i++)
   {
   	scanf("%lld %lld %lld",&a[i].n,&a[i].m,&a[i].a);
   	a[i].i=i;
   }
   sort(a+1,a+1+tiao,cmp);
   for (long long i=1;i<=100000;i++)
   {
   	for (long long j=1;j<=100000/i;j++)F[i*j].v+=i;
   	F[i].i=i;
   }
   sort(F+1,F+100001,tmp);
   u[1]=1;u[2]=-1;
   for (long long i=2;i<=100000;i++)
   {
   	if (zhi[i]==0){pp++;line[pp]=i;u[i]=-1;}
   	  for (long long j=1;j<=pp;j++)
   	  {
   	  	 if (i*line[j]>100000)break;
   	     if (i%line[j]==0){zhi[i*line[j]]=1;break;}
			zhi[i*line[j]]=1;u[i*line[j]]=u[i]*-1;
	  }
   }
   dang=0;
    for (long long i=1;i<=tiao;i++)
    {
    	while (dang<100000 && F[dang+1].v<=a[i].a)
    	{
    		for (long long j=1;j<=100000/F[dang+1].i;j++)jia(F[dang+1].i*j,u[j]*F[dang+1].v);
    		dang++;
		}
		zuo(a[i].n,a[i].m,a[i].i);
	}
	for (int i=1;i<=tiao;i++)printf("%lld\n",shuchu[i]);
}