hdu--4028--dp

这个dp我没做出来啊...其实不难..主要题意没理解好 fuck.

给你1-N这N个数  一共2^N-1个子集 每个子集的LCM值>=M的情况数有多少种

我也是醉了 这么个题目 给我套他那个题面 硬是没看懂 他在问什么 还是 英语太渣了

然后就是个 状态转移方程的考虑了

map<LL,LL>dp[size] 表示  前size个数 构成的lcm值为 it->first的情况为 it->second种

dp[ i ] = dp[ i-1 ]  //不添加第 I 个元素

dp[ i ][ i ] ++//第I个元素自身构成的集合

for it - >   begin() -  end()

dp[ i][ it->fist ] += it->second //第I个元素与前面的 I-1元素构成的集合组合出的情况

//给你1-N 这N个数  问有多少个子集 该集合的LCM >= M

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

const int size = 40;
typedef long long LL;
map<LL,LL>dp[size+5];//前 i 个数组合出的lcm值
map<LL,LL>::iterator it;

LL gcd( LL x , LL y )
{
    return x % y == 0 ? y : gcd( y , x%y );
}

void init( int n )
{
    dp[1][1] = 1;
    for( int i = 1 ; i<=n ; i++ )
    {
        dp[i] = dp[i-1];
        dp[i][i] ++;
        for( it=dp[i-1].begin() ; it!=dp[i-1].end(); it++ )
        {
            dp[i][ it->first*i/gcd(i,it->first) ] += it->second;
        }
    }
}

int main()
{
    cin.sync_with_stdio(false);
    int n , t;
    LL m , ans;
    init( size );
    cin >> t;
    for( int k = 1 ; k<=t ; k++ )
    {
        cin >> n >> m;
        ans = 0;
        for( it = dp[n].begin() ; it!=dp[n].end() ; it++ )
        {
            if( it->first>=m )
            {
                ans += it->second;
            }
        }
        cout << "Case #" << k << ": " << ans << endl;
    }
    return 0;
}