POJ P1069 [NOIP2009 普及组] 细胞分裂【整除的联想】【RE的原因】【打表法分解质因数】【向上取整】

最新推荐文章于 2024-07-18 07:40:05 发布

bei2002315

最新推荐文章于 2024-07-18 07:40:05 发布

阅读量173

点赞数

分类专栏：数论错误总结文章标签：算法

本文链接：https://blog.csdn.net/bei2002315/article/details/125961450

版权

细胞分裂试管分配质因数分解生物技术实验数学模型

关键词由CSDN通过智能技术生成

数论同时被 2 个专栏收录

23 篇文章 0 订阅

订阅专栏

错误总结

5 篇文章 0 订阅

订阅专栏

题目描述

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

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

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

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

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

输入格式

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

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

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

输出格式

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

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

输入输出样例

输入 #1复制

1 
2 1 
3

输出 #1复制

-1

输入 #2复制

2
24 1
30 12

输出 #2复制

说明/提示

【输入输出说明】

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

【输入输出样例22说明】

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

【数据范围】

对于 50%的数据，有m_1^{m_2} ≤ 30000m1m2≤30000。

对于所有的数据，有1 ≤N≤ 10000,1 ≤m_1 ≤ 30000,1 ≤m_2 ≤ 10000,1 ≤ S_i ≤ 2,000,000,0001≤N≤10000,1≤m1≤30000,1≤m2≤10000,1≤Si≤2,000,000,000。

NOIP 2009 普及组第三题

RE的原因：

1.一开始是因为数组开小了

2.后面是因为没有对s[i]==1的时候continue;

对1分解质因数毫无意义，并且1不可能分裂

导致测试用例有s[]中有1的情况，变成了死循环所以 RE了

#include <bits/stdc++.h>
#define x first
#define y second
using namespace std;
typedef long long LL;
typedef pair<int,int> PII;
const int N=1e4+10;
const int M=45010;
LL s[N];
int n,m1,m2;

int primes[M],cnt;
bool st[M];

vector<PII> re;
vector<PII> t;

void init(int n)
{
	for(int i=2;i<=n;i++)
	{
		if(!st[i])
		{
			primes[cnt++]=i;
		}
		
		for(int j=0;primes[j]*i<=n;j++)
		{
			st[primes[j]*i]=true;
			if(i%primes[j]==0) break;
		}
	}

}

void div(vector<PII> &re,LL a,int f)
{
	for(int i=0;i<cnt;i++)
	{
		int p=primes[i];
		int co=0;
		if(a%p==0)
		{

		while(a%p==0)
		{
			a/=p;
			co++;
		}
		if(f==1)
		re.push_back({p,co*m2});
		else re.push_back({p,co});
    	}
	}
	
	if(a>1)
	{
		if(f==1)
		re.push_back({a,1*m2});
		else re.push_back({a,1});
	}
}

int main()
{
	init(M-1);
	
	cin>>n>>m1>>m2;
	for(int i=1;i<=n;i++)
	cin>>s[i];  //输入细胞个数 
	
	if(m1==1)
	{
		puts("0");
		return 0;
	}
	
	div(re,m1,1);
	
	/*
	for(int i=0;i<re.size();i++)
	cout<<re[i].x<<" ";
	cout<<endl;
	*/
	
	int minv=0x3f3f3f3f;
	bool f=false;
	for(int i=1;i<=n;i++)
	{
	    t.clear();
	    if(s[i]==1) continue; //当某个si是1的时候，这个细胞永远分裂不了想要的结果 
		div(t,s[i],0);
		bool flag=true;

		int maxv=0;
		for(int j=0;j<re.size();j++)
		{
			int x=re[j].x,co=re[j].y;
			//cout<<"在m中，x是 "<<x<<" 它的数量是 "<<co<<endl; 
			int l=0,r=t.size()-1;
			while(l<r)
			{
				int mid=(l+r)>>1;
				if(t[mid].x>=x) r=mid;
				else l=mid+1;
			}
			
			if(t[r].x!=x)
			{
				flag=false;
				break;
			}
			else
			{
				//cout<<"在s[]也就是 "<<s[i]<<"  x的数量是 "<<t[r].y<<endl;
				if(co>t[r].y)
				{
					int sec=0;
					int sum=t[r].y;
					if(sum<co)
					{
					  sec=(int)ceil((double)co/(double)sum);
					  
					}
					maxv=max(maxv,sec);
				}
			}
		}
		if(flag) 
		{
		minv=min(maxv,minv);
		f=true;
    	}
		
	}
	
	if(!f) puts("-1");
	else cout<<minv<<endl;
	
	return 0;
	
}