#413564 (C++17) No.573 a^2[i] = a[i]

提出ソース

結果

問題	No.573 a^2[i] = a[i]
ユーザー	legosuke
提出日時	2020-01-02 21:22:22
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	15 ms / 2,000 ms
コード長	4,249 bytes
コンパイル時間	2,188 ms
コンパイル使用メモリ	204,380 KB
実行使用メモリ	6,944 KB
最終ジャッジ日時	2024-05-02 06:44:18
合計ジャッジ時間	3,721 ms
ジャッジサーバーID （参考情報）	judge4 / judge2

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	6 ms 6,812 KB
testcase_01	AC	7 ms 6,944 KB
testcase_02	AC	6 ms 6,940 KB
testcase_03	AC	6 ms 6,940 KB
testcase_04	AC	6 ms 6,940 KB
testcase_05	AC	7 ms 6,940 KB
testcase_06	AC	6 ms 6,940 KB
testcase_07	AC	6 ms 6,940 KB
testcase_08	AC	6 ms 6,940 KB
testcase_09	AC	7 ms 6,944 KB
testcase_10	AC	6 ms 6,944 KB
testcase_11	AC	7 ms 6,944 KB
testcase_12	AC	7 ms 6,940 KB
testcase_13	AC	7 ms 6,944 KB
testcase_14	AC	6 ms 6,940 KB
testcase_15	AC	6 ms 6,940 KB
testcase_16	AC	6 ms 6,944 KB
testcase_17	AC	6 ms 6,940 KB
testcase_18	AC	6 ms 6,940 KB
testcase_19	AC	7 ms 6,940 KB
testcase_20	AC	6 ms 6,944 KB
testcase_21	AC	7 ms 6,940 KB
testcase_22	AC	6 ms 6,944 KB
testcase_23	AC	6 ms 6,940 KB
testcase_24	AC	6 ms 6,940 KB
testcase_25	AC	6 ms 6,940 KB
testcase_26	AC	7 ms 6,944 KB
testcase_27	AC	6 ms 6,940 KB
testcase_28	AC	6 ms 6,944 KB
testcase_29	AC	6 ms 6,944 KB
testcase_30	AC	6 ms 6,940 KB
testcase_31	AC	6 ms 6,944 KB
testcase_32	AC	6 ms 6,940 KB
testcase_33	AC	6 ms 6,944 KB
testcase_34	AC	6 ms 6,940 KB
testcase_35	AC	7 ms 6,940 KB
testcase_36	AC	6 ms 6,944 KB
testcase_37	AC	6 ms 6,940 KB
testcase_38	AC	6 ms 6,944 KB
testcase_39	AC	6 ms 6,944 KB
testcase_40	AC	7 ms 6,940 KB
testcase_41	AC	7 ms 6,940 KB
testcase_42	AC	6 ms 6,940 KB
testcase_43	AC	6 ms 6,940 KB
testcase_44	AC	7 ms 6,940 KB
testcase_45	AC	7 ms 6,940 KB
testcase_46	AC	15 ms 6,940 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#include <bits/stdc++.h>
#define endl '\n'
#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 FORR(i, a, b) for(int i = (a); i > (b); --i)
#define REP(i, n) FOR(i, 0, n)
#define SZ(v) ((int)v.size())
#define ZERO(a) memset(a, 0, sizeof(a))
#define MINUS(a) memset(a, 0xff, sizeof(a))
#define FILLINF(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(begin(a), end(a)); }
template <typename T> inline void REV(T& a) { reverse(begin(a), end(a)); }
template <typename T> inline void UNI(T& a) { sort(a.begin(), a.end()); a.erase(unique(a.begin(), a.end()), a.end()); }
template <typename S, typename T> inline T LB(S& v, T a) { return *lower_bound(begin(v), end(v), a); }
template <typename S, typename T> inline int LBP(S& v, T a) { return lower_bound(begin(v), end(v), a) - begin(v); }
template <typename S, typename T> inline T UB(S& v, T a) { return *upper_bound(begin(v), end(v), a); }
template <typename S, typename T> inline int UBP(S& v, T a) { return upper_bound(begin(v), end(v), a) - begin(v); }
template <typename S, typename T> ostream& operator << (ostream& os, const pair<S, T>& p) { os << p.first << " " << p.second; return os; }
template <typename S, typename T> istream& operator >> (istream& is, pair<S, T>& p) { is >> p.first >> p.second; return is; }
template <typename T> ostream& operator << (ostream& os, const vector<T>& v) { for (int i = 0; i < v.size(); ++i) { 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 double PI = 3.14159265358979323846;
const double EPS = 1e-10;
const lint INF = 0x3f3f3f3f3f3f3f3f;
const lint MOD = 1000000007;

template <int MOD> class Combination {
    vector<long long> inv, fac, finv;

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

public:
    Combination(int n) : inv(n + 1), fac(n + 1), finv(n + 1) { init(n); }

    long long P(int n, int k) {
        if (n < k || n < 0 || k < 0) return 0;
        return fac[n] * finv[n - k] % MOD;
    }

    long long C(int n, int k) {
        if (n < k || n < 0 || k < 0) return 0;
        return fac[n] * finv[k] % MOD * finv[n - k] % MOD;
    }

    long long H(int n, int k) {
        if (n < 0 || k < 0) return 0;
        return C(n + k - 1, k);
    }
};

// O(log n) time
long long mod_pow(long long a, long long n, long long m) {
    a %= m;
    if (a < 0) a += m;
    long long res = 1;
    while (n > 0) {
        if (n & 1) res = res * a % m;
        a = a * a % m;
        n >>= 1;
    }
    return res;
}

int N;

void _main() {
    cin >> N;
    Combination<MOD> C(100000);
    int ans = 0;
    FOR (k, 1, N + 1) {
        (ans += C.C(N, k) * mod_pow(k, N - k, MOD) % MOD) %= MOD;
    }
    cout << ans << 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);
    }
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cout << fixed << setprecision(10);
    _main();
    return 0;
}

yukicoder

結果

テストケース

ソースコード