洛谷 【P1069 细胞分裂】 {求解质因子} 奋斗的珂珂~

题目描述

Hanks 博士是BT (Bio-Tech，生物技术) 领域的知名专家。现在，他正在为一个细胞实验做准备工作：培养细胞样本。

HanksHanks 博士手里现在有 N种细胞，编号从 1−N，一个第 i种细胞经过 1 秒钟可以分裂为Si个同种细胞（Si为正整数）。现在他需要选取某种细胞的一个放进培养皿，让其自由分裂，进行培养。一段时间以后，再把培养皿中的所有细胞平均分入M个试管，形成M份样本，用于实验。Hanks 博士的试管数M很大，普通的计算机的基本数据类型无法存储这样大的M值，但万幸的是，M 总可以表示为m1的m2次方，即M = m1^{m2}，其中 m1,m2均为基本数据类型可以存储的正整数。

注意，整个实验过程中不允许分割单个细胞，比如某个时刻若培养皿中有 4个细胞，

Hanks博士可以把它们分入 2 个试管，每试管内2 个，然后开始实验。但如果培养皿中有5个细胞，博士就无法将它们均分入2 个试管。此时，博士就只能等待一段时间，让细胞们继续分裂，使得其个数可以均分，或是干脆改换另一种细胞培养。

为了能让实验尽早开始，Hanks博士在选定一种细胞开始培养后，总是在得到的细胞“刚好可以平均分入 M个试管”时停止细胞培养并开始实验。现在博士希望知道，选择哪种细胞培养，可以使得实验的开始时间最早。

输入格式

第一行,有一个正整数 N，代表细胞种数。

第二行,有两个正整数 m1,m2 ，以一个空格隔开，即表示试管的总数 M = m1^{m2} .

第三行有 N 个正整数，第 i 个数 Si表示第 i 种细胞经过 1 秒钟可以分裂成同种细胞的个数。

输出格式

一个整数，表示从开始培养细胞到实验能够开始所经过的最少时间（单位为秒）。

如果无论Hanks博士选择哪种细胞都不能满足要求，则输出整数−1。

输入输出样例

输入 #1复制
1
2 1
3
输出 #1复制
-1
输入 #2复制
2
24 1
30 12
输出 #2复制
2

说明/提示

【输入输出说明】

经过 1秒钟，细胞分裂3 个，经过2秒钟，细胞分裂成9个，……，可以看出无论怎么分裂，细胞的个数都是奇数，因此永远不能分入 2个试管。

【输入输出样例2说明】

第 1 种细胞最早3 秒后才能均分入24 个试管，而第2种最早在2 秒后就可以均分（每试管144/(241)=6 个）。故实验最早可以在2 秒后开始。

【数据范围】

对于 50%的数据，有m1^{m2} ≤ 30000 ≤30000。

对于所有的数据，有1 ≤N≤ 10000,1 ≤m1 ≤ 30000,1 ≤m2 ≤ 10000,1 ≤ Si ≤ 2,000,000,000
NOIP 2009 普及组 第三题

解题思路

学习了别人的题解，具体见代码注解。

完整代码

#include<bits/stdc++.h>
using namespace std;

struct node{//使用结构体进行存储 
	int cnt,yin[30005],ci[30005];//cnt代表质因数的个数，yin代表质因数，ci代表每个质因数的个数 
}p,q;//p存储m1的分解，q存储m2的分解 

void fenjie(int x,node &t)//将x分解成质因子，并记录每一个质因子的个数 
{
	t.cnt=0;
	for(int i=2;i<=x/i;i++)
	{
		if(x%i==0)
		{
			t.yin[++t.cnt]=i;
			t.ci[t.cnt]=0;//还没有进行除法，先赋值为0 
			while(x%i==0)
			{
				x/=i;
				t.ci[t.cnt]++;
			}
		}
	}
	if(x>1)
	{
		t.yin[++t.cnt]=x;
		t.ci[t.cnt]=1;
	}
 } 
 
 
int main()
{
	int N;
	int m1,m2;
	scanf("%d",&N);
	scanf("%d %d",&m1,&m2);//m1的质因子一定存在在x的质因子中，x的质因子个数可以大于m1 
	
	if(m1==1) //如果是1，都可以直接放入，则所需要时间是0 
	 {
	  printf("0");
	  return 0;
     }
	
	fenjie(m1,p);
	
	int ans=-1,x;
	
	for(int i=1;i<=N;i++)
	{
		scanf("%d",&x);
		fenjie(x,q);
		
		int maxn=0,nxt=1;
		//maxn代表细胞分裂的次数，nxt代表当前处理的质因子的对应结构体下标
		
		bool  flag=false;
		
		for(int j=1;j<=p.cnt;j++)
		{
			while(q.yin[nxt]<p.yin[j]&&nxt<=q.cnt)  nxt++;
			//寻找对应的质因子是否存在 
			if(q.yin[nxt]>p.yin[j]||nxt>q.cnt) break;
			//x的质因子中不包含m1的质因子，就说明不满足
			
			int f=p.ci[j]*m2/q.ci[nxt];//*m2的原因是本来他就要m1^m2,那么它的质因子的次数也要相应的乘以m2 
			if(p.ci[j]*m2%q.ci[nxt])  f++;//（接上行），因为两数相乘，次幂是相加的 。 
			//求解细胞需要分裂的次数//除不尽就要多分裂一次
			if(maxn<f)  maxn=f;
			//因为x需要包含每一个质因数，且分裂的次数一定是分裂最多的那一次
			//才能满足
			
			if(j==p.cnt)
			   flag=true;//说明完全满足要求 ，既没有中途退出，又满足m1的质因子都在x的质因子中 
		
		 } 
		 if(flag&&(ans==-1||ans>maxn)) ans=maxn;
		 //因为是多次输入，可能出现这是第一次满足的情况 
		 // 也可能出现以前有情况，但是现在存在更优解 
	 } 
	 printf("%d",ans);
	 return 0; 
 }