[bzoj-1607][Usaco2008 Dec]Patting Heads 轻拍牛头 题解

题目传送门
题意解析：题目很明显，给你每一只cow的编号ai，然后问你每只奶牛是它的约数是有几只（不包括自己）。

---

My opinion：我还是小看了数据，一开始以为用O(n√n)可以卡过，没想到bzoj上的总时间只能有3秒，那么只能用O(nlgn)的算法了，那么很明显，对于每个a[i],它对答案的贡献是1，即对每个ans[j] (j%a[i]==0)答案加一，然后我们用类似筛法的写法就可以做到了。
总结：
1、输入，计算出当前a[i]的值有多少个。
2、计算每个数的倍数，并在这个数的倍数的ans上加上这个数的个数。
3、输出。

---

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define rep(i,a,n) for (int i=a;i<=n;i++)
#define per(i,a,n) for (int i=a;i>=n;i--)
#define Clear(a,x) memset(a,x,sizeof(a))
#define ll long long
#define db double
#define INF 2000000000
#define eps 1e-8
using namespace std;
int read(){
    int x=0,f=1;
    char ch=getchar();
    while (ch<'0'||ch>'9') f=ch=='-'?-1:f,ch=getchar();
    while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
const int maxn=100005,maxsum=1000005;
int n,num,Max;
int a[maxn],ans[maxsum],p[maxsum];
int main(){
    n=read();
    rep(i,1,n){
        a[i]=read();
        p[a[i]]++;
        Max=max(a[i],Max);
    }
    rep(i,1,Max)
        if (p[i])
            for (int j=i;j<=Max;j+=i)
                ans[j]+=p[i];
    rep(i,1,n)
        printf("%d\n",ans[a[i]]-1);
    return 0;
}