[HGOI]CSP2019训练模拟赛AM-题解01

佐手、 提交于 2019-12-02 19:47:48

C 小G坐电梯

题目描述

小G来到了著名的某大厦。大厦一共有 n 层,初始的时候小 G 在第 A 层。
小G特别想去 B 层小 M 的办公室看一看,然而因为安保原因,B 层已经被封锁无法进入。
但是小G既然来了,就想在大厦里面逛一逛。大厦里面有一部电梯,小G决定坐 k 次电梯。
因为小G比较无聊,他给自己设定了这样一个规矩:假如当前他在 x 层,则他要去的下一个楼层 y 和 x 的楼层差必须要小于 x 和 B 的楼层差,即 \(|x−y| < |x−B|\)
每到达一个楼层,小 G 都要记录下来其楼层号。
当小G转完一圈后,他也记录下了 \(k + 1\) 个楼层号(可能有重复)。
小G现在 想知道,按照他定下的规矩,一共有多少种可能的楼层号序列?

题解

定义f[i][j][k]为当前走了i步,距离终点为j的,方向为k(0向下,1向上)。
暴力DP+前缀和,滚动数组压掉一维。

代码

#include <cstdio>
#include <cmath>

using namespace std;

#define max(a,b) ((a>b)?a:b)

const int MOD = 1e9 + 7;
           //k     dlt  lft 0:dwn, 1:up
long long f[10005][2];
long long pre[20005][2];

int mx_dlt, dlt;

inline long long pls(long long a, long long b){
    return ((a + b >= MOD) ? (a + b - MOD) : (a + b));
}

inline void getPre(){
    pre[0][0] = pre[0][1] = 0;
    for (int j = 1; j <= mx_dlt; ++j){
        pre[j][0] = pls(pre[j - 1][0], f[j][0]);
        pre[j][1] = pls(pre[j - 1][1], f[j][1]);
    }
    for (int j = mx_dlt + 1; j <= mx_dlt * 2; ++j){ 
        pre[j][0] = pre[j - 1][0];
        pre[j][1] = pre[j - 1][1];
    }
}

int main(){
    freopen("lift.in", "r", stdin);
    freopen("lift.out", "w", stdout);
    int n, A, B, k; scanf("%d %d %d %d", &n, &A, &B, &k);
    mx_dlt = max(B * 2, abs(n - B) * 2); dlt = abs(A - B);
    f[dlt][0] = (A < B), f[dlt][1] = (A > B);
    getPre();
    for (int i = 1; i <= k; ++i){
        for (int j = 1; j <= mx_dlt; ++j){
            int tmp0 = f[j][0], tmp1 = f[j][1];
            if (B - j > 0)
                f[j][0] = ((pre[mx_dlt][0] - pre[j / 2][0] - tmp0) % MOD + MOD) % MOD;
            if (B + j <= n)
                f[j][1] = ((pre[mx_dlt][1] - pre[j / 2][1] - tmp1) % MOD + MOD) % MOD;
        }
        getPre();
    }
    long long ans = (pre[mx_dlt][0] + pre[mx_dlt][1] + MOD * 2) % MOD;
    printf("%lld", ans);
    return 0;
}
标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!