ACdream 1071

ACdream 1071
神奇的%系列一
Time Limit: 6000/3000MS (Java/Others) Memory Limit: 65536/32768KB (Java/Others)
Submit Statistic Next Problem
Problem Description

在计算机的世界里，%不是百分比，而是除法取余哟！
比如：
4 % 2 = 0
5 % 3 = 2
给你 2 ≤ N ≤ 100000 个数，a[1],a[2]…a[i]…a[n] (1 ≤ a[i] ≤ 100000)。
问有几个组合 (a[i], a[j])，(i != j, a[i] > a[j])，使得 a[i] % a[j] != 0。
Input

输入有多组数据。(<= 30)
对于每组数据：
第一行：N（表示 N 个数）
第二行：N 个元素 a[i]
Output

输出有几个组合 (a[i],a[j])，使得 a[i] % a[j] != 0
Sample Input

3
1 1 1
4
1 2 3 4
5
1 2 2 4 6
Sample Output

0
2
1

题意描述: 对于一个序列 如果存在有多少组 a[i],a[j] 使得 a[j] > a[i] 并且a[j] 不是 a[i] 的倍数。
我的思路：对于长度为 n的序列我们求出 a[j]>a[i] 并且 a[j] 是a[i] 的倍数的组数ans，然后用 sum = 1+2+3+…+n-1; sum - ans就是答案了。比如对于1 2 2 4 6这组数两两配对 共有 12 12 14 16 22 24 26 24 26 46 其中符合题意的是 46 这一组情况，我们的 ans 就是求出 诸如14 这种 a[j] 是 a[i] 倍数的的组数，其中要注意 22 也是满足 ans定义的。 好了下面贴代码

#include <string.h>#include <stdio.h>#include <cstdlib>using namespace std;const int N = 100000+100;typedef long long ll;int arr[N],book[N];int cmp(const void*a , const void*b){    return *(int*)a - *(int*)b;}int main(){   int i,j,n;   while(~scanf("%d",&n))   {       int inf = -1;       int nn = 0;       memset(book,0,sizeof(book));       for(i=1; i<=n; i++)       {           scanf("%d",&arr[i]);           book[arr[i]]++; //对原数组离散化，我们将用book数组求值           if(arr[i] > inf)             inf = arr[i];//inf 是你输入数据的最大值       }       ll ans = 0;       qsort(arr,n,sizeof(arr[0]),cmp);       for(i=1; i<=inf; i++)       {           if(book[i] == 0)//如果我们没有输入i，则不进行下去。              continue;           ll sum = 0;           for(j=i*2; j<=inf; j+=i)//这里对book数组进行遍历，如果遇到i的倍数，sum加一           {               if(book[j])                 sum += book[j];           }           if(book[i] > 1)//若输入了多个i，比如两个2，则还有22 这种情况             ans += (book[i]-1)*book[i]/2;           ans += sum*book[i];详见24 26 24 26       }       ll ans1 = -ans;       for(i=n-1; i>=1; i--)       {           ans1 += i;       }       printf("%lld\n",ans1);   }   return 0;}

这个实属个人瞎想，大家要是有更好的写法望联系本弱。