• NOIP 模拟 $87; m 肯德基$


    题解 (by;zjvarphi)

    推一下式子:

    [f(x)=mu(x)^2x\ mu(x)^2=sum_{d^2|x}mu(d)\ f(x)=xsum_{d^2|x}mu(d)\ ]

    原式就等于:

    [sum_{i=1}^nisum_{d^2|i}mu(d) ]

    调整枚举顺序:

    [sum_{d=1}^{sqrt n}sum_{i=1}^{frac{n}{d^2}}mu(d)d^2i\ sum_{d=1}^{sqrt n}mu(d)d^2sum_{i=1}^{frac{n}{d^2}}i\ ]

    注意,文中的 (frac{n}{d^2}) 都是向下取整。

    这里来证明一下为什么 (mu(x)^2=sum_{d^2|x}mu(d))

    如果 (x) 中没有平方因子,那么 (mu(x)^2)(1)(d=1) 时恰好贡献了 (1)

    如果有,那么考虑通过容斥的角度证明。

    假设 (x) 中有 (k) 个平方质因子,那么选一个的时候被多减去了,在选两个的时候补回来,可以画一个维恩图理解一下。

    其实也就是 (sum_{i=0}^n(-1)^iinom{n}{i}),其实就是二项式定理。

    化完的式子直接整除分块即可,复杂度 (sqrt n+Tsqrt[3] n)

    Code
    #include<bits/stdc++.h>
    #define ri signed
    #define pd(i) ++i
    #define bq(i) --i
    #define func(x) std::function<x>
    namespace IO{
        char buf[1<<21],*p1=buf,*p2=buf;
        #define gc() p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?(-1):*p1++
        #define debug1(x) std::cerr << #x"=" << x << ' '
        #define debug2(x) std::cerr << #x"=" << x << std::endl
        #define Debug(x) assert(x)
        struct nanfeng_stream{
            template<typename T>inline nanfeng_stream &operator>>(T &x) {
                bool f=false;x=0;char ch=gc();
                while(!isdigit(ch)) f|=ch=='-',ch=gc();
                while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=gc();
                return x=f?-x:x,*this;
            }
        }cin;
    }
    using IO::cin;
    namespace nanfeng{
        #define FI FILE *IM
        #define FO FILE *OUT
        template<typename T>inline T cmax(T x,T y) {return x>y?x:y;}
        template<typename T>inline T cmin(T x,T y) {return x>y?y:x;}
        using ull=unsigned long long;
        static const int N=1e7+7;
        int mu[N],prim[N],cnt,T;
        bool vis[N];
        ull sum[N],a;
        auto Getmu=[]() {
            const int n=N-7;
            mu[1]=1;
            for (ri i(2);i<=n;pd(i)) {
                if (!vis[i]) mu[prim[++cnt]=i]=-1;
                for (ri j(1);j<=cnt&&1ll*i*prim[j]<=n;pd(j)) {
                    const int cur=i*prim[j];
                    vis[cur]=true;
                    if (!(i%prim[j])) break;
                    mu[cur]=-mu[i];
                }
            }
            for (ri i(1);i<=n;pd(i)) sum[i]=sum[i-1]+1ull*mu[i]*i*i;
        };
        auto cal=[](ull n) {return n&1?n*((n+1)/2):(n/2)*(n+1);};
        auto calc=[](ull x) {
            ull res=0,lim=sqrt(x);
            for (ull l(1),r;l<=lim;l=r+1) {
                r=cmax(l,(ull)sqrt(x/(x/l/l)));
                res+=(sum[r]-sum[l-1])*cal(x/l/l);
            }
            return res;
        };
        inline int main() {
            FI=freopen("kfc.in","r",stdin);
            FO=freopen("kfc.out","w",stdout);
            Getmu();
            cin >> T;
            for (ri i(1);i<=T;pd(i)) cin >> a,printf("%llu
    ",calc(a));
            return 0;
        }
    }
    int main() {return nanfeng::main();}
    
  • 相关阅读:
    对json按某个键的值进行排序(转载)
    nrm -- NPM registry 管理工具
    easyui-filebox 文件上传
    Netty,Thrifty
    VS2015 framework4.5代码提示英文切换为中文
    .NET Runtime version 2.0.50727.8762
    js moment.js日期操作类 datetime,日期操作,dayjs
    Oracle.ManagedDataAccess.dll
    打造自己的JavaScript武器库(转)
    Vue.js实战 5.5章 购物车
  • 原文地址:https://www.cnblogs.com/nanfeng-blog/p/15494388.html
Copyright © 2020-2023  润新知