結果
問題 | No.1419 Power Moves |
ユーザー |
|
提出日時 | 2021-12-17 13:31:04 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 207 ms / 2,000 ms |
コード長 | 2,813 bytes |
コンパイル時間 | 4,240 ms |
コンパイル使用メモリ | 239,136 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-09-14 14:26:54 |
合計ジャッジ時間 | 9,885 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 31 |
ソースコード
#include <bits/stdc++.h>#include <atcoder/all>typedef long long int ll;typedef long long int ull;#define MP make_pairusing namespace std;using namespace atcoder;typedef pair<ll, ll> P;// const ll MOD = 998244353;const ll MOD = 1000000007;// using mint = modint998244353;using mint = modint1000000007;const double pi = 3.1415926536;const int MAX = 2000003;long long fac[MAX], finv[MAX], inv[MAX];// テーブルを作る前処理void COMinit() {fac[0] = fac[1] = 1;finv[0] = finv[1] = 1;inv[1] = 1;for (int 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;}}// 二項係数計算long long COM(int n, int k){if (n < k) return 0;if (n < 0 || k < 0) return 0;return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD;}void comp(vector<ll>&a){set<ll>s(a.begin(),a.end());map<ll,ll>d;ll cnt=0;for(auto x:s)d[x]=cnt++;for(auto&x:a)x=d[x];}ll gcd (ll x, ll y) {if (y == 0) return x;else if (y > x) {return gcd (y, x);}else return gcd(x % y, y);}ll lcm (ll x, ll y) {return x / gcd(x, y) * y;}ll my_sqrt(ll x) {ll m = 0;ll M = 1000000001;while (M - m > 1) {ll now = (M + m) / 2;if (now * now <= x) {m = now;}else {M = now;}}return m;}ll keta(ll n) {ll ret = 0;while (n) {n /= 10;ret++;}return ret;}ll ceil(ll n, ll m) {ll ret = n / m;if (n % m) ret++;return ret;}ll pow_ll(ll x, ll n) {if (n == 0) return 1;if (n % 2) {return pow_ll(x, n - 1) * x;}else {ll tmp = pow_ll(x, n / 2);return tmp * tmp;}}int main() {ll n, k;cin >> n >> k;ll now = 0;now = (now - pow_mod(2, k, n) + 1 + 2 * n) % n;mint dp[n];for (int i = 0; i < n; i++) dp[i] = 0;if (n % 2 == 0) {ll syuuki = n / 2;ll amari = pow_mod(2, k, syuuki);for (int i = 0; i < syuuki; i++) {dp[(now + i * 2) % n] = (pow_mod(2, k, MOD) - amari) * pow_mod(syuuki, MOD - 2, MOD);}for (int i = 0; i < amari; i++) {dp[(now + i * 2) % n] += 1;}}else {ll syuuki = n;ll amari = pow_mod(2, k, syuuki);for (int i = 0; i < syuuki; i++) {dp[(now + i * 2) % n] = (pow_mod(2, k, MOD) - amari) * pow_mod(syuuki, MOD - 2, MOD);}for (int i = 0; i < amari; i++) {dp[(now + i * 2) % n] += 1;}}for (int i = 0; i < n; i++) {mint ans = dp[i] * pow_mod(pow_mod(2, k, MOD), MOD - 2, MOD);cout << ans.val() << endl;}return 0;}