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;
}

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