BSOJ:3801--隐藏指令 递推

该博客介绍了3801号问题——【四校联考1】隐藏指令的解决策略。内容涉及递推方法,通过递推公式f[i][j]=sigma(f[i-1][k]C(2*n-2*k,2(j-k))C(2(j-k),j-k))来计算在限定步数内使状态变为幸福的指令数量,其中C表示组合数。预处理杨辉三角并在线质因数分解可能导致超时。" 132232103,19687618,R语言实战:DFFITS度量在回归模型中的应用,"['r语言', '统计分析', '回归模型']

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

3801 – 【四校联考1】隐藏指令
Description
这里写图片描述

Input
输入仅一行,两个用空格分开的非负整数d,N。

Output
输出仅一行,仅一个整数表示能够变得幸福的指令数mod 1 000 000 007。

Sample IO
这里写图片描述
这里写图片描述

数据范围
这里写图片描述

递推:f[i][j]表示前i维空间走2*j步的方案数。
方程:f[i][j]=sigma(f[i-1][k]C(2*n-2*k,2(j-k))C(2(j-k),j-k)); (0<=k<=j)
简单解释一下这个方程:
已经走了2*k步,还剩2*n-2*k步,在这么多步中选出2*(j-k)步,即C(2*n-2*k,2*(j-k));
由题目得知走的过程是对称的,即在2*(j-k)中选出(j-k)个,也即C(2*(j-k),j-k);

这样递推即可,C可以预处理杨辉三角,在线质因数分解亲测超时。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#define mod 1000000007ll
#define ll long long
using namespace std;
inline ll read()
{
    ll bj=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')bj=-1;
        ch=getchar();
    }
    ll ret=0;
    while(ch>='0'&&ch<='9')ret=ret*10+ch-'0',ch=getchar();
    return ret*bj;
}
ll f[505][505]={0},c[505][505]={0};
ll d,n;
void getc()
{
    c[0][0]=1;
    for(ll i=1;i<=500;i++)c[i][0]=c[i][i]=1;
    for(ll i=1;i<=500;i++)
    {
        for(ll j=1;j<i;j++)c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
    }
}
void work()
{
    d=read();n=read();
    f[0][0]=1;
    for(ll i=1;i<=d;i++)f[i][0]=1;
    for(ll i=1;i<=d;i++)
    {
        for(ll j=1;j<=n;j++)
        {
            for(ll k=0;k<=j;k++)
            {
                f[i][j]=(f[i][j]+f[i-1][k]*c[2*(j-k)][(j-k)]%mod*c[2*n-k*2][2*(j-k)]%mod)%mod;
            }
        }
    }
    cout<<f[d][n];
}
int main()
{
    getc();
    work();
return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值