結果

問題 No.502 階乗を計算するだけ
ユーザー coderkc
提出日時 2019-06-06 01:50:55
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 997 bytes
コンパイル時間 715 ms
コンパイル使用メモリ 65,520 KB
実行使用メモリ 7,368 KB
最終ジャッジ日時 2024-09-22 16:39:58
合計ジャッジ時間 2,601 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 2 WA * 50
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:25:23: warning: ‘N’ is used uninitialized [-Wuninitialized]
   25 |     cout << Factrial[N] << endl;
      |                       ^

ソースコード

diff # 

#include <iostream>
#include <math.h>
#include <vector>
#define ll long long
using namespace std;

//x!を10^9+7で割った余り
ll Factrial[400010];
//(x!)^-1=(x!)^mod-2
ll FactInv[100010];
ll mod = 1000000007;
vector<int> v;

void makeFactrial();
void makeFactInv();
ll powmod(int a);


int main(){

    int N;
    
    makeFactrial();

    cout << Factrial[N] << endl;
    return 0;
}


void makeFactrial(){
    Factrial[0] = 1;
    for(int i=1;i<400000;i++){
        Factrial[i] = (i * Factrial[i-1]) % mod;
    }
}
void makeFactInv(){
    for(int i=0;i<100000;i++){
        FactInv[i] = powmod(Factrial[i]);
    }
}

ll powmod(int a){
    int ret=1;
    for(int i=mod;i>0;i -= (i&-i)){
        v.push_back(log2(i&-i));
    }
    int vmax = v[v.size()-1];
    int a_2_mod[vmax];
    a_2_mod[0] = a%mod;
    for(int i=1;i<=vmax;i++){
        a_2_mod[i] = (a_2_mod[i-1]*a_2_mod[i-1])%mod;
    }
    for(int i=0;i<v.size();i++){
        ret = (ret * a_2_mod[v[i]])%mod;
    }
    return ret;
}
0