hdu4028 离散DP

最新推荐文章于 2022-07-15 07:41:41 发布

ls_0222

最新推荐文章于 2022-07-15 07:41:41 发布

阅读量515

点赞数

分类专栏： ACM 文章标签： iterator pair system

本文链接：https://blog.csdn.net/ls_0222/article/details/6784125

版权

ACM 专栏收录该内容

9 篇文章 0 订阅

订阅专栏

本文介绍了一个利用计数动态规划方法解决特定数学问题的实例：从40个数中选择若干个数，使得这些数的最小公倍数大于等于给定值m。文章通过C++代码详细展示了如何使用map来存储中间状态及其计数，并最终求得所有符合条件的选择方案的数量。

摘要由CSDN通过智能技术生成

以后一定要每天都写报告

题意：40个数，选出几个，使其lcm大于等于m,计数DP, map

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <map>
#define inf (1 << 30)
const double eps = 1e-12;
const double pi = acos(-1.0);
using namespace std;

int T;
__int64 n, m;
map<__int64, __int64> dp[41];
map<__int64, __int64> :: iterator iter;

__int64 gcd(__int64 a, __int64 b)
{
    return a % b ? gcd(b, a % b) : b;
}

__int64 lcm(__int64 a, __int64 b)
{
    return a * b / gcd(a, b);
}

int solve()
{
    int i;
    __int64 first, second, t1;
    //dp[1].insert(make_pair(1, 1));
    dp[1][1] = 1;
    for (i = 2; i <= 40; i++)
    {
        dp[i] = dp[i - 1];
        dp[i][i]++;
        for (iter = dp[i - 1].begin(); iter != dp[i - 1].end(); iter++)
        {
            first = iter -> first, second = iter -> second;
            t1 = lcm(first, i);
            dp[i][t1] += second;
        }
    }
    return 1;
}

int main()
{
    int i, j, h;
    __int64 ans, first, second;
    solve();
    scanf("%d", &T);
    for (h = 1; h <= T; h++)
    {
        scanf("%I64d%I64d", &n, &m);
        printf("Case #%d: ", h);
        ans = 0;
        for (iter = dp[n].begin(); iter != dp[n].end(); iter++)
        {
            first = iter -> first, second = iter -> second;
            if (first >= m) ans += second;
        }
        printf("%I64d\n", ans);
    }
    system("pause");
    return 0;
}