【AMPPZ2014】【BZOJ4146】Divisors

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],…,an。

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上传

计数器+枚举因子
调和级数复杂度
(但是我跑的好慢

#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#include<algorithm>#define MAXN 2000010#define LL long long#define GET (ch>='0'&&ch<='9')using namespace std;int n,maxn;int a[MAXN],c[MAXN],f[MAXN];LL ans;void in(int &x){    char ch=getchar();x=0;    while (!GET)    ch=getchar();    while (GET) x=x*10+ch-'0',ch=getchar();}int main(){    in(n);    for (int i=1;i<=n;i++)  in(a[i]),c[a[i]]++,maxn=max(maxn,a[i]);    for (int i=1;i<=maxn;i++)        for (int j=i;j<=maxn;j+=i)  f[j]+=c[i];    for (int i=1;i<=maxn;i++)   ans+=(LL)(c[i])*f[i];    ans-=n;cout<<ans<<endl;}