結果
問題 |
No.3178 free sort
|
ユーザー |
![]() |
提出日時 | 2025-06-13 21:54:45 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 101 ms / 2,000 ms |
コード長 | 3,860 bytes |
コンパイル時間 | 2,793 ms |
コンパイル使用メモリ | 281,976 KB |
実行使用メモリ | 27,204 KB |
最終ジャッジ日時 | 2025-06-13 21:54:56 |
合計ジャッジ時間 | 8,192 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 5 |
other | AC * 40 |
ソースコード
#include<bits/stdc++.h> //#include <atcoder/all> using namespace std; //using namespace atcoder; #define all(v) v.begin(),v.end() using ll = long long; using ull = unsigned long long; using lll = __int128; using vll=vector<ll>; using vvll = vector<vector<ll>>; using P = pair<ll,ll>; using vp=vector<pair<ll, ll>>; //using mint=modint1000000007; //using mint=modint998244353; const ll INF=1ll<<60; ll mod10=1e9+7; ll mod99=998244353; const double PI = acos(-1); #define rep(i,n) for (ll i=0;i<n;++i) #define per(i,n) for(ll i=n-1;i>=0;--i) #define rep2(i,a,n) for (ll i=a;i<n;++i) #define per2(i,a,n) for (ll i=a;i>=n;--i) template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return true; } return false; } template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return true; } return false; } template <long long MOD> struct modint { long long value; modint(long long x = 0) { if (x >= 0) { value = x % MOD; } else { value = MOD - (-x) % MOD; } } modint operator-() const { return modint(-value); } modint operator+() const { return modint(*this); } modint &operator+=(const modint &other) { value += other.value; if (value >= MOD) { value -= MOD; } return *this; } modint &operator-=(const modint &other) { value += MOD - other.value; if (value >= MOD) { value -= MOD; } return *this; } modint &operator*=(const modint other) { value = value * other.value % MOD; return *this; } modint &operator/=(modint other) { (*this) *= other.inv(); return *this; } modint operator+(const modint &other) const { return modint(*this) += other; } modint operator-(const modint &other) const { return modint(*this) -= other; } modint operator*(const modint &other) const { return modint(*this) *= other; } modint operator/(const modint &other) const { return modint(*this) /= other; } modint pow(long long x) const { modint ret(1), mul(value); while (x > 0) { if (x % 2 == 1) { ret *= mul; } mul *= mul; x /= 2; } return ret; } modint inv() const { return pow(MOD - 2); } bool operator==(const modint &other) const { return value == other.value; } bool operator!=(const modint &other) const { return value != other.value; } friend std::ostream &operator<<(std::ostream &os, const modint &x) { return os << x.value; } friend std::istream &operator>>(std::istream &is, modint &x) { long long v; is >> v; x = modint<MOD>(v); return is; } }; using mint = modint<998244353>; using mod107 = modint<1000000007>; vector<ll> fact,invfact,inv; void init(ll N,ll mod){ fact.resize(N+5); inv.resize(N+5); invfact.resize(N+5); fact[0]=1;fact[1]=1; invfact[0]=1;invfact[1]=1; inv[0]=1;inv[1]=1; for(int i=2;i<N+3;i++){ fact[i]=fact[i-1]*i%mod; inv[i]=mod-inv[mod%i]*(mod/i)%mod; invfact[i] = invfact[i-1]*inv[i]%mod; } } ll nCk(ll n,ll x,ll mod){ return fact[n]*invfact[n-x]%mod*invfact[x]%mod; } bool solve(){ string S;cin>>S; ll N=S.size(); vll cnt(10); rep(i,N) cnt[S[i]-'0']++; mint ans=fact[N];rep(i,10) ans/=fact[cnt[i]]; if(cnt[0]==0){ cout << ans << endl; return 0; } mint ans2=(mint)fact[N-1]/fact[cnt[0]-1];rep2(i,1,10) ans2/=fact[cnt[i]]; cout << ans-ans2 << endl; return 0; } int main(){ init(1000000,mod99); cin.tie(0); ios::sync_with_stdio(false); ll T=1;//cin>>T; rep(i,T) solve(); }