LYDSY模拟赛day2 Divisors

/*
注意分解质因数,如果i是,那么n/i也是,这样就可以解决分解质因数的时间问题了
当 k ≥ 1 时,只有这些数的约数才会对答案产生贡献。
求出 m 个数的所有不超过 n 的约数,去重后统计即可。
求出 k = 1 到 m 的所有答案后,剩下的数字个数就是 k = 0
的答案。
时间复杂度 O(m2√a)
*/
#include<cstdio>
#include<algorithm>
const int N=205;
int n,m,i,a[N],cnt,b[1000000],f[N];
void add(int n){for(int i=1;i<=n/i;i++)if(n%i==0)b[++cnt]=i,b[++cnt]=n/i;}
int cal(int d){
  int t=0;
  for(int i=1;i<=m;i++)if(a[i]%d==0)t++;
  return t;
}
int main(){
  freopen("div.in","r",stdin);freopen("div.out","w",stdout);
  scanf("%d%d",&n,&m);
  for(i=1;i<=m;i++)scanf("%d",&a[i]),add(a[i]);
  std::sort(b+1,b+cnt+1);
  for(i=1;i<=cnt;i++)if(b[i]<=n&&b[i]!=b[i-1])f[cal(b[i])]++;
  f[0]=n;
  for(i=1;i<=m;i++)f[0]-=f[i];
  for(i=0;i<=m;i++)printf("%d
",f[i]);
  fclose(stdin);fclose(stdout);
  return 0;
}
原文地址:https://www.cnblogs.com/hyfer/p/5928039.html