結果
問題 | No.2351 Butterfly in Summer |
ユーザー |
|
提出日時 | 2023-07-27 16:00:07 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 5,416 bytes |
コンパイル時間 | 3,315 ms |
コンパイル使用メモリ | 251,180 KB |
最終ジャッジ日時 | 2025-02-15 19:33:21 |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 21 |
ソースコード
#pragma region Macros#include <bits/stdc++.h>#include <atcoder/all>#if defined(LOCAL) || defined(_DEBUG)#include "debug.hpp"#else#define O(...)#define START()#define STOP()#define MEMORY()#endifusing namespace std;#define REP(i, n) for(int i=0, i##_len=(n); i<i##_len; ++i)#define REPR(i, n) for(int i=(n); i>=0; --i)#define FOR(i, n, m) for(int i=(m), i##_len=(n); i<i##_len; ++i)#define EACH(i, v) for(const auto& i : v)#define ALL(x) (x).begin(),(x).end()#define ALLR(x) (x).rbegin(),(x).rend()template<class T, class U>bool chmax(T &a, const U &b) { if (a<(T)b) { a=(T)b; return 1; } return 0; }template<class T, class U>bool chmin(T &a, const U &b) { if (b<(T)a) { a=(T)b; return 1; } return 0; }#define vec vector#define umap unordered_map#define uset unordered_setusing ll = long long;using ld = long double;using P = pair<ll, ll>;using Tup = tuple<ll, ll, ll>;using vl = vec<ll>;#define fi first#define se second#define el endlconstexpr ll INF = numeric_limits<ll>::max()/2-1;template<class T> istream &operator>>(istream &stream, vec<T>& o){REP(i, o.size())stream >> o[i];return stream;}template<class T, class U> istream &operator>>(istream &stream, pair<T, U>& o){cin >> o.fi >> o.se; return stream;}namespace myinput {void input() {}template<class T, class... Ts> void input(T&& o, Ts&&... args){cin >> o;input(forward<Ts>(args)...);}void assign_vl(size_t siz) {};template<class T, class... Ts> void assign_vl(size_t siz, T&& o, Ts&&... args){o.resize(siz);assign_vl(siz, forward<Ts>(args)...);}void input_horizon_sub(size_t index) {};template<class T, class... Ts> void input_horizon_sub(size_t index, T&& o, Ts&&... args) {cin>>o[index];input_horizon_sub(index, forward<Ts>(args)...);}template<class... Ts> void input_horizon(size_t siz, Ts&&... args){REP(i, siz) input_horizon_sub(i, forward<Ts>(args)...);}}#define _I(T, ...) ;myinput::input(__VA_ARGS__);#define I(T, ...) ;T __VA_ARGS__;_I(T, __VA_ARGS__);#define _Iv(T, siz, ...) ;myinput::assign_vl(siz, __VA_ARGS__);myinput::input(__VA_ARGS__);#define Iv(T, siz, ...) ;vec<T> __VA_ARGS__;_Iv(T, siz, __VA_ARGS__);#define _Ih(T, siz, ...) ;myinput::assign_vl(siz, __VA_ARGS__);myinput::input_horizon(siz, __VA_ARGS__);#define Ih(T, siz, ...) ;vec<T> __VA_ARGS__;_Ih(T, siz, __VA_ARGS__);#pragma endregionvoid Main();int main(){std::cin.tie(nullptr);std::cout << std::fixed << std::setprecision(15);Main();MEMORY();return 0;}// mod Pのための型#pragma region modinttemplate<ll m> struct modint{protected:ll v=0;public://コンストラクタmodint(ll x=0){while(x<0)x+=m;v=x%m;}modint(const modint &r){v=r.val();}//算術演算子(単項)modint operator -(){return modint(-v);}modint& operator ++(){if(++v>=m)v-=m; return *this;}modint operator ++(int){modint t = *this; ++*this; return t;}modint& operator --(){if(--v<0)v+=m; return *this;}modint operator --(int){modint t = *this; --*this; return t;}//算術演算子(二項)friend modint operator +(const modint &l, const modint &r){return modint(l)+=r;}friend modint operator -(const modint &l, const modint &r){return modint(l)-=r;}friend modint operator *(const modint &l, const modint &r){return modint(l)*=r;}friend modint operator /(const modint &l, const modint &r){return modint(l)/=r;}// キャスト演算子explicit operator int() const {return val();}explicit operator ll() const {return val();}//代入演算子modint &operator +=(const modint &r){v+=r.val();if(v>=m)v-=m;return *this;}modint &operator -=(const modint &r){if(v<r.val())v+=m;v-=r.val();return *this;}modint &operator *=(const modint &r){v=v*r.val()%m;return *this;}modint &operator /=(const modint &r){v=v*r.inv().val()%m;return *this;}//等価比較演算子friend bool operator ==(const modint &l, const modint& r){return l.val()==r.val();}friend bool operator <(const modint &l, const modint& r){return l.val()<r.val();}friend bool operator >(const modint &l, const modint& r){return l.val()>r.val();}friend bool operator !=(const modint &l, const modint& r){return l.val()!=r.val();}// 各種関数int val() const{return v;}constexpr static ll mod(){return m;}static modint raw(ll x){modint r;r.v=x;return r;}modint inv() const{ll a=v,b=m,u=1,v=0;while(b){ll t=a/b;a-=t*b;swap(a,b);u-=t*v;swap(u,v);}u %= m;if (u < 0) u += m;return u;}modint pow(ll n) const{modint ret = 1, x(v);while (n > 0) {if (n & 1) ret *= x;x *= x;n >>= 1;}return ret;}//入出力ストリームfriend istream &operator>>(istream &is,modint& x){ll t;is>>t;x=t;return (is);}friend ostream &operator<<(ostream &os,const modint& x){return os<<x.val();}};using modint1000000007 = modint<1000000007>;using modint998244353 = modint<998244353>;#pragma endregionusing mint = modint998244353;void Main(){I(ll, n, k);mint x = (k * (k-1)) / 2;x *= n;x *= 2;mint y = mint(k).pow(n);O(x, y);cout << x/y << el;}