1621：轻拍牛头

1621：轻拍牛头原题(嵌入代码可能效果不好，请点链接看原题)

【题目描述】
原题来自：USACO 2008 Dec. Silver

今天是贝茜的生日，为了庆祝自己的生日，贝茜邀你来玩一个游戏。

贝茜让 NN 头奶牛坐成一个圈。除了 11 号与 NN 号奶牛外，ii 号奶牛与 i−1i−1 号和 i+1i+1 号奶牛相邻，NN 号奶牛与 11 号奶牛相邻。农夫约翰用很多纸条装满了一个桶，每一张包含了一个 11 到 106106 的数字。

接着每一头奶牛 i 从桶中取出一张纸条 Ai ，每头奶牛轮流走一圈，同时拍打所有「编号是 AiAi 的约数」的牛，然后走回到原来的位置。牛们希望你帮助他们确定，每一头奶牛需要拍打的牛。

【输入】
第一行包含一个整数 NN；

接下来第二到第 N+1N+1 行每行包含一个整数 AiAi 。

【输出】
第一到第 N 行，第 i 行的输出表示第 i 头奶牛要拍打的牛数量。

【输入样例】
5
2
1
2
3
4
【输出样例】
2
0
2
1
3
【提示】
数据范围与提示：

对于全部数据，1≤N≤1051≤N≤105 。

【题意说明】

　　对这个题目的描述，我觉得有点描述不清。”每头奶牛轮流走一圈，同时拍打所有[编号是Ai的约数]的牛“这里的编号指向不明确，如果是最初的1-N编号，显然样例数据不符。因为2号奶牛拿的数字是1，那它至少应该拍1号奶牛，但样例输出是0。查看另一篇博文https://www.luogu.com.cn/problem/P2926，描述为”每头奶牛轮流走上一圈，同时拍打所有手上数字能整除在自己纸条上的数字的牛的头“，这个描述很清楚，也符合样例数据，以下内容均按后者理解。

【分析】

　　这个题目看起来并不难，判断一个数是否是另一个数的约数取余就可以了，要计算拍牛的次数，加一个计数器就OK。所以轻松搞定代码。

#include<iostream>
using namespace std;
int const N=1e6+5;
int a[N],ans[N],n; 
int main(){
    cin>>n;
    for(int i=1;i<=n;i++)cin>>a[i];
    for(int i=1;i<=n;i++)//枚举每一头牛 
        for(int j=1;j<=n;j++)//其他牛的编号是否为当前牛编号的约数 
            if(a[i]%a[j]==0)ans[i]++;
    for(int i=1;i<=n;i++)cout<<ans[i]-1<<endl;//ans[i]-1是要减去本身 
    return 0;
}

　　样例顺利通过，提交网站结果应该可以想到。这样都能通过是不是也太没水平了。简单算一个，对于百分之百的数据n<=10^5，一个双循环10^10，超时是必然的。

　　在之前一个题1619：用到了筛选法可以优化时间复杂度，那我们也可以考虑用筛选法优化。我们不用把每一个编号都去验证，只把当前编号的倍数标记一次。那么每一头牛的编号由原来的验证其余牛的时间复杂度n可降为最大编号值/当前编号值。代码实现

//1621：轻拍牛头 
#include<iostream>
using namespace std;
int const N=1e6+5;
int a[N],b[N],n,c[N],maxn=0; 
int main(){
    cin>>n;
    for(int i=1;i<=n;i++)
    {
        cin>>a[i];
        if(maxn<a[i])maxn=a[i];
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=a[i];j<=maxn;j+=a[i])
            b[j]++;
    }
    for(int i=1;i<=n;i++)
        cout<<b[a[i]]-1<<endl;
    return 0;
}

　　再次提交，46分，其余超时。想想一下，这个题数据量达到10^5，输入和输出应该占时较大，所以cin/cout的效率让我们不得不放入考虑范围，换吧，scanf/printf的效率可以提高不少，在前面博文中也有提及，特别是1205题就是一个很好例子。得以下代码

//1621：轻拍牛头 
#include<cstdio>
int const N=1e6+5;
int a[N],b[N],n,maxn=0; 
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",a+i);
        if(maxn<a[i])maxn=a[i];
    }
    for(int i=1;i<=n;i++)
        for(int j=a[i];j<=maxn;j+=a[i])
            b[j]++;
    for(int i=1;i<=n;i++)
        printf("%d
",b[a[i]]-1);
    return 0;
}

　　再次提交，84分。这还能优化？必须啊，别人能过啊，那题肯定没问题啊。思考良久，想到一个地方：编号如果有重复，那我们验证时就会重复验证，特别编号又小，重复量再大一点，这一环节的效果更明显，再挤一挤吧。把重复项合并，只筛选一次就好。代码如下

//1621：轻拍牛头 
#include<cstdio>
int const N=1e6+5;
int a[N],b[N],c[N],d[N][2],cnt,n,maxn,x; 
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",a+i);
        c[a[i]]++;
        if(maxn<a[i])maxn=a[i];
    }
    for(int i=1;i<=maxn;i++)
        if(c[i])d[++cnt][0]=i,d[cnt][1]=c[i];
    for(int i=1;i<=cnt;i++)
        for(int j=d[i][0];j<=maxn;j+=d[i][0])
            b[j]+=d[i][1];
    for(int i=1;i<=n;i++)
        printf("%d
",b[a[i]]-1);
    return 0;
}

　　再次提交：100分。终于AC了！

如果大家有更好看法，欢迎评论！！！