結果

問題 No.1348 Split Tile
ユーザー kk
提出日時 2021-04-18 18:35:33
言語 C++17(1z)
(gcc 10.2.0 + boost 1.73.0)
結果
AC  
実行時間 57 ms / 2,000 ms
コード長 2,737 Byte
コンパイル時間 2,070 ms
使用メモリ 18,568 KB
最終ジャッジ日時 2021-04-18 18:35:38
合計ジャッジ時間 4,864 ms
ジャッジサーバーID
(参考情報)
judge9 / judge8
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
使用メモリ
testcase_00 AC 1 ms
3,292 KB
testcase_01 AC 2 ms
3,444 KB
testcase_02 AC 4 ms
3,812 KB
testcase_03 AC 49 ms
16,908 KB
testcase_04 AC 36 ms
13,084 KB
testcase_05 AC 12 ms
5,808 KB
testcase_06 AC 6 ms
4,212 KB
testcase_07 AC 54 ms
18,504 KB
testcase_08 AC 38 ms
13,724 KB
testcase_09 AC 35 ms
12,516 KB
testcase_10 AC 49 ms
16,976 KB
testcase_11 AC 7 ms
4,640 KB
testcase_12 AC 29 ms
10,632 KB
testcase_13 AC 40 ms
14,324 KB
testcase_14 AC 30 ms
10,936 KB
testcase_15 AC 42 ms
15,152 KB
testcase_16 AC 27 ms
10,280 KB
testcase_17 AC 16 ms
7,020 KB
testcase_18 AC 36 ms
13,616 KB
testcase_19 AC 26 ms
9,512 KB
testcase_20 AC 27 ms
10,052 KB
testcase_21 AC 23 ms
8,740 KB
testcase_22 AC 54 ms
17,036 KB
testcase_23 AC 57 ms
17,612 KB
testcase_24 AC 49 ms
17,168 KB
testcase_25 AC 53 ms
18,280 KB
testcase_26 AC 50 ms
17,572 KB
testcase_27 AC 55 ms
18,568 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>

using namespace std;

template<typename T,T MOD>
struct mod_int{
  T v;
  mod_int(): v(0) {}
  mod_int(T t) { v=t%MOD; if(v<0) v+=MOD; }
  
  mod_int inv() const {
    T a=v, b=MOD, u=1, v=0;
		while (b>0) { T t=a/b; a-=t*b; swap(a,b); u-=t*v; swap(u,v); }
		return u;
  }
  
  mod_int pow(long long p) const {
    mod_int ret = 1, base(v);
    while (p > 0) {
      if (p & 1)
        ret *= base;
      base *= base;
      p >>= 1;
    }
    return ret;
  }
  
  mod_int& operator+=(const mod_int &a) { v+=a.v; if(v>=MOD)v-=MOD; return *this; }
  mod_int& operator-=(const mod_int &a) { v+=MOD-a.v; if(v>=MOD)v-=MOD; return *this; }
  mod_int& operator*=(const mod_int &a) { v=v*a.v%MOD; return *this; }
  mod_int& operator/=(const mod_int &a) { return (*this)*=a.inv(); }
  
  mod_int operator+(const mod_int &a) const { return mod_int(v)+=a; }
  mod_int operator-(const mod_int &a) const { return mod_int(v)-=a; }
  mod_int operator*(const mod_int &a) const { return mod_int(v)*=a; }
  mod_int operator/(const mod_int &a) const { return mod_int(v)/=a; }
  
  friend mod_int operator+(T x, const mod_int& m) { return mod_int(x)+m; }
  friend mod_int operator-(T x, const mod_int& m) { return mod_int(x)-m; }
  friend mod_int operator*(T x, const mod_int& m) { return mod_int(x)*m; }
  friend mod_int operator/(T x, const mod_int& m) { return mod_int(x)/m; }
  
  mod_int operator-() const{ return v?mod_int(MOD-v):mod_int(v); }
  
  bool operator==(const mod_int &a) const { return v==a.v; }
  bool operator!=(const mod_int &a) const { return v!=a.v; }
  bool operator<(const mod_int &a) const { return v <a.v; }
  
	friend istream& operator>>(istream& is, mod_int& m){ T t; is>>t; m=mod_int(t); return is; }
	friend ostream& operator<<(ostream& os,const mod_int& m){ return os<<m.v; }
};

template<typename T>
class chooser {
  vector<T> f, g;
public:
  chooser(int n) {
    f.resize(n + 1);
    g.resize(n + 1);
    f[0] = 1;
    for (int i = 1; i <= n; i++) f[i] = i * f[i-1];
    g[n] = 1 / f[n];
    for (int i = n - 1; i >= 0; i--) g[i] = g[i+1] * (i+1);
  }
  
  T perm(int n, int k) {
    return k > n ? 0 : f[n] * g[n-k];
  }
  
  T comb(int n, int k) {
    return k > n ? 0 : f[n] * g[k] * g[n-k];
  }
  
  T homo(int n, int k) {
    return comb(n+k-1, k);
  }
};

using mint = mod_int<long long, 998244353>;

int main() {
  ios_base::sync_with_stdio(0);
  cin.tie(0);

  int n;
  cin >> n;

  mint ret;
  chooser<mint> choo(n);
  
  // 先頭
  for (int i = 1; i < n; i++)
    ret += choo.perm(n-1, i) * choo.perm(n-i, n-i);
  
  // 先頭以外
  for (int i = 1; i < n; i++)
    ret += choo.comb(i, 1) * choo.perm(n-2, i-1) * choo.perm(n-i, n-i) * (n-1);

  cout << ret << endl;
  
  return 0;
}
0