P1414 又是毕业季II

思路：

一开始很容易想到枚举n个数取k个的所有组合，然后分别用辗转相除法求最大公约数，但是复杂度明显不符合要求，于是必须换一种思路。

我们想到，k个数的公约数含义就是这k个数均含有某个因数，如果我们把所有数的因数全部求出来，发现有k个数均含有某个因数，那么这个数必然是这k个数的公约数。其中找出最大的就是它们的最大公约数。但是要如何高效的做到这点呢？考虑到对于k=1，2……，n都要求出，我们可以这么做：

* 1、 求出每个因数出现的次数。

* 2、 对于每个次数记录最大的因数。

* 3、 根据f[k]=max(f[k],f[k+1])逆向递推。（如果已经知道k个数的最大公约数是m，那么l(l<k)个数的最大公约数一定大于等于m）。

我们想到，k个数的公约数含义就是这k个数均含有某个因数，如果我们把所有数的因数全部求出来，发现有k个数均含有某个因数，那么这个数必然是这k个数的公约数。其中找出最大的就是它们的最大公约数。
每个数分解因数，c[i]表示i作为因子的次数。 对于答案i，c[x]>i的p可以作为答案。
可以发现i的答案一定大于等于i+1的答案。
扫一遍找答案行了。

ACocde:

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+10;
int cnt[N],n,x,t;
void solve() {
  cin>>n;
  for(int i=1;i<=n;i++){
  	cin>>x;
  	t=max(t,x);//小优化，找出最大数，这样子后面就不用从最大范围开始找了
	 //把x的所有约数出现的次数都加一  
  	for(int j=1;j<=sqrt(x);j++){//枚举标记x的约数 
  		if(x%j==0){
  			 cnt[j]++;
            if(j*j!=x) cnt[x/j]++;//如果J*j==x，那么只能算一个约数 
		  }
  		
	  } 
  }
  for(int i=1;i<=n;i++){
  	while(cnt[t]<i)t--;
  	cout<<t<<"\n"; 
  }
}
int main() {

	ios::sync_with_stdio(false);
	int tt=1;
	//cin>>tt;
	while(tt--)
		solve();
	return 0;
}

over~