結果

問題 No.1102 Remnants
ユーザー Enjapma_kyoproEnjapma_kyopro
提出日時 2020-07-03 22:32:50
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 644 ms / 2,000 ms
コード長 3,739 bytes
コンパイル時間 1,741 ms
コンパイル使用メモリ 172,556 KB
実行使用メモリ 18,220 KB
最終ジャッジ日時 2024-09-17 04:11:58
合計ジャッジ時間 8,659 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <unistd.h>
#include <sys/time.h>
#include <stdlib.h>

using namespace std;

typedef long long ll;
typedef pair<ll, ll> P;
typedef pair<P, ll> T;
typedef pair<long double , ll> Ps;
typedef pair<ll, bool> Pb;

const ll INF = 1e18;
const ll fact_table = 101000000;
long double Pi = 3.1415926535897932384626;

priority_queue <ll> pql;
priority_queue <P> pqp;
priority_queue <P> bag;
//big priority queue
priority_queue <ll, vector<ll>, greater<ll> > pqls;
priority_queue <P, vector<P>, greater<P> > pqps;
//small priority queue
//top pop

ll dx[8] = {1, 0, -1, 0, 1, 1, -1, -1};
ll dy[8] = {0, 1, 0, -1, 1, -1, -1, 1};
char dir[] = "DRUL";
//↓,→,↑,←

/*
#define endl "\n"
#ifdef ENJAPMA
#undef endl
#endif
*/
#define p(x) cout<<x<<endl;
#define el cout<<endl;
#define pe(x) cout<<x<<" ";
#define ps(x) cout<<fixed<<setprecision(25)<<x<<endl;
#define pu(x) cout<<(x);
#define pb push_back
#define lb lower_bound
#define ub upper_bound
#define CLEAR(a) a = decltype(a)();
#define pc(x) cout << x << ",";
#define rep(i, n) for (ll i = 0; i < (n); i ++)

// const ll mod = 998244353ll;
const ll mod = 1000000007ll;

ll mypow(ll a, ll b, ll mod) {ll x = 1; while (b) {while (!(b & 1)) {(a *= a) %= mod; b >>= 1;}(x *= a) %= mod; b--;} return x;}
void YES(bool cond) { if (cond) { p("YES");} else { p("NO");} return;}
void Yes(bool cond) { if (cond) { p("Yes");} else { p("No");} return;}
void line() { p("--------------------"); return;}


ll fact[fact_table + 5], rfact[fact_table + 5];

/*
void c3_init() {
    fact[0] = rfact[0] = 1;
    for (ll i = 1; i <= fact_table; i++) {
        fact[i] = (fact[i - 1] * i) % mod;
    }
    rfact[fact_table] = mypow(fact[fact_table], mod - 2, mod);
    for (ll i = fact_table; i >= 1; i--) {
        rfact[i - 1] = rfact[i] * i;
        rfact[i - 1] %= mod;
    }
    return;
}
ll c3(ll n, ll r) {
    return (((fact[n] * rfact[r]) % mod ) * rfact[n - r]) % mod;
}
*/

bool multicase = false;

ll n, k, num, a, b, c, d, e, h, q, ans, sum;

typedef vector<ll> vec;
typedef vector<vector<ll>> mat;

struct Timer {
    int64_t start;
    const int64_t CYCLES_PER_SEC = 2800000000;
    Timer() { reset(); }
    void reset() { start = getCycle(); }
    inline double get() { return (double)(getCycle() - start) / CYCLES_PER_SEC; }
    inline int64_t getCycle() {
        uint32_t low, high;
        __asm__ volatile ("rdtsc" : "=a" (low), "=d" (high));
        return ((int64_t)low) | ((int64_t)high << 32);
    }
};

Timer timer;

ll x[500005], y[500005], z[500005];
ll fact2[500005];
ll krfact, kfact;

ll calc(ll m){
    ll res = fact2[m-1];
    ll bunbo = kfact * fact[m - 1];
    bunbo %= mod;
    res *= mypow(bunbo, mod - 2, mod);
    res %= mod;
    return res;
}

void solve() {
    cin >> n >> k;
    for(ll i=0;i<n;i++){
        cin >> x[i];
    }
    fact[0] = 1, rfact[0] = 1;
    for(ll i=1;i<=n+10;i++){
        fact[i] = (fact[i-1] * i) % mod;
    }
    rfact[n+10] = mypow(fact[n+10], mod - 2, mod);
    kfact = 1ll;
    for(ll i=1;i<=k;i++){
        kfact *= i;
        kfact %= mod;
    }
    krfact = mypow(kfact, mod - 2, mod);
    fact2[0] = kfact;
    for(ll i=1;i<=n+10;i++){
        fact2[i] = (fact2[i-1] * (i+k)) % mod;
    }


    ll ans = 0;
    for(ll i=0;i<n;i++){
        ll res = calc(i + 1) * calc(n - i);
        res %= mod;
        ans += res * x[i];
        ans %= mod;
    }
    p(ans);
    return;
}

int main() {
    // init();
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    ll q, testcase = 1;
    if (multicase) {
        cin >> q;
    } else {
        q = 1;
    }
    while (q--) {
        // pu("Case ");pu("#");pu(testcase);pu(": ");
        solve();
        testcase++;
    }
    // solve();
    return 0;
}
0