【HDU 4028】The time of a day

1.题目链接。题目大意：n个指针，他们的周期分别是[1,n]问这n个数中任选k个,(0<=k<=n)最小公倍数大于等于m的方案书有多少种。

2.显然DP。简单推导一下DP转移：DP[i][j]代表前i个数中，lcm为j的方案数。(注意这里我们是没有考虑选区几个数的，因为没有必要，我们考虑的对象只是lcm，而不是几个数组成了这个lcm).然后就成了一个计数类的DP。

DP[i][j]+=DP[i-1][k][lcm(i,k)==j].这个转移应该是显然的，因为DP[i][j]从他上一个转移过来的时候，只要是dp[i-1][k]并且lcm(i,k)等于j即可。但是我们这里式没有办法这样写的，因为j很大，数组显然是开不下的。我们可以用map代替数组，然后就可以完成了。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#pragma warning(disable:4996)
map<ll, ll>dp[50];
ll n, m;
ll gcd(ll a, ll b)
{
	return b == 0 ? a : gcd(b, a%b);
}
ll lcm(ll a, ll b)
{
	return a / gcd(a, b)*b;
}
void init()
{
	dp[1][1] = 1;
	for (int i = 2; i <= 40; i++)
	{
		dp[i] = dp[i - 1];//这里表示不选i
		dp[i][i]++;//这里表示只选i
		for (auto iter = dp[i - 1].begin(); iter != dp[i - 1].end(); iter++)
		{
			dp[i][lcm(i, iter->first)]+= iter->second;//这里就是我们的DP转移式了。
		}
	}
}
int main()
{
	init();
	int t, cas=1;
	scanf("%d", &t);
	while (t--)
	{
		scanf("%lld%lld", &n, &m);
		ll ans = 0;
		for (auto iter = dp[n].begin(); iter != dp[n].end(); iter++)
		{
			if (iter->first >= m)
				ans+=iter->second;
		}
		printf("Case #%d: %lld\n", cas++, ans);
	}
}