题目背景

“叮铃铃铃”,随着高考最后一科结考铃声的敲响,三年青春时光顿时凝固于此刻。毕业的欣喜怎敌那离别的不舍,憧憬着未来仍毋忘逝去的歌。1000多个日夜的欢笑和泪水,全凝聚在毕业晚会上,相信,这一定是一生最难忘的时刻!

题目描述

彩排了一次,老师不太满意。当然啦,取每位同学的号数来找最大公约数显然不太合理。于是老师给每位同学评了一个能力值。于是现在问题变为,从n个学生中挑出k个人使得他们的默契程度(即能力值的最大公约数)最大。但因为节目太多了,而且每个节目需要的人数又不知道。老师想要知道所有情况下能达到的最大默契程度是多少。这下子更麻烦了,还是交给你吧~

PS:一个数的最大公约数即本身。

输入输出格式

输入格式:

第一行一个正整数n。

第二行为n个空格隔开的正整数,表示每个学生的能力值。

输出格式:

总共n行,第i行为k=i情况下的最大默契程度。

输入输出样例

输入样例#1: 复制

4
1 2 3 4
输出样例#1: 复制

4
2
1
1

说明

【题目来源】

lzn原创

【数据范围】

记输入数据中能力值的最大值为inf。

对于20%的数据,n<=5,inf<=1000

对于另30%的数据,n<=100,inf<=10

对于100%的数据,n<=10000,inf<=1e6

思路:求k个数的最大公约数,就是找一个最大因子它在n个数中出现了k次及以上,因此我们需要分解因子,统计个数。

代码:

 #include"bits/stdc++.h"
#define db double
#define ll long long
#define vec vector<ll>
#define Mt vector<vec>
#define ci(x) scanf("%d",&x)
#define cd(x) scanf("%lf",&x)
#define cl(x) scanf("%lld",&x)
#define pi(x) printf("%d\n",x)
#define pd(x) printf("%f\n",x)
#define pl(x) printf("%lld\n",x)
//#define rep(i, x, y) for(int i=x;i<=y;i++)
#define rep(i,n) for(int i=0;i<n;i++)
const int N = 1e6 + ;
const int mod = 1e9 + ;
const int MOD = mod - ;
const int inf = 0x3f3f3f3f;
const db PI = acos(-1.0);
const db eps = 1e-;
using namespace std;
int n,cnt;
int p;
int c[N],pri[N];
bool v[N];
int sum[N];
int a[],b[];
int s[N];
void init()//素数表
{
cnt = ;
for(int i = ; i*i < N; i++)
{
if(!v[i]) pri[cnt++] = i;
for(int j = ; j < cnt && pri[j] <= N/i; j++)
{
v[i*pri[j]] = ;
if(i % pri[j]==) break;
}
}
} void dfs(int id,int res){//dfs求每个因子的出现次数
if(id>=p) return;
dfs(id+,res);
for(int i=;i<b[id];i++){
res*=a[id];
sum[res]++;
dfs(id+,res);
}
}
void work(int x){//分解数字
p=;
memset(a,,sizeof(a));
memset(b,,sizeof(b));
for(int i=;i*i<=x;i++){
if(x%i==){
a[p]=i;
while(x%i==) b[p]++,x/=i;
p++;
}
}
if(x!=) a[p]=x,b[p]++,p++;
dfs(,);
} int main()
{
//freopen("data.in","r",stdin);
//freopen("data.out","w",stdout);
init();
ci(n);
int mx=-;
for(int i=;i<n;i++) ci(c[i]),mx=max(mx,c[i]),work(c[i]);
sum[]=n;
for(int i=;i<N;i++) s[i]=;
for(int i=;i<=mx;i++){
int x=sum[i];
if(s[x]<i) s[x]=i;
}
mx=;
for(int i=n;i>=;i--){
if(mx>s[i]) s[i]=mx;
mx=max(s[i],mx);
}
for(int i=;i<=n;i++) printf("%d\n",s[i]);
return ;
}

最新文章

  1. 如何挂载阿里云Linux服务器的“数据盘”(新购买)
  2. http详解
  3. 转 :meta name的含义:&lt;META http-equiv=Content-Type content=&quot;text/html; charset=gb2312&quot;&gt;
  4. Vmware web client 5.5 控制台连接不上:Connection timed out
  5. iOS JS交互
  6. FZU2143Board Game(最小费用流)
  7. Application,Session和Cookie
  8. BZOJ_1012_[JSOI2008]_最大数maxnumber_(线段树/树状数组+RMQ)
  9. 【动态规划】【最长公共子序列】Vijos P1111 小胖的水果
  10. Java实践 — SSH远程执行Shell脚本(转)
  11. Amaze UI 是一个移动优先的跨屏前端框架。 http://amazeui.org/
  12. C++\virtual 虚函数、纯虚函数
  13. 基于tcp的云盘上传下载的模拟
  14. thinkPHP3.2使用__ROOT__显示为空,为什么
  15. Centos 6.5使用vsftpd配置FTP服务器教程
  16. puppet(2)-资源介绍
  17. APP(通信)接口定义
  18. ElasticStack系列之二十 &amp; 数据均衡、迁移、冷热分离以及节点自动发现原理与机制
  19. ABP之事件总线(1)
  20. git learn

热门文章

  1. 基于vue-cli 将webpack3 升级到 webpack4 配置
  2. [RabbitMQ]Windows环境下rabbitmqclt(Command Line Tools)出现Erlang distribution failed错误的解决方法
  3. jQuery(四)(DOM一)
  4. 2013应届毕业生各大IT公司待遇整理汇总篇(转)
  5. 深入了解MongoDB
  6. MySQL安装和简单操作
  7. Android(java)学习笔记67:Android Studio新建工程中的build.gradle、settings.gradle
  8. 2017.11.2 JavaWeb----第六章 Servlet技术
  9. 移动页面请使用CSS3动画
  10. js中的变量提升(Hoisting)