【bzoj1607】【Usaco2008 Dec】Patting Heads (数学) 题解&代码

题目链接：
http://www.lydsy.com/JudgeOnline/problem.php?id=1607
题解：
逆向思考，要求这些数字中有多少个数字是其中一个数字x的约数，不如对于每个给出的数字，在最大数字即1000000的范围内枚举它的倍数，更新枚举到的倍数的答案，这样的复杂度大概是(1000000/a1)+(1000000/a2)+…….+(1000000/an)，注意到它会给出重复数字，如果全给2,2,2,2,2………2,2，很明显是超时的，那么我们就用cnt数组统计某一个数字出现的个数，枚举1到1000000中出现过的数，这样的话，最坏情况的数据是1,2,3,………100000,用程序运算一下，发现在这样的最坏情况下复杂度也只有12041067还是可以接受的。
代码：

#include<iostream>
#include<algorithm>
#include<stdio.h>
#define maxn (100005)
#define maxm (1000005)
using namespace std;
int n,m,a[maxn],vis[maxm],cnt[maxm];
int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        cnt[a[i]]++;
    }
    for (int i=1;i<=maxm;i++)
    if (cnt[i])
    {
        for (int j=i;j<=maxm;j+=i)
        vis[j]+=cnt[i];
    }   
    for (int i=1;i<=n;i++)
    {
        printf("%d\n",vis[a[i]]-1);
//减一是因为运算过程中多加了他自己
    }
}