zzuli 1902 985的因子对难题

Description

985有n个正整数，他想知道存在多少个不同的因子对(a[i], a[j])使得

1 <= i, j <= n && i != j && a[j] % a[i] == 0，其中i和j是元素的下标。

特别地，他认为（a[i]，a[j]）与（a[j]，a[i]）是一样的因子对。

Input

第一行输入一个整数t，代表有t组测试数据。

每组数据占两行，第一行输入一个n代表元素个数，下面一行输入n个整数a[]。

注：1 <= t <= 30，1 <= n <= 1e5，1 <= a[] <= 1e6。

Output

一个整数代表最后的答案。

Sample Input

251 2 3 4 552 2 2 2 2

Sample Output

510

题解：

因子对个数 = 该数个数 * 该数倍数的个数 + 该数的个数的组合数；

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <stack>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#include <time.h>
using namespace std;

#define LL long long
#define N 100010
#define mod 1000000007
int a[N];
int main()
{
   int T,n,x;
   scanf("%d",&T);
   while(T--)
   {
       scanf("%d",&n);
       int maxx=0,ans=0;
       memset(a,0,sizeof(a));
       for(int i=1;i<=n;i++)
       {
           scanf("%d",&x);
           a[x]++;
           if(x>maxx) maxx=x;
       }
       for(int i=1;i<=maxx;i++)
       {
           if(a[i])
           {
               for(int j=i*2;j<=maxx;j+=i)
                 ans+=a[i]*a[j];
               ans+=a[i]*(a[i]-1)/2;
           }
       }
       printf("%d\n",ans);
   }
}