結果

問題 No.573 a^2[i] = a[i]
ユーザー asugen0402
提出日時 2019-06-14 15:01:12
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 38 ms / 2,000 ms
コード長 3,666 bytes
コンパイル時間 245 ms
コンパイル使用メモリ 32,512 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-07 19:35:23
合計ジャッジ時間 1,823 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 47
権限があれば一括ダウンロードができます

ソースコード

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

#include <float.h>
#include <limits.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <time.h>
// 内部定数
#define D_MOD           1000000007                              // 除数(10の9乗+7)
#define D_FACT_MAX      100005                                  // 最大階乗数
// 内部変数
static FILE *szpFpI;                                            // 入力
static int si1ModFact[D_FACT_MAX];                              // 階乗リスト
static int si1ModFactR[D_FACT_MAX];                             // 階乗リスト(逆元)
// 内部変数 - テスト用
#ifdef D_TEST
    static int siRes;
    static FILE *szpFpA;
    static int siTNo;
#endif
// 出力
int
fOut(
    char *pcpLine               // <I> １行
)
{
    char lc1Buf[1024];
#ifdef D_TEST
    fgets(lc1Buf, sizeof(lc1Buf), szpFpA);
    if (strcmp(lc1Buf, pcpLine)) {
        siRes = -1;
    }
#else
    printf("%s", pcpLine);
#endif
    return 0;
}
// べき乗 - 取得
int
fGetModPower(
    int piBase                  // <I> 基数
    , int piIdx                 // <I> 指数
)
{
    // 平方リスト - 作成
    int li1Val[100];
    li1Val[0] = piBase;
    int liCnt = 1;
    int liIdx = 1;
    while (piIdx > liIdx) {
        li1Val[liCnt] = (int)((long long)li1Val[liCnt - 1] * (long long)li1Val[liCnt - 1] % (long long)D_MOD);
        liCnt++;
        liIdx += liIdx;
    }
    // べき乗 - 取得
    int liVal = 1;
    while (piIdx > 0) {
        if (piIdx >= liIdx) {
            piIdx -= liIdx;
            liVal = (int)((long long)liVal * (long long)li1Val[liCnt - 1] % (long long)D_MOD);
        }
        liCnt--;
        liIdx /= 2;
    }
    return liVal;
}
// 階乗リスト - 作成
int
fMakeModFact(
    int piMax                   // <I> 最大値
)
{
    int i;
    si1ModFact[0] = 1;
    si1ModFact[1] = 1;
    for (i = 2; i <= piMax; i++) {
        si1ModFact[i] = (int)((long long)si1ModFact[i - 1] * (long long)i % (long long)D_MOD);
    }
    return 0;
}
// 階乗リスト(逆元) - 作成
int
fMakeModFactR(
    int piMax                   // <I> 最大値
)
{
    int i;
    for (i = 0; i <= piMax; i++) {
        si1ModFactR[i] = fGetModPower(si1ModFact[i], D_MOD - 2);
    }
    return 0;
}
// nCk - 取得
int
fGetnCk(
    int piN                     // <I> N
    , int piK                   // <I> K
)
{
    int liCnt = (int)((long long)si1ModFact[piN] * (long long)si1ModFactR[piN - piK] % (long long)D_MOD);
    return (int)((long long)liCnt * (long long)si1ModFactR[piK] % (long long)D_MOD);
}
// 実行メイン
int
fMain(
)
{
    int i;
    char lc1Buf[1024];
    // N - 取得
    int liN;
    fgets(lc1Buf, sizeof(lc1Buf), szpFpI);
    sscanf(lc1Buf, "%d", &liN);
    // 階乗リスト - 作成
    fMakeModFact(liN);
    fMakeModFactR(liN);
    // 使用数字数でループ
    long long llPCnt = 0;
    for (i = 1; i <= liN; i++) {
        // 配列位置 - 使用数字
        long long llUCnt = fGetnCk(liN, i);
        // 配列位置 - 使用しない数字
        long long llOCnt = fGetModPower(i, liN - i);
        // 加算
        llPCnt += llUCnt * llOCnt;
        llPCnt %= D_MOD;
    }
    return (int)llPCnt;
}
// １回実行
int
fOne(
)
{
    int liRet;
    char lc1Buf[1024];
    // 入力 - セット
#ifdef D_TEST
    sprintf(lc1Buf, ".\\Test\\T%d.txt", siTNo);
    szpFpI = fopen(lc1Buf, "r");
    sprintf(lc1Buf, ".\\Test\\A%d.txt", siTNo);
    szpFpA = fopen(lc1Buf, "r");
    siRes = 0;
#else
    szpFpI = stdin;
#endif
    // 実行メイン
    liRet = fMain();
    // 出力
    sprintf(lc1Buf, "%d\n", liRet);
    fOut(lc1Buf);
    // 残データ有無
#ifdef D_TEST
    lc1Buf[0] = '\0';
    fgets(lc1Buf, sizeof(lc1Buf), szpFpA);
    if (strcmp(lc1Buf, "")) {
        siRes = -1;
    }
#endif
    // テストファイルクローズ
#ifdef D_TEST
    fclose(szpFpI);
    fclose(szpFpA);
#endif
    // テスト結果
#ifdef D_TEST
    if (siRes == 0) {
        printf("OK %d\n", siTNo);
    }
    else {
        printf("NG %d\n", siTNo);
    }
#endif
    return 0;
}
// プログラム開始
int
main()
{
#ifdef D_TEST
    int i;
    for (i = D_TEST_SNO; i <= D_TEST_ENO; i++) {
        siTNo = i;
        fOne();
    }
#else
    fOne();
#endif
    return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0