hdu 2082 找单词

该博客探讨了一种使用组合数学的母函数和动态规划(DP)方法来解决HDU 2082编程题目——找单词。文章提供了一个母函数的模板,并展示了如何通过二维DP实现状态转移方程,以计算使用特定字母数量形成特定价值的单词种类数。最终,博主给出了获得AC(正确答案)的代码片段。

这道题,可以用组合数学的母函数,也可以用dp背包。

先来一个母函数模板:

!!!转自Tanky Woo !!!

#include <iostream>  
using namespace std;  
// Author: Tanky Woo    
// www.wutianqi.com  
const int _max = 10001;   
// c1是保存各项质量砝码可以组合的数目  
// c2是中间量,保存每一次的情况  
int c1[_max], c2[_max];     
int main()  
{   //int n,i,j,k;  
    int nNum;   //   
    int i, j, k;  
  
    while(cin >> nNum)  
    {  
        for(i=0; i<=nNum; ++i)   // ---- ①  
        {  
            c1[i] = 1;  
            c2[i] = 0;  
        }  
        for(i=2; i<=nNum; ++i)   // ----- ②  
        {  
  
            for(j=0; j<=nNum; ++j)   // ----- ③  
                for(k=0; k+j<=nNum; k+=i)  // ---- ④  
                {  
                    c2[j+k] += c1[j];  
                }  
                for(j=0; j<=nNum; ++j)     // ---- ⑤  
                {  
                    c1[j] = c2[j];  
                    c2[j] = 0;  
                }  
        }  
        cout << c1[nNum] << endl;  
    }  
    return 0;  
}  
①  、首先对c1初始化,由第一个表达式(1+x+x^2+..x^n)初始化,把质量从0到n的所有砝码都初始化为1.
②  、 i从2到n遍历,这里i就是指第i个表达式,上面给出的第二种母函数关系式里,每一个括号括起来的就是一个表达式。
③、j 从0到n遍历,这里j就是(前面i個表达式累乘的表达式)里第j个变量,(这里感谢一下seagg朋友给我指出的错误,大家可以看下留言处的讨论)。
如(1+x)(1+x^2)(1+x^3),j先指示的是1和x的系数,i=2执行完之后变为(1+x+x^2+x^3)(1+x^3),
这时候j应该指示的是合并后的第一个括号的四个变量的系数。
④ 、 k表示的是第j个指数,所以k每次增i(因为第i个表达式的增量是i)。
⑤  、把c2的值赋给c1,而把c2初始化为0,因为c2每次是从一个表达式中开始的。

AC代码:

#include<bits/stdc++.h>
#define N 55
#define LL long long
using namespace std;

int a[33];
LL c1[N],c2[N];
int main()
{
    int T,i;
    scanf("%d",&T);
    while(T--)
    {
        memset(c1,0,sizeof(c1));
        for(i=1;i<=26;i++)
            scanf("%d",&a[i]);
        for(i=1;i<=26;i++)
        {
            if(a[i])
            {
                for(int j=0;j<=a[i];j++)
                {
                    if(j*i>50) break;
                    c1[j*i]=1;
                }
                break;
            }
        }

        for(i=i+1;i<=26;i++)
        {
            for(int j=0;j<=50;j++)
                if(c1[j])
                    for(int k=0;k<=a[i];k++)
                    {
                        if(j+i*k>50) break;
                        c2[j+k*i]+=c1[j];
                    }
            for(int j=0;j<=50;j++)
            {
                c1[j]=c2[j];
                c2[j]=0;
            }

        }
        LL sum=0;
        for(i=1;i<=50;i++)
           sum+=c1[i];
        printf("%lld\n",sum);
    }
    return 0;
}

接下来用dp 

二维DP    

dp[i][j]   分别表示使用i个字母 价值为j时的种类有多少。

 首先  dp[i][j]=dp[i-1][j];

   然后 状态转移方程  dp[i][j]=dp[i-1][j-k*i];

  最后sum=dp[26][1->50];

  即为所求.

ACd代码:

#include<bits/stdc++.h>
#define LL long long
using namespace std;
LL dp[30][60];
int a[30];
int main()
{
    int T;
    LL sum;
    scanf("%d",&T);
    while(T--)
    {
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=26;i++)
            scanf("%d",&a[i]);
         for(int i=0; i<=26; i++)
            dp[i][0]=1;
        for(int i=1;i<=26;i++)
        {
            for(int j=0;j<=50;j++)
            {
                dp[i][j]=dp[i-1][j];
                for(int k=1;k<=a[i];k++)
                {
                    if(j-i*k<0) break;
                    dp[i][j]+=dp[i-1][j-i*k];
                }
            }
        }
        sum=0;
        for(int i=1;i<=50;i++)
            sum+=dp[26][i];
        printf("%lld\n",sum);
    }
}



评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值