QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#396667#1882. DrunkardsdeepthoughtRE 0ms0kbC++233.0kb2024-04-22 23:28:272024-04-22 23:28:28

Judging History

你现在查看的是最新测评结果

  • [2024-04-22 23:28:28]
  • 评测
  • 测评结果:RE
  • 用时:0ms
  • 内存:0kb
  • [2024-04-22 23:28:27]
  • 提交

answer

#include <bits/stdc++.h>
using namespace std;
template <int MOD> struct modint{
	static const int mod = MOD;
	int v;
	operator int() const {return v;} 
	modint():v(0) {}
	modint(long long _v):v(int(_v%MOD)) {v += (v<0)*MOD;}
	modint& operator+=(modint o) {if((v += o.v) >= MOD) v-= MOD; return *this;}
	modint& operator-=(modint o) {if((v -= o.v) < 0) v+= MOD; return *this;}
	modint& operator*=(modint o) {v = int((long long)v*o.v%MOD); return *this;}
	modint& operator/=(modint o) {return *this *= inv(o);}
	friend modint pow(modint a, long long p) {
		modint ans = 1; assert(p >= 0);
		for (; p; p /= 2, a *= a) if (p&1) ans *= a;
		return ans; }	friend modint inv(modint a) {assert(a.v != 0); return pow(a, MOD-2);}
	friend modint operator+(modint a, modint b) {return a += b;}
	friend modint operator-(modint a, modint b) {return a -= b;}
	friend modint operator*(modint a, modint b) {return a *= b;}
	friend modint operator/(modint a, modint b) {return a /= b;}
};
using mint = modint<(int)998244353>;
const int mx = 5005;
mint dp[2 * mx + 1][2 * mx + 1][2 * mx + 1];

int main() {
    
    int n, pp;
    cin >> n >> pp;
    mint stay = (mint)pp * inv((mint)100);
    mint move = (mint)(100 - pp) * inv((mint)100);
    
    vector <int> arr(n + 1);
    for(int i = 1; i <= n; i++){
        cin >> arr[i];
    }
    
    mint home_prob = inv(mint(2 * n + 1));

    for(int i = 0; i < 2 * n + 1; i++){
        for(int j = 0; j < 2 * n + 1; j++) {
            for(int k = 0; k < 2 * n + 1; k++) {
                dp[i][j][k] = 0;
            }
        }
        dp[i][n][0] = home_prob;
    }
    // dp[home][idx][time]

    mint answer = home_prob; // if home = n

    dp[n][n][0] = 0;


    for(int time = 1; time <= n; time++) {
        for(int home = 0; home <= 2 * n; home++) {
            for(int idx = 0; idx <= 2 * n; idx++) {
                dp[home][idx][time] += dp[home][idx][time - 1] * stay;
                if(idx > 0) {
                    if(arr[time] == 1) {
                        dp[home][idx][time] += dp[home][idx - 1][time - 1] * move;
                    }
                }
                if(idx < 2 * n) {
                    if(arr[time] == -1) {
                        dp[home][idx][time] += dp[home][idx + 1][time - 1] * move;
                    }
                }
                if(home == idx) {
                    // cout << idx << " " << dp[home][idx][time] << endl;
                    answer += dp[home][idx][time];
                    // if(home == 6 and idx == 6 and time == 2) cout << dp[home][idx][time] << endl;
                    // if(dp[home][idx][time] != 0) cout << dp[home][idx][time] << " home:" << home << " time:" << time << endl;
                    dp[home][idx][time] = 0;
                }
            }
        } 
    }

    cout << ((int)answer) << endl;
    // cout << home_prob << endl;
    // cout << "prob:" << (home_prob * stay * stay * move) << endl;
}

/*

d.in: 340310575

a.in: 702764025

*/



詳細信息

Test #1:

score: 0
Runtime Error

input:

2 28
1 1

output:


result: