my ability is poor!

原创 2012年03月23日 20:08:21
#include <iostream>
#include <cmath>
#include <queue>

using namespace std;
#define  M 10000

int prime[M];
bool is[M];
int factor[20];

struct E
{
	int a,b;
	E(int a=0,int b=0){this->a=a;this->b=b;};
};
void creatprime()
{
	memset(is,0,sizeof(is));
	int k=0;
	prime[k++]=2;
	for(int i=4;i<M;i+=2)
		is[i]= true;
	for(int i=3;i<M;i+=2)
	{
		if(!is[i])
		{
			prime[k++] = i;
			for(int j=i+i;j<M;j+=i)
				is[j] = true;
		}
	}
	return ;
}

void creatfactor(int n)
{
	memset(factor,0,sizeof(factor));
	int k=0;int m=n;
	for(int i=0;prime[i]*prime[i]<=m && n!=1;i++)
		if(n%prime[i]==0)
		{
			factor[k++] = prime[i];
			while(n%prime[i]==0)
				n/=prime[i];
		}
		if(n!=1) factor[k++] = n;
		return ;
}
int howmany(int n)
{
	queue<E> G;
	for(int i=0;i<20&&factor[i]!=0;i++)
		G.push(E(factor[i],factor[i]));
	int sum = 0;
	while(!G.empty())
	{
		E ele = G.front();G.pop();
		sum += n/ele.a;
		for(int i=0;i<20&&factor[i]!=0;i++)
		{   
			if(factor[i]>ele.b)
			{
				int temp = ele.a;
				temp *= -factor[i];
				if(n/abs(temp)!=0)
				{
					G.push(E(temp,factor[i]));
				}
			}
		}
	}
	return sum;
	
}
bool hasfactor(int n)
{
	for(int i=0;i<20&&factor[i]!=0;i++)
		if(n%factor[i]==0)return true;
	return false;
}
int main()
{ 
	int k,m;
	
	creatprime();

	while(scanf("%d%d",&k,&m)!=EOF)
	{
		if(m==1){puts("1");continue;}
		if(k==1){printf("%d\n",m);continue;}
		creatfactor(k);

		int temp = howmany(m);
		temp = m-temp;
		for(int i=m+1;;i++)
		{
			if(!hasfactor(i))temp++;
			if(temp==m) {printf("%d\n",i);break;}
		}
	}
    return 0;
}

版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

this is my first article

  • 2011-02-24 11:03
  • 1.31MB
  • 下载

is my recording

  • 2010-11-24 06:04
  • 121KB
  • 下载

#ASM 翻译系列第二十六弹:ASM 高级知识 Where is my data

原文: Where is my data 作者: Bane Radulovic 译者:邱大龙,沃趣科技数据库技术专家,主要参与公司产品实施、测试、维护以及优化。对SQL优化、数据迁移、备份...

this is my SMSC document2

  • 2010-01-25 14:34
  • 102KB
  • 下载

this is my SMSC document4

  • 2010-01-25 14:35
  • 146KB
  • 下载

udev is renaming my ethernet devices

一台 HP 的 centos 6 机器,原来有 onboard 的 Broadcom 网卡在 eth2/eth3,另外还有一张另外加上的 intel 82576 网卡 eth0/eth1。最近更换了那...

this is my SMSC document3

  • 2010-01-25 14:35
  • 1.13MB
  • 下载

The ultimate SO_LINGER page, or: why is my tcp not reliable

转自:http://blog.csdn.net/cpp_chen/article/details/29864509 This post is about an obscure corner of T...

Who Is Using My WIFI

  • 2014-03-19 19:50
  • 262KB
  • 下载
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)