CodeChef - CODIE

%%%

参考自:yts1999’s blog

题目链接 : Codie Bird

sol:

  1. 考虑朴素的 d p dp dp方程

d p [ i ] [ j ] = d p [ i − 1 ] [ j − 1 ] + d p [ i − 1 ] [ j ] + d p [ i ] [ j + 1 ] dp[i][j] = dp[i-1][j-1]+dp[i-1][j]+dp[i][j+1] dp[i][j]=dp[i1][j1]+dp[i1][j]+dp[i][j+1]

递推矩阵大概长这样。
T = [ 1 1 0 0 0 0 … 0 0 1 1 1 0 0 0 … 0 0 0 1 1 1 0 0 … 0 0 … … … … … … … … … 0 0 0 0 0 0 … 1 1 ] T=\begin{bmatrix} 1 & 1 & 0 & 0 & 0 & 0 & \dots & 0 & 0 \\ 1 & 1 & 1 & 0 & 0 & 0 & \dots & 0 & 0 \\ 0 & 1 & 1 & 1 & 0 & 0 & \dots & 0 & 0 \\ \dots & \dots & \dots & \dots & \dots & \dots & \dots & \dots & \dots \\ 0 & 0 & 0 & 0 & 0 & 0 & \dots & 1 & 1\\ \end{bmatrix} T=11001110011000100000000000010001

有障碍的时候,将 [ 1 , a ] [1,a] [1,a] [ b , m ] [b,m] [b,m]列的矩阵元素全部置为0。将这个转移矩阵设为 S S S.
v i = [ d p [ i ] [ 1 ] d p [ i ] [ 2 ] ⋯ d p [ i ] [ k ] ] v_i=\begin{bmatrix} dp[i][1] & dp[i][2] & \cdots & dp[i][k]\end{bmatrix} vi=[dp[i][1]dp[i][2]dp[i][k]],则 v n = v 1 ( ∑ i = 1 n − 2 T i − 1 S T n − 2 − i ) T v_n=v_1(\sum\limits_{i=1}^{n-2}T^{i-1}ST^{n-2-i})T vn=v1(i=1n2Ti1STn2i)T
P n = ∑ i = 1 n T i − 1 S T n − i P_n = \sum\limits_{i=1}^{n} T^{i-1}ST^{n-i} Pn=i=1nTi1STni,则 v n = v 1 P n − 2 T v_n = v_1P_{n-2}T vn=v1Pn2T
实现上可以用单位矩阵代替 v 1 v_1 v1 v n v_n vn就是最终矩阵的第 k 2 \frac{k}{2} 2k行。

  1. P n P_n Pn可以用类似倍增的做法。

P 1 = S P_1 = S P1=S
若n为偶数, P n = P n 2 ⋅ T n 2 + T n 2 ⋅ P n 2 P_n = P_\frac{n}{2} \cdot T^\frac{n}{2} + T^\frac{n}{2} \cdot P_\frac{n}{2} Pn=P2nT2n+T2nP2n
若n为奇数, P n = P n − 1 ⋅ T + T n − 1 ⋅ S P_n = P_{n-1} \cdot T + T^{n-1} \cdot S Pn=Pn1T+Tn1S

实现上可以考虑 101 0 ( 2 ) = 1 0 ( 10 ) 1010_{(2)} = 10_{(10)} 1010(2)=10(10)这样的求解过程。

code:

#include<iostream>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>

using namespace std;
typedef long long ll;
const int maxn = 52;
const int mod = 1e9+7;
int n,m;

struct Mat{
    ll a[maxn][maxn];
    void init(int k){
        memset(a,0,sizeof(a));
        for(int i = 1;i<=k;i++) a[i][i] = 1;
    }

    void reset(){
        memset(a,0,sizeof(a));
    }

    Mat operator + (Mat b){
        Mat c;
        c.reset();
        for(int i = 1;i<=m;i++)
        for(int j = 1;j<=m;j++)
            c.a[i][j] = (a[i][j] + b.a[i][j]) % mod;
        return c;
    }

    Mat operator * (Mat b){
        Mat c; 
        for(int i = 1;i<=m;i++)
        for(int j = 1;j<=m;j++){
            c.a[i][j] = 0;
            for(int k = 1;k<=m;k++){
                c.a[i][j] += a[i][k] * b.a[k][j] % mod;
            }
            c.a[i][j] %= mod;
        }
        return c;
    }

    Mat pow(int b){
        Mat ret; ret.init(m);
        Mat a = (*this);
        while(b){
            if(b&1) ret = ret * a;
            a = a * a;
            b>>=1;
        }
        return ret;
    }
};

ll qpow(ll a,int b){
    ll ret = 1;
    while(b){
        if(b&1) ret = ret * a % mod;
        a = a*a % mod;
        b>>=1;
    }
    return ret;
}

int main(){
    int a,b;
    scanf("%d%d%d%d",&n,&m,&a,&b);
    n -= 2;
    Mat T,S;
    T.reset(); S.reset();
    for(int i = 1;i<=m;i++) 
        T.a[i][i] = T.a[i][i-1] = T.a[i][i+1] = 1;
    
    S = T;
    for(int i = 1;i<=m;i++){
        if(i<=a || i>=b){
            for(int j = 1;j<=m;j++){
                S.a[j][i] = 0;
            }
        }
    }
    
    int bit = 0;
    for(int i = n;i>0;i>>=1) bit++;
    Mat P = S; Mat Ts = T;
    for(int i = bit - 2;i>=0;i--){
        P = P * Ts + Ts * P;
        Ts = Ts * Ts;
        if((n>>i)&1){
            P = P * T + Ts * S;
            Ts = Ts * T;
        }
    }
    P = P * T;
    ll ans = P.a[m/2][m/2];
    ans = ans * qpow(n,mod-2) % mod;
    printf("%lld\n",ans);
    return 0;
}

1. 用户与身体信息管理模块 用户信息管理: 注册登录:支持手机号 / 邮箱注册,密码加密存储,提供第三方快捷登录(模拟) 个人资料:记录基本信息(姓名、年龄、性别、身高、体重、职业) 健康目标:用户设置目标(如 “减重 5kg”“增肌”“维持健康”)及期望周期 身体状态跟踪: 体重记录:定期录入体重数据,生成体重变化曲线(折线图) 身体指标:记录 BMI(自动计算)、体脂率(可选)、基础代谢率(根据身高体重估算) 健康状况:用户可填写特殊情况(如糖尿病、过敏食物、素食偏好),系统据此调整推荐 2. 膳食记录与食物数据库模块 食物数据库: 基础信息:包含常见食物(如米饭、鸡蛋、牛肉)的名称、类别(主食 / 肉类 / 蔬菜等)、每份重量 营养成分:记录每 100g 食物的热量(kcal)、蛋白质、脂肪、碳水化合物、维生素、矿物质含量 数据库维护:管理员可添加新食物、更新营养数据,支持按名称 / 类别检索 膳食记录功能: 快速记录:用户选择食物、输入食用量(克 / 份),系统自动计算摄入的营养成分 餐次分类:按早餐 / 午餐 / 晚餐 / 加餐分类记录,支持上传餐食照片(可选) 批量操作:提供常见套餐模板(如 “三明治 + 牛奶”),一键添加到记录 历史记录:按日期查看过往膳食记录,支持编辑 / 删除错误记录 3. 营养分析模块 每日营养摄入分析: 核心指标计算:统计当日摄入的总热量、蛋白质 / 脂肪 / 碳水化合物占比(按每日推荐量对比) 微量营养素分析:检查维生素(如维生素 C、钙、铁)的摄入是否达标 平衡评估:生成 “营养平衡度” 评分(0-100 分),指出摄入过剩或不足的营养素 趋势分析: 周 / 月营养趋势:用折线图展示近 7 天 / 30 天的热量、三大营养素摄入变化 对比分析:将实际摄入与推荐量对比(如 “蛋白质摄入仅达到推荐量的 70%”) 目标达成率:针对健
1. 用户管理模块 用户注册与认证: 注册:用户填写身份信息(姓名、身份证号、手机号)、设置登录密码(需符合复杂度要求),系统生成唯一客户号 登录:支持账号(客户号 / 手机号)+ 密码登录,提供验证码登录、忘记密码(通过手机验证码重置)功能 身份验证:注册后需完成实名认证(模拟上传身份证照片,系统标记认证状态) 个人信息管理: 基本信息:查看 / 修改联系地址、紧急联系人、邮箱等非核心信息(身份证号等关键信息不可修改) 安全设置:修改登录密码、设置交易密码(用于转账等敏感操作)、开启 / 关闭登录提醒 权限控制:普通用户仅能操作本人账户;管理员可管理用户信息、查看系统统计数据 2. 账户与资金管理模块 账户管理: 账户创建:用户可开通储蓄卡账户(默认 1 个主账户,支持最多 3 个子账户,如 “日常消费账户”“储蓄账户”) 账户查询:查看各账户余额、开户日期、状态(正常 / 冻结)、交易限额 账户操作:挂失 / 解挂账户、申请注销账户(需余额为 0) 资金操作: 转账汇款:支持同行转账(输入对方账户号 / 手机号),需验证交易密码,可添加常用收款人 存款 / 取款:模拟存款(输入金额增加余额)、取款(输入金额减少余额,需不超过可用余额) 交易记录:按时间、类型(转入 / 转出 / 存款 / 取款)查询明细,显示交易时间、金额、对方账户(脱敏显示)、交易状态 3. 账单与支付模块 账单管理: 月度账单:自动生成每月收支明细,统计总收入、总支出、余额变动 账单查询:按月份、交易类型筛选账单,支持导出为 Excel 格式 还款提醒:若有贷款(简化版可模拟),系统在还款日 3 天前发送提醒 快捷支付: 绑定支付方式:添加银行卡(系统内账户)作为支付渠道 模拟消费:支持输入商户名称和金额,完成支付(从账户余额扣减) 支付记录:保存所有消费记录,包含商户、时间、金额、支付状态 4.
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值