結果
| 問題 | 
                            No.3178 free sort
                             | 
                    
| コンテスト | |
| ユーザー | 
                             iomir
                         | 
                    
| 提出日時 | 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();
}
    
            
            
            
        
            
iomir