hdu 4028 2011上海赛区网络赛H dp+map离散
一开始用搜索直接超时，看题解会的
1 #include<iostream> 2 #include<cstdio> 3 #include<map> 4 #include<cstring> 5 #include<cmath> 6 #include<vector> 7 #include<queue> 8 #include<algorithm> 9 #include<set> 10 #define inf 110000 11 #define M 10005 12 #define N 10005 13 #define Min(a,b) ((a)<(b)?(a):(b)) 14 #define Max(a,b) ((a)>(b)?(a):(b)) 15 #define pb(a) push_back(a) 16 #define mem(a,b) memset(a,b,sizeof(a)) 17 #define eps 1e-9 18 #define zero(a) fabs(a)<eps 19 #define LL long long 20 #define MOD 1000000007 21 using namespace std; 22 map<LL,LL>dp[45]; 23 map<LL,LL>::iterator it; 24 LL gcd(LL a,LL b){ 25 return b==0?a:gcd(b,a%b); 26 } 27 LL lcm(LL a,LL b){ 28 return a/gcd(a,b)*b; 29 } 30 void DP(){ 31 dp[1][1]=1; 32 for(int i=2;i<=40;i++){ 33 dp[i]=dp[i-1]; //不取第i个数的所有情况，先复制过来 34 dp[i][i]++; //只取第i个数，不能落下 35 for(it=dp[i-1].begin();it!=dp[i-1].end();it++) 36 dp[i][lcm(i,it->first)]+=it->second; //然后考虑在前i-1个数的基础上加入第i个数 37 } 38 } 39 LL n,m; 40 int main(){ 41 DP(); 42 int t,cas=0; 43 scanf("%d",&t); 44 while(t--){ 45 scanf("%I64d%I64d",&n,&m); 46 LL ans=0; 47 //遍历一遍 48 for(it=dp[n].begin();it!=dp[n].end();it++) 49 if(it->first>=m) 50 ans+=it->second; 51 printf("Case #%d: %I64d\n",++cas,ans); 52 } 53 return 0; 54 }
