ACdream 1071 神奇的%系列一 筛选法

题目：http://acdream.info/problem?pid=1071

题意：

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

思路：用类似用筛素数法的筛选法求出对于每个a[j]的a[i] % a[j] == 0的个数，然后大于a[j]的区间长度减去a[i] % a[j] == 0的个数，这就是每个对于每个a[j]来说a[i] % a[j] != 0的个数，然后再乘以a[j]的个数，最后累加一下就好了，注意会爆int

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

typedef long long ll;
const int N = 100010;
int arr[N], mark[N];
int main()
{
    int n;
    while(~ scanf("%d", &n))
    {
        memset(mark, 0, sizeof mark);
        for(int i = 1; i <= n; i++)
            scanf("%d", arr + i);
        sort(arr + 1, arr + 1 + n); 
        for(int i = 1; i <= n; i++)
            mark[arr[i]]++; //统计每个数字出现的个数

        ll res = 0;
        for(int i = 1; i <= n; i++)
            if(arr[i] != arr[i-1]) //对于每个数字只筛选一次，相同的数字不再筛选
            {
                int tmp = 0; //tmp保存大于arr[i]且是其倍数的数的个数
                for(int j = arr[i] * 2; j <= arr[n]; j += arr[i])
                        tmp += mark[j];
                res += (n - i - mark[arr[i]] + 1 - tmp) * mark[arr[i]]; //大于arr[i]的区间长度减去是其倍数的数的个数，再乘以arr[i]的个数
            }
        printf("%lld\n", res);
    }
    return 0;
}