結果

問題 No.1138 No Bingo!
ユーザー mfbgjsczmfbgjscz
提出日時 2020-07-19 13:01:08
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 71 ms / 3,000 ms
コード長 1,950 bytes
コンパイル時間 1,977 ms
コンパイル使用メモリ 203,196 KB
実行使用メモリ 6,284 KB
最終ジャッジ日時 2023-09-16 21:13:02
合計ジャッジ時間 4,107 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,376 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 36 ms
4,640 KB
testcase_04 AC 71 ms
6,224 KB
testcase_05 AC 2 ms
4,376 KB
testcase_06 AC 1 ms
4,376 KB
testcase_07 AC 2 ms
4,380 KB
testcase_08 AC 1 ms
4,380 KB
testcase_09 AC 40 ms
4,904 KB
testcase_10 AC 65 ms
5,896 KB
testcase_11 AC 70 ms
6,284 KB
testcase_12 AC 23 ms
4,380 KB
testcase_13 AC 6 ms
4,380 KB
testcase_14 AC 70 ms
5,892 KB
testcase_15 AC 30 ms
4,376 KB
testcase_16 AC 63 ms
6,016 KB
testcase_17 AC 62 ms
5,760 KB
testcase_18 AC 61 ms
5,636 KB
testcase_19 AC 42 ms
4,908 KB
testcase_20 AC 62 ms
6,028 KB
testcase_21 AC 60 ms
5,756 KB
testcase_22 AC 56 ms
5,496 KB
testcase_23 AC 3 ms
4,376 KB
testcase_24 AC 41 ms
4,852 KB
testcase_25 AC 26 ms
4,376 KB
testcase_26 AC 22 ms
4,380 KB
testcase_27 AC 36 ms
4,632 KB
testcase_28 AC 22 ms
4,380 KB
testcase_29 AC 23 ms
4,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
/*#include <boost/multiprecision/cpp_dec_float.hpp>
#include <boost/multiprecision/cpp_int.hpp>
*/#pragma GCC target("avx")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")    
using namespace std;
/*namespace mp = boost::multiprecision;
using Bint = mp::cpp_int;
using Real = mp::number<mp::cpp_dec_float<1024>>;
*/#define lli long long int
#define uli unsigned long long int
#define INF 999999999999
#define rep(i,m,n) for(lli i = m;i < n;i++)
#define rrep(i,m,n) for(lli i=m-1;i>=n;i--)
#define pb(n) push_back(n)
#define UE(N) N.erase(unique(N.begin(),N.end()),N.end());
#define Sort(n) sort(n.begin(), n.end())
#define Rev(n) reverse(n.begin(),n.end())
#define Out(S) cout << S << endl
#define NeOut(S) cout << S
#define HpOut(S) cout << setprecision(50) << S << endl
#define Vec(K,L,N,S) vector<L> K(N,S)
#define DV(K,L,N,M,S) vector<vector<L>> K(N,vector<L>(M,S))
#define TV(K,L,N,M,R,S) vector<vector<vector<L>>> K(N,vector<vector<L>>(M,vector<L>(R,S)))
#define pint pair<lli,lli>
#define paf(L,R) pair<L,R>
#define mod 998244353
#define MAX 10000000
#define ALL(a)  a.begin(),a.end()
#define chmax(a, b) a = (((a)<(b)) ? (b) : (a))
#define chmin(a, b) a = (((a)>(b)) ? (b) : (a))
long long modinv(long long a, long long m) {
    long long b = m, u = 1, v = 0;
    while (b) {
        long long t = a/b;
        a -= t*b; swap(a, b);
        u -= t*v; swap(u, v);
    }
    u %= m;
    if (u < 0) u += m;
    return u;
}
int main(){
  ios::sync_with_stdio(false);
  cin.tie(0);
  long long N,sum=0;
  cin >> N;
  vector<long long>dp1(N+1,1);
  vector<long long>dp2(N+1,0);
  dp2[0]=1;
  rep(i,2,N+1)dp1[i]=dp1[i-1]*i%mod;
  rep(i,2,N+1)sum=(sum+mod+(long long)pow(-1,i%2)*(dp1[N]*modinv(dp1[i],mod)%mod))%mod;
  rep(i,4,N+1){
    if(i&1)dp2[i]=((i-1)*dp2[i-1]%mod+2*(i-1)*dp2[i-2]%mod)%mod;
    else dp2[i]=((i-1)*dp2[i-1]%mod+2*(i-2)*dp2[i-4]%mod)%mod;
  }
  Out((dp1[N]-2*sum+dp2[N]+10LL*mod)%mod);
}
0