LOJ6165

• 题目
  LOJ6165
• 分析
  运用$gcd$的性质：
  $a=p_1^{a_1}\ast p_2^{a_2}...p_k^{a_k}$
  $b=p_1^{b_1}\ast p_2^{b_2}...p_k^{b_k}$
  $lcm(a,b)=p_1^{max(a_1,b_1)}\ast p_2^{max(a_2,b_2)}...p_k^{max(a_k,b_k)}$
  推广到$n$个数
  本题就是对$1\sim n$求最小公倍数。那么我们考虑对于每个$p$最大会有多少个，分析可知对于$n$越大，$p$的$c$就会越大，所以就可以直接考虑等于$n$的情况，看小于$n$的质数，对于每一个质数分别最多会有多少个，最后把这些数累加起来就可以了。
• 代码

/*
  独立思考
  一个题不会做，收获5%，写了代码10%，提交对了30%，总结吃透了这个题才是100%.
*/
#include
using namespace std;
template <typename T>
void read(T &x)
{
	x = 0;
	char c = getchar();
	int sgn = 1;
	while (c < '0' || c > '9') {if (c == '-')sgn = -1; c = getchar();}
	while (c >= '0' && c <= '9')x = x * 10 + c - '0', c = getchar();
	x *= sgn;
}
template <typename T>
void out(T x)
{
	if (x < 0) {putchar('-'); x = -x;}
	if (x >= 10)out(x / 10);
	putchar(x % 10 + '0');
}
typedef long long ll;
typedef unsigned long long ull;
ll gcd(ll a, ll b) { return b ? gcd(b, a % b) : a;}
const int N = 1e8 + 5;
const ll mod = 100000007;
bitset<N> v;
int prime[7000005];
ll ans = 1;
int tot = 0;
int n;
void sieve()
{
	v[0] = v[1] = 1;
	for (register int i = 2; i <= n; i++)
	{
		if (!v[i]) {
			prime[++tot] = i;
			for (register ll k = i; k <= n; k *= i)//求对于每一个质数p最多有多少个。
			{
				ans = ans * 1ll * i % mod;
			}
		}
		for (register int j = 1; j <= tot && i * prime[j] <= N; j++)
		{
			v[i * prime[j]] = 1;
			if (i % prime[j] == 0)
				break;
		}
	}
}
int main ()
{
	cin >> n;
	sieve();
	cout << ans << endl;
	return 0 ;
}

• 方法
  $gcd,lcm$性质