zoukankan      html  css  js  c++  java
  • CodeForces

    http://codeforces.com/problemset/problem/803/F

    #include <iostream>
    #include <cstdio>
    #include <queue>
    #include <algorithm>
    #include <cmath>
    #include <cstring>
    #define inf 2147483647
    #define N 1000010
    #define mod 1000000007
    #define p(a) putchar(a)
    #define For(i,a,b) for(long long i=a;i<=b;++i)
    //by war
    //2019.8.21
    using namespace std;
    long long n,a[N],m=1e5,tot[N],num,cnt,ans;
    long long prime[N],mu[N];
    bool vis[N];
    void in(long long &x){
        long long y=1;char c=getchar();x=0;
        while(c<'0'||c>'9'){if(c=='-')y=-1;c=getchar();}
        while(c<='9'&&c>='0'){ x=(x<<1)+(x<<3)+c-'0';c=getchar();}
        x*=y;
    }
    void o(long long x){
        if(x<0){p('-');x=-x;}
        if(x>9)o(x/10);
        p(x%10+'0');
    }
    
    void Euler(){
        mu[1]=1;
        For(i,2,1e5){
            if(!vis[i]) prime[++cnt]=i,mu[i]=-1;
            for(long long j=1;j<=cnt&&i*prime[j]<=1e5;j++){
                vis[i*prime[j]]=1;
                if(i%prime[j]==0){
                    mu[i*prime[j]]=0;
                    break;
                }
                mu[i*prime[j]]=-mu[i];
            }
        }
    }
    
    long long ksm(long long a,long long b){
        long long r=1;
        while(b>0){
            if(b&1)
                r=r*a%mod;
            a=a*a%mod;
            b>>=1;
        }
        return r;
    }
    
    signed main(){
        in(n);
        Euler();
        For(i,1,n)
            in(a[i]),tot[a[i]]++;
        For(i,1,m){
            num=0;
            for(long long k=1;k*i<=m;k++)
                num+=tot[k*i];
            ans+=mu[i]*(ksm(2,num)-1);
            ans%=mod; 
        }
        o((ans+mod)%mod);
        return 0;
    }
  • 相关阅读:
    JS 数组
    JS 模拟彩票
    C++ 动态内存
    计算机网络--OSI七层模型
    C++ 异常处理
    C++ 文件和流
    数据库学习教程网站
    数据结构--哈夫曼树
    数据结构--红黑树
    数据结构--伸展树
  • 原文地址:https://www.cnblogs.com/war1111/p/11390900.html
Copyright © 2011-2022 走看看