BZOJ3309: DZY Loves Math

最新推荐文章于 2020-01-19 19:07:51 发布

Hillan_

最新推荐文章于 2020-01-19 19:07:51 发布

阅读量521

点赞数

CC 4.0 BY-SA版权

分类专栏：数学莫比乌斯反演

本文链接：https://blog.youkuaiyun.com/liutian429073576/article/details/50912446

数学同时被 2 个专栏收录

38 篇文章

订阅专栏

莫比乌斯反演

13 篇文章

订阅专栏

本文介绍了一种用于计算G(x)的高效算法，并提供了详细的题解和代码实现。算法利用筛法和预处理技巧，显著提高了计算速度。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

推G(x)的时候好神啊！！！

题解：http://blog.youkuaiyun.com/popoqqq/article/details/42122413

#include<cstdio>
#include<iostream>
using namespace std;
#define ll long long
char c;
inline void read(int&a)
{
    a=0;do c=getchar();while(c<'0'||c>'9');
    while(c<='9'&&c>='0')a=(a<<3)+(a<<1)+c-'0',c=getchar();
}
const int Maxn=10000001; 
int Prime[Maxn],tot,G[Maxn],Pow[Maxn],Maxp[Maxn];
bool Check[Maxn];
ll ans; 
inline int min(int a,int b){return a<b?a:b;} 
int main()
{ 
  ll k;
  int T,n,m,i,j;
  for(i=2;i<Maxn;i++)
   { 
    if(!Check[i])Prime[++tot]=i,G[i]=Pow[i]=1,Maxp[i]=i; 
    for(j=1;j<=tot;j++)
    { 
      k=i*Prime[j];
      if(k>=Maxn)break; 
      Check[k]=1; 
      if(i%Prime[j])
      { 
        Pow[k]=1,Maxp[k]=Prime[j]; 
        if(Pow[i]==1)G[k]=-G[i]; 
      }
       else
       { 
         Pow[k]=Pow[i]+1,Maxp[k]=Maxp[i]*Prime[j],n=i/Maxp[i]; 
          if(n==1)G[k]=1;
         else G[k]=Pow[n]==Pow[k]?-G[n]:0; 
        break; 
       } 
    } 
  } 
  for(i=2;i<Maxn;i++)
    G[i]+=G[i-1]; 
  read(T);
  while(T--)
  {  
    read(n),read(m);
    for(ans=0,i=1;i<=n&&i<=m;i=j+1)
     j=min(n/(n/i),m/(m/i)),ans+=(ll)(G[j]-G[i-1])*(n/i)*(m/i); 
    printf("%lld\n",ans); 
  }

  return 0; 
}