結果

問題 No.978 Fibonacci Convolution Easy
ユーザー tarattata1tarattata1
提出日時 2019-07-11 17:50:36
言語 C++11
(gcc 13.3.0)
結果
AC  
実行時間 24 ms / 2,000 ms
コード長 1,213 bytes
コンパイル時間 905 ms
コンパイル使用メモリ 71,244 KB
実行使用メモリ 18,432 KB
最終ジャッジ日時 2024-09-18 20:50:11
合計ジャッジ時間 1,580 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 21
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main(int, char**)’:
main.cpp:41:10: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   41 |     scanf("%d%d", &n, &p);
      |     ~~~~~^~~~~~~~~~~~~~~~

ソースコード

diff # 

#include <stdio.h>
#include <string>
#include <cstring>
#include <stdlib.h>
#include <math.h>
#include <algorithm>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <list>
#include <iterator>
#include <assert.h>
#pragma warning(disable:4996) 

typedef long long ll;
#define MIN(a, b) ((a)>(b)? (b): (a))
#define MAX(a, b) ((a)<(b)? (b): (a))
#define LINF 9223300000000000000
#define INF 2140000000
const long long MOD = 1000000007;
using namespace std;

ll mpow(ll x, ll n){ //x^n(mod M)
    ll ans = 1;
    while(n != 0){
        if(n&1) ans = ans*x % MOD;
        x = x*x % MOD;
        n = n >> 1;
    }
    return ans;
}

ll minv(ll x){
    return mpow( x, MOD-2 );
}

int main(int argc, char* argv[])
{
    int n,p;
    scanf("%d%d", &n, &p);

    ll sum1=0, sum2=0;
    vector<ll> a(n+1);
    int i;
    for(i=1; i<=n; i++) {
        if(i==1) {
            a[i]=0;
        }
        else if(i==2) {
            a[i]=1;
        }
        else {
            a[i]=(p*a[i-1]+a[i-2])%MOD;
        }
        sum1=(sum1+a[i])%MOD;
        sum2=(sum2+a[i]*a[i])%MOD;
    }
    ll ans = (sum1*sum1 + sum2 + MOD)%MOD;
    ans = ans * minv(2) % MOD;
    printf("%lld\n", ans);

    return 0;
}
0