zoukankan      html  css  js  c++  java
  • P2257 莫比乌斯+整除分块

    #include<bits/stdc++.h>
    #define ll long long
    using namespace std;
    const int maxn=1e7+10;
    int vis[maxn];
    int mu[maxn];
    int prime[maxn];
    int tot=0;
    int sum1[maxn];
    int sum2[maxn];
    void get_mu()
    {
        mu[1]=1; vis[1]=1;
        for(int i=2;i<maxn;i++)
        {
            if(!vis[i]) {mu[i]=-1; prime[++tot]=i; }
            for(int j=1;j<=tot && i*prime[j]<maxn;j++)
            {
                vis[i*prime[j]]=1;
                if(i%prime[j]==0) break;
                mu[i*prime[j]]=-mu[i];
            }
        }
        for(int i=1;i<=tot;i++)
            for(int j=prime[i];j<maxn;j+=prime[i])
                sum1[j]+=mu[j/prime[i]];
        for(int i=1;i<=maxn;i++)
            sum2[i]=sum2[i-1]+sum1[i];
    }
    int main()
    {
        get_mu();
        int T; cin>>T;
        while(T--)
        {
            int n,m; cin>>n>>m;
            ll ans=0;
            for(int l=1,r;l<=min(n,m);l=r+1)
            {
               // r=min(n,m)/(min(n,m)/l);  // l-r;
               r=min( n/(n/l),m/(m/l));
                ans+=1ll*(n/l)*(m/l)*(sum2[r]-sum2[l-1]);
            }
            /*int t=0;
            for(int i=1;i<=n;i++)
            {
                for(int j=1;j<=m;j++)
                {
                    if(vis[__gcd(i,j)]==0) t++;
                }
            }*/
            cout<<ans<<endl;
        }
    }
  • 相关阅读:
    wm
    usual
    itk_option remove
    Label Options
    imosflm controller
    set font
    namespace code
    git 之五分钟教程
    git 之五分钟教程
    学习perl正则表达式
  • 原文地址:https://www.cnblogs.com/Andromeda-Galaxy/p/10702417.html
Copyright © 2011-2022 走看看