ZZULI OJ 1902: 985的因子对难题

1902: 985的因子对难题

Time Limit: 2 Sec   Memory Limit: 128 MB
Submit: 116   Solved: 23

Submit Status Web Board

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

2
5
1 2 3 4 5
5
2 2 2 2 2

Sample Output

5
10

HINT

Source

hpu

。。。。。。暴力来一遍。。。。。

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
int prim[1000010]={1,1};
int t,n,a,p[1001000],b[1001000];
int main()
{
	int i,j;
	scanf("%d",&t);
	while(t--)
	{
		memset(p,0,sizeof(p));
		memset(b,0,sizeof(b));
		scanf("%d",&n);
		getprim();
		for(i=0;i<n;i++)
		{
			scanf("%d",&a);
			p[a]++;//统计每个数出现的次数 
		}
		for(i=1;i<=1000000;i++)//打印出质因数的个数 
		{
			if(p[i])
			{
				for(j=2*i;j<=1000000;j=j+i)
				b[j]+=p[i];
			}
		}
		long long ans=0;
		for(i=1;i<=1000000;i++)
		{
			if(p[i])
			{
					if(p[i]>1)
				    ans+=(p[i]*(p[i]-1))/2;
				    ans+=b[i]*p[i];
			}
		}
		printf("%lld\n",ans);
	}
	return 0;
}