今天是贝茜的生日，为了庆祝自己的生日，贝茜邀你来玩一个游戏．
贝茜让N(1≤N≤100000)头奶牛坐成一个圈．除了1号与N号奶牛外，i号奶牛与i-l号和i+l号奶牛相邻．N号奶牛与1号奶牛相邻．农夫约翰用很多纸条装满了一个桶，每一张包含了一个独一无二的1到1,000,000的数字．
接着每一头奶牛i从柄中取出一张纸条Ai．每头奶牛轮流走上一圈，同时拍打所有编号能整除在纸条上的数字的牛的头，然后做回到原来的位置．牛们希望你帮助他们确定，每一头奶牛需要拍打的牛．
比较水。。
维护出所有的颜色。
然后按埃氏筛的方法挨个筛就行了
根据调和级数，时间复杂度为$O(nlogn$
#include<cstdio> #include<algorithm> #include<cstring> using namespace std; const int MAXN = 1e5 + 10, INF = 1e9 + 10; inline int read() { char c = getchar(); int x = 0, f = 1; while(c < ‘0‘ || c > ‘9‘) {if(c == ‘-‘) f = -1; c = getchar();} while(c >= ‘0‘ && c <= ‘9‘) x = x * 10 + c - ‘0‘, c = getchar(); return x * f; } int a[1000001], color[MAXN], ans[1000001]; main() { #ifdef WIN32 freopen("a.in", "r", stdin); #endif int N = read(); for(int i = 1; i <= N; i++) color[i] = read(), a[color[i]]++; for(int i = 1; i <= 1000000; i++) { if(!a[i]) continue; for(int j = i; j <= 1000000; j += i) ans[j] += a[i]; } for(int i = 1; i <= N; i++) printf("%d\n", ans[color[i]] - 1); }
BZOJ1607: [Usaco2008 Dec]Patting Heads 轻拍牛头(模拟 调和级数)
原文：https://www.cnblogs.com/zwfymqz/p/9226499.html