結果

問題 No.1419 Power Moves
ユーザー au7777au7777
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <bits/stdc++.h>
#include <atcoder/all>
typedef long long int ll;
typedef long long int ull;
#define MP make_pair
using 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;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0