結果

問題 No.1049 Zero (Exhaust)
ユーザー walkrewalkre
提出日時 2020-05-08 22:06:59
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 27 ms / 2,000 ms
コード長 3,518 bytes
コンパイル時間 1,386 ms
コンパイル使用メモリ 170,944 KB
実行使用メモリ 18,900 KB
最終ジャッジ日時 2024-07-04 00:45:17
合計ジャッジ時間 2,511 ms
ジャッジサーバーID
(参考情報)
judge3 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 22
権限があれば一括ダウンロードができます

ソースコード

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

#include <bits/stdc++.h>
using namespace std;
using i64 = int64_t;
#define rep(i, x, y) for (i64 i = i64(x), i##_max_for_repmacro = i64(y); i < i##_max_for_repmacro; ++i)
#define debug(x) #x << "=" << (x)
#ifdef DEBUG
#define _GLIBCXX_DEBUG
#define print(x) std::cerr << debug(x) << " (L:" << __LINE__ << ")" << std::endl
#else
#define print(x)
#endif
template <i64 p>
class fp {
public:
i64 x;
fp() : x(0) {}
fp(i64 x_) : x((x_ % p + p) % p) {}
fp operator+() const { return fp(x); }
fp operator-() const { return fp(-x); }
fp& operator+=(const fp& y) {
x += y.x;
if (x >= p) x -= p;
return *this;
}
fp& operator-=(const fp& y) { return *this += -y; }
fp& operator*=(const fp& y) {
x = x * y.x % p;
return *this;
}
fp& operator/=(const fp& y) { return *this *= fp(inverse(y.x)); }
fp operator+(const fp& y) const { return fp(x) += y; }
fp operator-(const fp& y) const { return fp(x) -= y; }
fp operator*(const fp& y) const { return fp(x) *= y; }
fp operator/(const fp& y) const { return fp(x) /= y; }
bool operator==(const fp& y) const { return x == y.x; }
bool operator!=(const fp& y) const { return !(*this == y); }
i64 extgcd(i64 a, i64 b, i64& x, i64& y) {
i64 d = a;
if (b != 0) {
d = extgcd(b, a % b, y, x);
y -= (a / b) * x;
} else {
x = 1;
y = 0;
}
return d;
}
i64 inverse(i64 a) {
i64 x, y;
extgcd(a, p, x, y);
return (x % p + p) % p;
}
};
template <i64 p>
i64 abs(const fp<p>& x) { return x.x; }
template <i64 p>
istream& operator>>(istream& is, fp<p>& x) {
is >> x.x;
return is;
}
template <i64 p>
ostream& operator<<(ostream& os, const fp<p>& x) {
os << x.x;
return os;
}
template <typename T, typename U>
ostream& operator<<(ostream& os, const pair<T, U>& p) {
os << "(" << p.first << ", " << p.second << ")";
return os;
}
template <typename T>
ostream& operator<<(ostream& os, const vector<T>& vec) {
os << "[";
for (const auto& v : vec) {
os << v << ",";
}
os << "]";
return os;
}
template <typename T>
bool chmin(T& a, const T& b) {
if (a > b) {
a = b;
return true;
}
return false;
}
template <typename T>
bool chmax(T& a, const T& b) {
if (a < b) {
a = b;
return true;
}
return false;
}
template <typename A, typename T, size_t size>
void fill(A (&ary)[size], const T& val) {
fill((T*)ary, (T*)(ary + size), val);
}
constexpr int inf = 1.01e9;
constexpr i64 inf64 = 4.01e18;
constexpr long double eps = 1e-9;
// double(64bit)n(250). long double(80 x87 ?),
    2, .
constexpr i64 max_loop = 100;
void solve() {
constexpr i64 mod = 1'000'000'007;
i64 P,K;
cin >> P >> K;
vector<fp<mod>> all(K),cnt(K);
all[0]=2*P;
cnt[0]=1+P;
rep(i,1,K){
all[i]=all[i-1]*(2*P);
cnt[i]+=cnt[i-1]*(1+P);
cnt[i]+=(all[i-1]-cnt[i-1])*(1+1);
}
cout << cnt[K-1] << endl;
}
int main() {
std::cin.tie(0);
std::ios::sync_with_stdio(false);
cout.setf(ios::fixed);
cout.precision(16);
solve();
return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0