最小公倍数(nyoj 517)

最小公倍数

时间限制： 1000 ms  |  内存限制： 65535 KB

难度： 3

描述

为什么1小时有60分钟，而不是100分钟呢？这是历史上的习惯导致。

但也并非纯粹的偶然：60是个优秀的数字，它的因子比较多。

事实上，它是1至6的每个数字的倍数。即1,2,3,4,5,6都是可以除尽60。

我们希望寻找到能除尽1至n的的每个数字的最小整数m.

输入

多组测试数据（少于500组）。
每行只有一个数n(1<=n<=100).

输出

输出相应的m。

样例输入

2
3
4

样例输出

2
6
12

这一题运用  最小公倍数就是所有质数的相应幂的乘积。

例如n = 10，小于10的质数有2，3，5，7，对应的最大幂是：3，2，1，1（质数的对应幂不能大于n），则最小公倍数是：2^3x3^2x5^1x7^1  =  2520。

#include <stdio.h>
#include <string.h>

int pri[100] = {1, 1, 0};
int ans[101][45];//ans[i]存储1-i的最小公倍数

void prime()
{//筛选法
	int i, j;
	for(i = 2;i <= 100; i++)
	{
		if(!pri[i])
		{
			for(j = i + i; j <= 100; j += i)
                   pri[j] = 1;
		}
	}
}

void answer()
{
	int i, j, k, l;
	for(i = 3; i <= 100; i++)
	{
		ans[i][0] = 1;//因为后面都是相乘，所以要置为1
		for(j = 2; j <= i; j++)
		{
			if(!pri[j])//是素数
			{
				k = j;
				while(k <= i)
				{//质数的对应幂不能大于n
					ans[i][0] *= j;
					for(l = 1; l < 45; l++)
					{//每乘一个数进行一次进位处理
						ans[i][l] = ans[i][l] * j + ans[i][l - 1] / 10;
						ans[i][l - 1] %= 10;
					}
					k *= j;
				}
			}
		}
	}
}

int main (void)
{
	int n, i;
	prime();
	answer();
	while(scanf("%d", &n) != EOF)
	{
		if(n == 1 || n == 2)
			printf("%d\n", n);
		else
		{
			for(i = 44; i >= 0; i --)
				if(ans[n][i] != 0)//找到首个不是0的数
					break;
			for(; i >= 0; i--)
				printf("%d", ans[n][i]);
			printf("\n");
		}
	}
	return 0;
}