結果
| 問題 | No.3753 Certainly a Cretan |
| コンテスト | |
| ユーザー |
 edon8618
|
| 提出日時 | 2026-10-03 01:55:58 |
| 言語 | C++23 (gcc 15.3.0 + boost 1.92.0 + ACL) |
| 結果 |
AC
不安定
|
| 実行時間 | 618 ms / 2,500 ms |
| + 969µs | |
| コード長 | 8,784 bytes |
| 記録 | |
| コンパイル時間 | 5,162 ms |
| コンパイル使用メモリ | 382,076 KB |
| 実行使用メモリ | 116,836 KB |
| 最終ジャッジ日時 | 2026-10-03 01:56:42 |
| 合計ジャッジ時間 | 15,415 ms |
|
ジャッジサーバーID （参考情報） |
judge1_0 / judge4_0 |
（要ログイン）
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 46 |
ソースコード
// #pragma GCC optimize("O3,unroll-loops")
// #pragma GCC target("avx2")
#include <bits/stdc++.h>
using namespace std;
#include <atcoder/all>
using namespace atcoder;
// #include <boost/multiprecision/cpp_int.hpp>
// using namespace boost::multiprecision;
#define ll long long
#define ld long double
#define rep(i, n) for (ll i = 0; i < (ll)(n); ++i)
#define vi vector<int>
#define vl vector<ll>
#define vd vector<double>
#define vb vector<bool>
#define vs vector<string>
#define vc vector<char>
#define ull unsigned long long
#define all(a) (a).begin(), (a).end()
#define rall(a) (a).rbegin(), (a).rend()
template<class T, class U>
inline bool chmax(T &a, const U &b) {
if (a < b) {
a = b;
return true;
}
return false;
}
template<class T, class U>
inline bool chmin(T &a, const U &b) {
if (a > b) {
a = b;
return true;
}
return false;
}
// #define ll int
// #define ll int128_t
// #define ll int256_t
// #define ll cpp_int
constexpr ll inf = (1ll << 60);
// constexpr ll inf = (1 << 30);
// const double PI=3.1415926535897932384626433832795028841971;
uint32_t xor_x = 123456789, xor_y = 362436069, xor_z = 521288629, xor_w = 88675123;
inline uint32_t xor_next() {
uint32_t t = xor_x ^ (xor_x << 11);
xor_x = xor_y; xor_y = xor_z; xor_z = xor_w;
return xor_w = (xor_w ^ (xor_w >> 19)) ^ (t ^ (t >> 8));
}
inline int rnd(int max_val) { return xor_next() % max_val; }
struct Timer {
std::chrono::steady_clock::time_point start_time;
Timer() {
reset();
}
// 測定の起点リセット用
void reset() {
start_time = std::chrono::steady_clock::now();
}
// スタートからの経過時間をミリ秒(msec)で返す
long long get_ms() const {
auto now = std::chrono::steady_clock::now();
return std::chrono::duration_cast<std::chrono::milliseconds>(now - start_time).count();
}
};
// ll rui(ll a,ll b){
// if(b==0)return 1;
// if(b%2==1) return a*rui(a*a,b/2);
// return rui(a*a,b/2);
// }
// vl fact;
// ll kai(ll n){
// fact.resize(n,1);
// rep(i,n-1)fact[i+1]=fact[i]*(i+1);
// }
// using mint = ld;
using mint = modint998244353;//static_modint<998244353>
// using mint = modint1000000007;//static_modint<1000000007>
// using mint = static_modint<922267487>; // 多分落とされにくい NOT ntt-friendly
// using mint = static_modint<469762049>; // ntt-friendly
// using mint = static_modint<167772161>; // ntt-friendly
// using mint = modint;//mint::set_mod(mod);
// ll const mod=1000000007ll;
// ll const mod=998244353ll;
// ll modrui(ll a,ll b,ll mod){
// a%=mod;
// if(b==0)return 1;
// if(b%2==1) return a*modrui(a*a%mod,b/2,mod)%mod;
// return modrui(a*a%mod,b/2,mod)%mod;
// }
// void incr(vl &v,ll n){// n進法
// ll k=v.size();
// v[k-1]++;
// ll now=k-1;
// while (v[now]>=n)
// {
// v[now]=0;
// if(now==0)break;
// v[now-1]++;
// now--;
// }
// return;
// }
vector<mint> fact,invf;
void init_modfact(ll sz){
fact.resize(sz);
invf.resize(sz);
fact[0]=1;
rep(i,sz-1){
fact[i+1]=fact[i]*(i+1);
}
invf[sz-1]=1/fact[sz-1];
for(ll i=sz-2; i>=0; i--){
invf[i]=invf[i+1]*(i+1);
}
}
mint choose(ll n,ll r){
if(n<r || r<0 || n<0)return 0;
return fact[n]*invf[r]*invf[n-r];
}
// vector<mint> modpow,invpow;
// void init_modpow(ll x,ll sz){
// mint inv=1/mint(x);
// modpow.assign(sz,1);
// invpow.assign(sz,1);
// rep(i,sz-1){
// modpow[i+1]=modpow[i]*x;
// invpow[i+1]=invpow[i]*inv;
// }
// }
// long long phi(long long n) {// O(sqrt(n))
// long long res = n;
// for (long long i = 2; i * i <= n; i++) {
// if (n % i == 0) {
// res -= res / i;
// while (n % i == 0) n /= i;
// }
// }
// if (n > 1) res -= res / n;
// return res;
// }
mint ctl(ll n){
ll k=n/2;
return choose(n,k)-choose(n,k-1);
}
struct S{
bool lY,rY,suc;
mint cnt;
ll lcnt,rcnt;
bool e;
};
S op(S l,S r){
if(l.e)return r;
if(r.e)return l;
S ret;
ret.e=0;
ret.lY=l.lY;
ret.rY=r.rY;
ret.cnt=l.cnt*r.cnt;
ret.suc=l.suc && r.suc && (l.rY==r.lY);
if(ret.suc){
ret.lcnt=ret.rcnt=l.lcnt+r.rcnt;
return ret;
}
if(l.suc && r.suc){
ret.lcnt=l.lcnt;
ret.rcnt=r.rcnt;
}
else if(l.suc && !r.suc){
if(l.rY==r.lY){
ret.lcnt=l.lcnt+r.lcnt;
ret.rcnt=r.rcnt;
}
else{
ret.lcnt=l.lcnt;
ret.rcnt=r.rcnt;
ret.cnt*=ctl(r.lcnt);
}
}
else if(!l.suc && r.suc){
if(l.rY==r.lY){
ret.lcnt=l.lcnt;
ret.rcnt=l.rcnt+r.rcnt;
}
else{
ret.lcnt=l.lcnt;
ret.rcnt=r.rcnt;
ret.cnt*=ctl(l.rcnt);
}
}
else{
if(l.rY==r.lY){
ret.lcnt=l.lcnt;
ret.rcnt=r.rcnt;
ret.cnt*=ctl(l.rcnt+r.lcnt);
}
else{
ret.lcnt=l.lcnt;
ret.rcnt=r.rcnt;
ret.cnt*=ctl(l.rcnt);
ret.cnt*=ctl(r.lcnt);
}
}
return ret;
}
S e(){
S s;
s.e=1;
return s;
}
void solve(){
ll n,q;
string s;
cin >> n >> q >> s;
init_modfact(n*2+9);
// rep(i,n+1)cout << ctl(i).val() << " \n"[i==n];
// return;
segtree<S,op,e> seg(n);
rep(i,n){
if(s[i]=='Y'){
seg.set(i,S{1,1,1,1,1,1,0});
}
else{
seg.set(i,S{0,0,1,1,1,1,0});
}
}
ll cnt=0;
rep(i,s.size()/2){
if(s[i*2]!=s[i*2+1])cnt++;
}
while(q--){
ll cmd;
cin >> cmd;
if(cmd==1){
ll i;
cin >> i;
i--;
if(s[i]=='Y')s[i]='N';
else s[i]='Y';
if(s[i]=='Y'){
seg.set(i,S{1,1,1,1,1,1,0});
}
else{
seg.set(i,S{0,0,1,1,1,1,0});
}
ll a=i/2;
if(a*2+1==n)continue;
if(s[a*2]==s[a*2+1])cnt--;
else cnt++;
}
else{
ll k;
cin >> k;
// vector<mint> dp(1,1);
// // cout << 1 << endl;
// rep(k,n){
// vector<mint> ep(dp.size()+1,0);
// rep(t,dp.size()){
// if(s[k]=='Y'){
// if(2*t>=k)ep[t+1]+=dp[t];
// if(2*t>=k+1)ep[t]+=dp[t];
// }
// else{
// if(2*t<k)ep[t+1]+=dp[t];
// if(2*t<k+1)ep[t]+=dp[t];
// }
// }
// swap(ep,dp);
// // rep(i,dp.size())cout << dp[i].val() <<" \n"[i==dp.size()-1];
// }
// cout << dp[n-k].val() << endl;
// continue;
// Y only と N onlyをとけばよい
if(cnt!=0){
cout << 0 << endl;
continue;
}
S s=seg.all_prod();
if(s.suc){
if(s.rY){
if(k>n/2)cout << 0 << endl;
else{
cout << (choose(n,k)-choose(n,k-1)).val() << endl;
}
}
else{
ll t=n-k;
if(t>n/2)cout << 0 << endl;
else{
cout << (choose(n,t)-choose(n,t-1)).val() << endl;
}
}
}
else{
s.cnt*=ctl(s.lcnt);
ll m=s.rcnt;
ll u=k-(n-m)/2;
if(s.rY){
if(u>m/2)cout << 0 << endl;
else{
cout << (s.cnt*(choose(m,u)-choose(m,u-1))).val() << endl;
}
}
else{
ll t=m-u;
if(t>m/2)cout << 0 << endl;
else{
cout << (s.cnt*(choose(m,t)-choose(m,t-1))).val() << endl;
}
}
}
}
}
}
int main(){
ios::sync_with_stdio(false);
std::cin.tie(nullptr);
// ll mx=1234;
// vc fl(mx+1,0);
// for(ll d=2;d<=mx;d++){
// if(fl[d])continue;
// ll x=d;
// ps.push_back(x);
// while(x<=mx){
// fl[x]=1;
// x+=d;
// }
// }
ll t = 1;
// cin >> t;
while (t--){
solve();
cout << flush;
}
}