結果

問題 No.573 a^2[i] = a[i]
ユーザー Lepton_s
提出日時 2017-11-22 20:47:09
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 14 ms / 2,000 ms
コード長 1,889 bytes
コンパイル時間 783 ms
コンパイル使用メモリ 92,248 KB
実行使用メモリ 5,888 KB
最終ジャッジ日時 2024-11-26 11:22:57
合計ジャッジ時間 2,385 ms
ジャッジサーバーID
(参考情報)
judge2 / judge5
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 47
権限があれば一括ダウンロードができます

ソースコード

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

#include <algorithm>
#include <climits>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <ctime>
#include <iostream>
#include <sstream>
#include <functional>
#include <map>
#include <string>
#include <cstring>
#include <vector>
#include <queue>
#include <stack>
#include <deque>
#include <set>
#include <list>
#include <numeric>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<ll,ll> P;
const double PI = 3.14159265358979323846;
const double EPS = 1e-12;
const ll INF = 1LL<<29;
const ll mod = 1e9+7;
#define rep(i,n) for(int (i)=0;(i)<(ll)(n);++(i))
#define repd(i,n,d) for(ll (i)=0;(i)<(ll)(n);(i)+=(d))
#define all(v) (v).begin(), (v).end()
#define pb(x) push_back(x)
#define mp(x,y) make_pair((x),(y))
#define mset(m,v) memset((m),(v),sizeof(m))
#define chmin(X,Y) ((X)>(Y)?X=(Y),true:false)
#define chmax(X,Y) ((X)<(Y)?X=(Y),true:false)
#define fst first
#define snd second
#define UNIQUE(x) (x).erase(unique(all(x)),(x).end())
template<class T> ostream &operator<<(ostream &os, const vector<T> &v){int n=v.size();rep(i,n)os<<v[i]<<(i==n-1?"":" ");return os;}
#define N 100010
ll inv[N], fact[N], ifact[N];
void init_fact(ll n = N){
inv[1] = 1;
for(int i = 2; i < n;i++) inv[i] = inv[mod%i] * (mod - mod/i) % mod;
fact[0] = ifact[0] = 1;
for(int i = 1; i < n; i++){
fact[i] = (fact[i-1]*i)%mod;
ifact[i]=(ifact[i-1]*inv[i])%mod;
}
}
ll comb(ll a, ll b){
return fact[a+b]*ifact[a]%mod*ifact[b]%mod;
}
ll comb_nk(ll n, ll k){
return comb(n-k, k);
}
ll pow_mod(ll a, ll r, ll m){
ll x = 1;
while(r){
if(r&1) (x*=a)%=m;
(a*=a)%=m;
r>>=1;
}
return x;
}
int main(){
ll x;
cin>>x;
init_fact();
ll res = 0;
rep(i, x){
res += comb(x-i, i)*pow_mod(x-i, i, mod)%mod;
res %= mod;
}
cout<<res<<endl;
return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0