【BZOJ】【4146】 【AMPPZ2014】Divisors

暴力

---

　　由于值的范围很小（$ \leq 2*10^6$），所以用一个cnt数组统计每个值有多少个数，然后从小到大，统计每个数的倍数即可。

　　根据调和数？的神奇性质= =这样是$O(nlogn)$的……

 /**************************************************************
     Problem: 4146
     User: Tunix
     Language: C++
     Result: Accepted
     Time:4644 ms
     Memory:9084 kb
 ****************************************************************/

 //BZOJ 4146
 #include<vector>
 #include<cstdio>
 #include<cstring>
 #include<cstdlib>
 #include<iostream>
 #include<algorithm>
 #define rep(i,n) for(int i=0;i<n;++i)
 #define F(i,j,n) for(int i=j;i<=n;++i)
 #define D(i,j,n) for(int i=j;i>=n;--i)
 #define pb push_back
 using namespace std;
 inline int getint(){
     int v=,sign=; char ch=getchar();
     while(ch<''||ch>''){ if (ch=='-') sign=-; ch=getchar();}
     while(ch>=''&&ch<=''){ v=v*+ch-''; ch=getchar();}
     return v*sign;
 }
 const int N=2e6+,INF=~0u>>;
 typedef long long LL;
 /******************tamplate*********************/
 int n,c[N];
 int main(){
 #ifndef ONLINE_JUDGE
     freopen("4146.in","r",stdin);
     freopen("4146.out","w",stdout);
 #endif
     n=getint();
     int mx=;
     F(i,,n){
         int x=getint();
         c[x]++;
         mx=max(mx,x);
     }
     LL ans=;
     F(i,,mx) if (c[i]){
         ans+=(LL)c[i]*(c[i]-);
         for(int j=i+i;j<=mx;j+=i) ans+=(LL)c[i]*c[j];
     }
     printf("%lld\n",ans);
     return ;
 }
 

4146: [AMPPZ2014]Divisors

Time Limit: 20 Sec  Memory Limit: 256 MB
Submit: 12  Solved: 6
[Submit][Status][Discuss]

Description

给定一个序列a[1],a[2],...,a[n]。求满足i!=j且a[i]|a[j]的二元组(i,j)的个数。

Input

第一行包含一个正整数n(1<=n<=2000000)，表示序列长度。

第二行包含n个正整数，依次表示a[1],a[2],...,a[n](1<=a[i]<=2000000)。

Output

一个整数，即满足条件的二元组的个数。

Sample Input

5
2 4 5 2 6

Sample Output

6

HINT

满足条件的6组分别为(1,2),(1,4),(1,5),(4,1),(4,2),(4,5)。

Source

鸣谢Claris上传

[Submit][Status][Discuss]