HDU 4028

题目：http://acm.hdu.edu.cn/showproblem.php?pid=4028

线性离散化DP。。。表示不会。。如果直接用数组存放会爆掉内存的。。所以用map。

DP[i][j]是以第i个指针为结束的最小公倍数j的方案数。

typedef    map<LL,LL>mp;     第一个表示第i个指针为结束的最小公倍数j，第二个为以第i个指针为结束的最小公倍数j的方案数

下面是AC代码：

#include<iostream>
#include<map>
using namespace std;

#define LL __int64 
LL m;
int n;
typedef map<LL,LL> mp;
mp dp[45];                        //离散化DP，DP[i][j]是以第i个指针为结束的最小公倍数j的方案数。

LL gcd(LL a,LL b){
	return b==0?a:gcd(b,a%b);
}

LL lcm(LL a, LL b){
	return a*b/gcd(a,b);
}
void DP(){
	int i;
	dp[1][1]=1;
	
	for(i=2;i<=40;i++){
		dp[i]=dp[i-1];
		dp[i][i]+=1;
		
		for(mp::iterator it=dp[i-1].begin();it!=dp[i-1].end();it++){
			LL t=lcm(it->first,i);
			
			dp[i][t]+=it->second;
		}
		
	}
	
	
}

int main()
{
	int T,t;
	DP();
	scanf("%d",&T);
	for(t=1;t<=T;t++)
	{
		scanf("%d%I64d",&n,&m);
		printf("Case #%d: ",t);
		LL ans=0;
		for(mp::iterator i=dp[n].begin();i!=dp[n].end();i++)
			if(i->first>=m)ans+=i->second;
			printf("%I64d\n",ans);
	}
	return 0;
}