URAL 1036 Lucky Tickets (高精度，数位dp)

最新推荐文章于 2018-03-04 21:11:45 发布

My_ACM_Dream

最新推荐文章于 2018-03-04 21:11:45 发布

阅读量779

点赞数

CC 4.0 BY-SA版权

分类专栏：动态规划—数位dp 动态规划

本文链接：https://blog.youkuaiyun.com/My_ACM_Dream/article/details/44609921

动态规划同时被 2 个专栏收录

246 篇文章

订阅专栏

动态规划—数位dp

30 篇文章

订阅专栏

博客探讨了一种求解特定数值条件的计数问题，即找到所有前半部分数字之和等于后半部分数字之和，且该和等于给定值S的2*N位数字的个数。通过使用高精度存储的动态规划方法，可以解决这个问题。当总和为奇数时，由于题目条件无法满足，答案为0。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题意：

求2*N位的数，前N位数和等于后N位数和并且这个和等于S对应的数的个数。

题解：

高精度存dp，dp[i][j]表示前i位和为j的个数，状态转移很容易。注意和为奇数时无解所以结果是0.

#include<iostream>
#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<vector>
#include<queue>
#include<map>
#include<set>
using namespace std;
#define B(x) (1<<(x))
typedef long long ll;
void cmax(int& a,int b){ if(b>a)a=b; }
void cmin(int& a,int b){ if(b<a)a=b; }
const int oo=0x3f3f3f3f;
const int MOD=1000000007;
const int maxn=105;
struct BigInt
{
    const static int mod = 10;
    int a[maxn],len;
    BigInt()
    {
        memset(a,0,sizeof(a));
        len = 1;
    }
    BigInt(int v)
    {
        memset(a,0,sizeof(a));
        len = 0;
        do
        {
            a[len++] = v%mod;
            v /= mod;
        }while(v);
    }
    BigInt operator +(const BigInt &b)const
    {
        BigInt res;
        res.len = max(len,b.len);
        for(int i = 0;i <= res.len;i++)
            res.a[i] = 0;
        for(int i = 0;i < res.len;i++)
        {
            res.a[i] += ((i < len)?a[i]:0)+((i < b.len)?b.a[i]:0);
            res.a[i+1] += res.a[i]/mod;
            res.a[i] %= mod;
        }
        if(res.a[res.len] > 0)res.len++;
        return res;
    }

    BigInt operator *(const BigInt &b)const
    {
        BigInt res;
        for(int i = 0; i < len;i++)
        {
            int up = 0;
            for(int j = 0;j < b.len;j++)
            {
                int temp = a[i]*b.a[j] + res.a[i+j] + up;
                res.a[i+j] = temp%mod;
                up = temp/mod;
            }
            if(up != 0)
                res.a[i + b.len] = up;
        }
        res.len = len + b.len;
        while(res.a[res.len - 1] == 0 &&res.len > 1)res.len--;
        return res;
    }

    void output()
    {
        printf("%d",a[len-1]);
        for(int i = len-2;i >=0 ;i--)
            printf("%d",a[i]);
        printf("\n");
    }

    void init(){
        memset(a,0,sizeof a);
        len=1;
    }

}dp[55][505];

void CL(){
    for(int i=0;i<55;i++){
        for(int j=0;j<505;j++){
            dp[i][j].init();
        }
    }
}


int main(){
    //freopen("E:\\read.txt","r",stdin);
    int n,s;
    while(scanf("%d %d",&n,&s)!=EOF){
        if(s&1){
            printf("0\n");
            continue;
        }
        CL();
        s/=2;
        dp[0][0]=1;
        for(int i=1;i<=n;i++){
            for(int j=0;j<=s;j++){
                for(int k=0;k<=9;k++){
                    if(j-k>=0)
                        dp[i][j]=dp[i][j]+dp[i-1][j-k];
                }
            }
        }
        BigInt ans=dp[n][s]*dp[n][s];
        ans.output();
    }
	return 0;
}