結果

問題 No.391 CODING WAR
ユーザー legosukelegosuke
提出日時 2019-09-28 00:13:00
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 23 ms / 2,000 ms
コード長 3,993 bytes
コンパイル時間 3,896 ms
コンパイル使用メモリ 167,660 KB
実行使用メモリ 5,916 KB
最終ジャッジ日時 2023-10-25 08:03:47
合計ジャッジ時間 3,062 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 5 ms
5,916 KB
testcase_01 AC 5 ms
5,916 KB
testcase_02 AC 5 ms
5,916 KB
testcase_03 AC 5 ms
5,916 KB
testcase_04 AC 5 ms
5,916 KB
testcase_05 AC 5 ms
5,916 KB
testcase_06 AC 5 ms
5,916 KB
testcase_07 AC 5 ms
5,916 KB
testcase_08 AC 6 ms
5,916 KB
testcase_09 AC 23 ms
5,916 KB
testcase_10 AC 22 ms
5,916 KB
testcase_11 AC 5 ms
5,916 KB
testcase_12 AC 5 ms
5,916 KB
testcase_13 AC 20 ms
5,916 KB
testcase_14 AC 18 ms
5,916 KB
testcase_15 AC 19 ms
5,916 KB
testcase_16 AC 14 ms
5,916 KB
testcase_17 AC 16 ms
5,916 KB
testcase_18 AC 12 ms
5,916 KB
testcase_19 AC 13 ms
5,916 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#define int long long
#define lint long long
#define pii pair<int,int>
#define FOR(i,a,b) for(int i=(a);i<(b);++i)
#define REP(i,n) FOR(i,0,n)
#define ALL(c) (c).begin(),(c).end()
#define ZERO(a) memset(a,0,sizeof(a))
#define MINUS(a) memset(a,0xff,sizeof(a))
#define MINF(a) memset(a,0x3f,sizeof(a))
#define POW(n) (1LL<<(n))
#define POPCNT(n) (__builtin_popcount(n))
#define IN(i,a,b) (a <= i && i <= b)
using namespace std;
template <typename T> inline bool CHMIN(T& a,T b) { if(a>b) { a=b; return 1; } return 0; }
template <typename T> inline bool CHMAX(T& a,T b) { if(a<b) { a=b; return 1; } return 0; }
template <typename T> inline void SORT(T& a) { sort(ALL(a)); }
template <typename T> inline void REV(T& a) { reverse(ALL(a)); }
template <typename T> inline void UNI(T& a) { sort(ALL(a)); a.erase(unique(ALL(a)),a.end()); }
template <typename T> inline T LB(vector<T>& v, T a) { return *lower_bound(ALL(v),a); }
template <typename T> inline int LBP(vector<T>& v, T a) { return lower_bound(ALL(v),a) - v.begin(); }
template <typename T> inline T UB(vector<T>& v, T a) { return *upper_bound(ALL(v),a); }
template <typename T> inline int UBP(vector<T>& v, T a) { return upper_bound(ALL(v),a) - v.begin(); }
template <typename T1, typename T2> ostream& operator<< (ostream& os, const pair<T1,T2>& p) { os << p.first << " " << p.second; return os; }
template <typename T1, typename T2> istream& operator>> (istream& is, pair<T1,T2>& p) { is >> p.first >> p.second; return is; }
template <typename T> ostream& operator<< (ostream& os, const vector<T>& v) { REP(i,v.size()) { if (i) os << " "; os << v[i]; } return os; }
template <typename T> istream& operator>> (istream& is, vector<T>& v) { for(T& in : v) is >> in; return is; }
template <typename T = int> vector<T> make_v(size_t a) { return vector<T>(a); }
template <typename T, typename... Ts> auto make_v(size_t a, Ts... ts) { return vector<decltype(make_v<T>(ts...))>(a,make_v<T>(ts...)); }
template <typename T, typename V> typename enable_if<is_class<T>::value == 0>::type fill_v(T &t, const V &v) { t = v; }
template <typename T, typename V> typename enable_if<is_class<T>::value != 0>::type fill_v(T &t, const V &v) { for(auto &e : t) fill_v(e,v); }
const lint MOD = 1000000007;
const lint INF = 0x3f3f3f3f3f3f3f3f;
const double EPS = 1e-10;

const int MAX = 100010;
long long inv[MAX];
long long fac[MAX];
long long finv[MAX];

void init() {
    fac[0] = fac[1] = 1;
    finv[0] = finv[1] = 1;
    inv[1] = 1;
    for (long long i = 2; i < MAX; i++) {
        fac[i] = fac[i - 1] * i % MOD;
        inv[i] = MOD - inv[MOD % i] * (MOD / i) % MOD;
        finv[i] = finv[i - 1] * inv[i] % MOD;
    }
}

// 累乗 a^n
long long pow_mod(long long a, long long n) {
    long long res = 1;
    while (n > 0) {
        if (n & 1) res = res * a % MOD;
        a = a * a % MOD;
        n >>= 1;
    }
    return res;
}

// nCk (階乗 + 逆元)
// n < 10^7, k < 10^7, n,k < MOD, MOD: 素数
// 前計算 O(n) 時間
// 本計算 O(1) 時間
long long nCk(int n, int k) {
    if (n < k) return 0;
    if (n < 0 || k < 0) return 0;
    return fac[n] * finv[k] % MOD * finv[n - k] % MOD;
}

// 包除原理 (n: 玉の数, k: 箱の数)
// O(k log n)
long long PIE(int n, int k) {
    if (k > n) return 0;
    if (n < 0 || k < 0) return 0;
    long long res = 0;
    for (int i = 0; i <= k; i++) {
        long long tmp = nCk(k, i) * pow_mod(i, n) % MOD;
        if ((k - i) & 1) res = (res + MOD - tmp) % MOD;
        else res = (res + tmp) % MOD;
    }
    return res;
}

void _main() {
    init();
    int N, M;
    cin >> N >> M;
    cout << PIE(N, M) << endl;
}

signed main(signed argc, char **argv) {
    if (argc > 1) {
        if (strchr(argv[1], 'i'))
            freopen("input.txt", "r", stdin);
        if (strchr(argv[1], 'o'))
            freopen("output.txt", "w", stdout);
    }
    cin.tie(nullptr);
    ios::sync_with_stdio(false);
    cout << fixed << setprecision(10);
    _main();
    return 0;
}
0