結果
問題 | No.1419 Power Moves |
ユーザー |
![]() |
提出日時 | 2022-02-16 10:12:54 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 202 ms / 2,000 ms |
コード長 | 1,993 bytes |
コンパイル時間 | 919 ms |
コンパイル使用メモリ | 87,160 KB |
実行使用メモリ | 7,296 KB |
最終ジャッジ日時 | 2024-06-29 07:12:42 |
合計ジャッジ時間 | 6,387 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 31 |
ソースコード
#include <iostream>#include <cstdio>#include <string>#include <algorithm>#include <utility>#include <cmath>#include <vector>#include <stack>#include <queue>#include <deque>#include <set>#include <map>#include <tuple>#include <numeric>#include <functional>using namespace std;typedef long long ll;typedef vector<ll> vl;typedef vector<vector<ll>> vvl;typedef pair<ll, ll> P;#define rep(i, n) for(ll i = 0; i < n; i++)#define exrep(i, a, b) for(ll i = a; i <= b; i++)#define out(x) cout << x << endl#define exout(x) printf("%.10f\n", x)#define chmax(x, y) x = max(x, y)#define chmin(x, y) x = min(x, y)#define all(a) a.begin(), a.end()#define rall(a) a.rbegin(), a.rend()#define pb push_back#define re0 return 0const ll mod = 1000000007;const ll INF = 1e16;const ll MAX = 100010;ll fact[MAX]; // fact[i] : iの階乗のmodll inv[MAX]; // inv[i] : iの逆数のmodll invfact[MAX]; // invfact[i] : iの階乗の逆数のmodvoid init() {fact[0] = 1;inv[0] = inv[1] = 1;invfact[0] = 1;for(ll i = 1; i < MAX; i++) {fact[i] = i * fact[i-1] % mod;if(i >= 2) {inv[i] = mod - inv[mod % i] * (mod / i) % mod;}invfact[i] = invfact[i-1] * inv[i] % mod;}}// a^n (mod.MOD)を求める。計算量はO(logn)ll modpow(ll a, ll n, ll MOD = mod) {if(n == 0) {return 1;}if(n%2 == 1) {return (a * modpow(a, n-1, MOD)) % MOD;}return (modpow(a, n/2, MOD) * modpow(a, n/2, MOD)) % MOD;}int main() {ll n, k;cin >> n >> k;init();vl b(n);exrep(i, 1, n) {b[(2*i - 1) % n]++;}ll x = modpow(2LL, k-1, n);vl a(n);rep(i, n) {a[i] = b[i] * (modpow(2LL, k-1) + mod - x) % mod * inv[n] % mod;}exrep(i, 1, x) {a[(2*i - 1) % n]++;a[(2*i - 1) % n] %= mod;}rep(i, n) {cout << (a[i] + a[(n - i) % n]) * modpow(inv[2], k) % mod << "\n";}re0;}