#1027062 (C++23) No.2964 Obstruction Bingo

提出ソース

結果

問題	No.2964 Obstruction Bingo
ユーザー	hato336
提出日時	2024-11-16 16:18:36
言語	C++23 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	1,731 ms / 2,468 ms
コード長	7,446 bytes
コンパイル時間	6,447 ms
コンパイル使用メモリ	327,120 KB
実行使用メモリ	6,820 KB
最終ジャッジ日時	2024-11-16 16:20:26
合計ジャッジ時間	57,485 ms
ジャッジサーバーID （参考情報）	judge1 / judge3

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 6,816 KB
testcase_01	AC	2 ms 6,816 KB
testcase_02	AC	10 ms 6,820 KB
testcase_03	AC	2 ms 6,820 KB
testcase_04	AC	2 ms 6,816 KB
testcase_05	AC	396 ms 6,820 KB
testcase_06	AC	777 ms 6,816 KB
testcase_07	AC	79 ms 6,816 KB
testcase_08	AC	783 ms 6,816 KB
testcase_09	AC	452 ms 6,816 KB
testcase_10	AC	1,061 ms 6,820 KB
testcase_11	AC	183 ms 6,816 KB
testcase_12	AC	1,189 ms 6,816 KB
testcase_13	AC	547 ms 6,816 KB
testcase_14	AC	937 ms 6,820 KB
testcase_15	AC	144 ms 6,816 KB
testcase_16	AC	219 ms 6,816 KB
testcase_17	AC	137 ms 6,820 KB
testcase_18	AC	369 ms 6,816 KB
testcase_19	AC	176 ms 6,820 KB
testcase_20	AC	844 ms 6,820 KB
testcase_21	AC	858 ms 6,816 KB
testcase_22	AC	87 ms 6,820 KB
testcase_23	AC	32 ms 6,820 KB
testcase_24	AC	20 ms 6,816 KB
testcase_25	AC	1,678 ms 6,820 KB
testcase_26	AC	1,656 ms 6,820 KB
testcase_27	AC	1,706 ms 6,816 KB
testcase_28	AC	1,651 ms 6,816 KB
testcase_29	AC	1,672 ms 6,816 KB
testcase_30	AC	1,694 ms 6,820 KB
testcase_31	AC	1,669 ms 6,816 KB
testcase_32	AC	1,668 ms 6,816 KB
testcase_33	AC	1,685 ms 6,816 KB
testcase_34	AC	1,653 ms 6,820 KB
testcase_35	AC	1,666 ms 6,816 KB
testcase_36	AC	1,642 ms 6,816 KB
testcase_37	AC	1,693 ms 6,816 KB
testcase_38	AC	1,731 ms 6,820 KB
testcase_39	AC	1,660 ms 6,820 KB
testcase_40	AC	1,673 ms 6,820 KB
testcase_41	AC	1,648 ms 6,816 KB
testcase_42	AC	1,668 ms 6,816 KB
testcase_43	AC	1,683 ms 6,816 KB
testcase_44	AC	1,686 ms 6,816 KB
testcase_45	AC	747 ms 6,816 KB
testcase_46	AC	647 ms 6,816 KB
testcase_47	AC	649 ms 6,816 KB
testcase_48	AC	961 ms 6,820 KB
testcase_49	AC	934 ms 6,820 KB
testcase_50	AC	1,193 ms 6,816 KB
testcase_51	AC	1,629 ms 6,820 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
#define all(...) std::begin(__VA_ARGS__), std::end(__VA_ARGS__)
#define rall(...) std::rbegin(__VA_ARGS__), std::rend(__VA_ARGS__)
#define OVERLOAD_REP(_1, _2, _3, _4, name, ...) name
#define REP1(n) for(ll i=0;i<n;i++)
#define REP2(i, n) for (ll i=0;i<n;i++)
#define REP3(i, a, n) for (ll i=a;i<n;i++)
#define REP4(i, a, b, n) for(ll i=a;i<n;i+=b)
#define rep(...) OVERLOAD_REP(__VA_ARGS__, REP4, REP3, REP2, REP1)(__VA_ARGS__)
#define OVERLOAD_RREP(_1, _2, _3, _4, name, ...) name
#define RREP1(n) for(ll i=n-1;i>=0;i--)
#define RREP2(i, n) for(ll i=n-1;i>=0;i--)
#define RREP3(i, a, n) for(ll i=n-1;i>=a;i--)
#define RREP4(i, a, b, n) for(ll i=n-1;i>=a;i-=b)
#define rrep(...) OVERLOAD_RREP(__VA_ARGS__, RREP4, RREP3, RREP2, RREP1)(__VA_ARGS__)
#define foa(a,v)  (auto& a : (v))
#define uniq(a) sort(all(a));a.erase(unique(all(a)),end(a))
#define len(n) (long long)(n).size()
#define pb push_back
using ll = long long;
using ld = long double;
using ull = unsigned long long;
using vi = vector<int>;
using vvi = vector<vi>;
using vvvi = vector<vvi>;
using vll = vector<ll>;
using vvll = vector<vll>;
using vvvll = vector<vvll>;
using vs = vector<string>;
using vvs = vector<vs>;
using vvvs = vector<vvs>;
using vld = vector<ld>;
using vvld = vector<vld>;
using vvvld = vector<vvld>;
using vc = vector<char>;
using vvc = vector<vc>;
using vvvc = vector<vvc>;
using pll = pair<ll,ll>;
using vpll = vector<pll>;
template<class... T>
constexpr auto min(T... a){
    return min(initializer_list<common_type_t<T...>>{a...});
}
template<class... T>
constexpr auto max(T... a){
    return max(initializer_list<common_type_t<T...>>{a...});
}
template<class... T>
void input(T&... a){
    (cin >> ... >> a);
}
ll POW(ll a,ll b){
    ll ans = 1;
    while (b){
        if (b & 1){
            ans *= a;
        }
        a *= a;
        b /= 2;
    }
    return ans;
}
ll MODPOW(ll a,ll b,ll c){
    ll ans = 1;
    while (b){
        if (b & 1){
            ans *= a;
            ans %= c;
        }
        a *= a;
        a %= c;
        b /= 2;
    }
    return ans;
}
#define OVERLOAD_POW(_1, _2, _3, name, ...) name
#define pow(...) OVERLOAD_POW(__VA_ARGS__, MODPOW, POW)(__VA_ARGS__)
#define INT(...) int __VA_ARGS__; input(__VA_ARGS__)
#define LL(...) ll __VA_ARGS__; input(__VA_ARGS__)
#define ULL(...) ull __VA_ARGS__; input(__VA_ARGS__)
#define LD(...) ld __VA_ARGS__; input(__VA_ARGS__)
#define STR(...) string __VA_ARGS__; input(__VA_ARGS__)
#define CHA(...) char __VA_ARGS__; input(__VA_ARGS__)
#define VLL(name,length) vll name(length);rep(i,length){cin >> name[i];}
#define VVLL(name,h,w) vvll name(h,vll(w));rep(i,h)rep(j,w){cin >> name[i][j];}
#define VVVLL(name,a,b,c) vvvll name(a,vvll(b,vll(c)));rep(i,a)rep(j,b)rep(k,c){cin >> name[i][j][k];}
#define VI(name,length) vi name(length);rep(i,length){cin >> name[i];}
#define VVI(name,h,w) vvi name(h,vi(w));rep(i,h)rep(j,w){cin >> name[i][j];}
#define VVVI(name,a,b,c) vvvi name(a,vvll(b,vi(c)));rep(i,a)rep(j,b)rep(k,c){cin >> name[i][j][k];}
#define VLD(name,length) vld name(length);rep(i,length){cin >> name[i];}
#define VVLD(name,h,w) vvld name(h,vld(w));rep(i,h)rep(j,w){cin >> name[i][j];}
#define VVVLD(name,a,b,c) vvvld name(a,vvld(b,vld(c)));rep(i,a)rep(j,b)rep(k,c){cin >> name[i][j][k];}
#define VC(name,length) vc name(length);rep(i,length){cin >> name[i];}
#define VVC(name,h,w) vvc name(h,vc(w));rep(i,h)rep(j,w){cin >> name[i][j];}
#define VVVC(name,a,b,c) vvvc name(a,vvc(b,vc(c)));rep(i,a)rep(j,b)rep(k,c){cin >> name[i][j][k];}
#define VS(name,length) vs name(length);rep(i,length){cin >> name[i];}
#define VVS(name,h,w) vvs name(h,vs(w));rep(i,h)rep(j,w){cin >> name[i][j];}
#define VVVS(name,a,b,c) vvvs name(a,vvs(b,vs(c)));rep(i,a)rep(j,b)rep(k,c){cin >> name[i][j][k];}
#define PLL(name) pll name;cin>>name.first>>name.second;
#define VPLL(name,length) vpll name(length);rep(i,length){cin>>name[i].first>>name[i].second;}

void print(){cout << "\n";}
template<class T, class... Ts>
void print(const T& a, const Ts&... b){cout << a;(cout << ... << (cout << ' ', b));cout << '\n';}
void print(vll x){rep(i,len(x)){cout << x[i];if(i!=len(x)-1){cout << " ";}else{cout << '\n';}}}
void print(vvll x){rep(i,len(x))rep(j,len(x[i])){cout << x[i][j];if(j!=len(x[i])-1){cout << " ";}else{cout << '\n';}}}
void print(vi x){rep(i,len(x)){cout << x[i];if(i!=len(x)-1){cout << " ";}else{cout << '\n';}}}
void print(vvi x){rep(i,len(x))rep(j,len(x[i])){cout << x[i][j];if(j!=len(x[i])-1){cout << " ";}else{cout << '\n';}}}
void print(vvvi x){rep(i,len(x))rep(j,len(x[i]))rep(k,len(x[i][j])){cout << x[i][j][k];if(k!=len(x[i][j])-1){cout << " ";}else if(j!=len(x[i])-1){cout << " | ";}else{cout << '\n';}}}
void print(vld x){rep(i,len(x)){cout << x[i];if(i!=len(x)-1){cout << " ";}else{cout << '\n';}}}
void print(vvld x){rep(i,len(x))rep(j,len(x[i])){cout << x[i][j];if(j!=len(x[i])-1){cout << " ";}else{cout << '\n';}}}
void print(vvvld x){rep(i,len(x))rep(j,len(x[i]))rep(k,len(x[i][j])){cout << x[i][j][k];if(k!=len(x[i][j])-1){cout << " ";}else if(j!=len(x[i])-1){cout << " | ";}else{cout << '\n';}}}
void print(vc x){rep(i,len(x)){cout << x[i];if(i!=len(x)-1){cout << " ";}else{cout << '\n';}}}
void print(vvc x){rep(i,len(x))rep(j,len(x[i])){cout << x[i][j];if(j!=len(x[i])-1){cout << " ";}else{cout << '\n';}}}
void print(vvvc x){rep(i,len(x))rep(j,len(x[i]))rep(k,len(x[i][j])){cout << x[i][j][k];if(k!=len(x[i][j])-1){cout << " ";}else if(j!=len(x[i])-1){cout << " | ";}else{cout << '\n';}}}
void print(vs x){rep(i,len(x)){cout << x[i];if(i!=len(x)-1){cout << " ";}else{cout << '\n';}}}
void print(vvs x){rep(i,len(x))rep(j,len(x[i])){cout << x[i][j];if(j!=len(x[i])-1){cout << " ";}else{cout << '\n';}}}
void print(vvvs x){rep(i,len(x))rep(j,len(x[i]))rep(k,len(x[i][j])){cout << x[i][j][k];if(k!=len(x[i][j])-1){cout << " ";}else if(j!=len(x[i])-1){cout << " | ";}else{cout << '\n';}}}
void print(pll x){cout << x.first << x.second << '\n';}
void print(vpll x){rep(i,len(x)){cout << x[i].first << x[i].second << '\n';}}
using mint = modint998244353;
int main(){
    LL(l,k);
    STR(s);
    STR(t);
    VLL(a,26);
    ll su = 0;
    rep(i,26){
        su += a[i];
    }
    mint inv = mint(su).inv();
    vector<vector<mint>> dp(k+2,vector<mint>(k+2,mint(0)));
    dp[0][0] = 1;
    mint ans1 = 0;
    mint ans2 = 0;
    rep(i,k+1){
        vector<vector<mint>> ndp(k+2,vector<mint>(k+2,mint(0)));
        rep(j,k+1){
            rep(m,k+1){
                if(abs(j-m) == l){
                    if(j > m){
                        ans1 += dp[j][m];
                    }
                    else{
                        ans2 += dp[j][m];
                    }
                    continue;
                }
                if(abs(j-m) > l){
                    continue;
                }
                if(s[j%l] != t[m%l]){
                    ndp[j][m] += dp[j][m] * (su - (a[s[j%l] - 'a']) - (a[t[m%l] - 'a'])) * inv;
                    ndp[j+1][m] += dp[j][m] * (a[s[j%l] - 'a']) * inv;
                    ndp[j][m+1] += dp[j][m] * (a[t[m%l] - 'a']) * inv;
                }
                else{
                    ndp[j][m] += dp[j][m] * (su - (a[s[j%l] - 'a'])) * inv;
                    ndp[j+1][m+1] += dp[j][m] * (a[t[m%l] - 'a']) * inv;

                }
            }
            
        }
        dp = ndp;
    }
    print(ans1.val(),ans2.val());
}

yukicoder

結果

テストケース

ソースコード