結果
| 問題 |
No.978 Fibonacci Convolution Easy
|
| ユーザー |
|
| 提出日時 | 2021-09-06 00:00:10 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 83 ms / 2,000 ms |
| コード長 | 909 bytes |
| コンパイル時間 | 1,737 ms |
| コンパイル使用メモリ | 170,184 KB |
| 実行使用メモリ | 19,608 KB |
| 最終ジャッジ日時 | 2024-12-22 19:02:20 |
| 合計ジャッジ時間 | 3,422 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 21 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
#define all(A) A.begin(),A.end()
using vll = vector<ll>;
#define rep(i, n) for (long long i = 0; i < (long long)(n); i++)
using Graph = vector<vector<ll>>;
ll mod=1e9+7;
vll DP(1e5+7,-1);
ll DF(ll M){
if(DP[M]!=-1)return DP[M];
ll res=0;
for(ll i=1;i*i<=M;i++){
if(M%i==0){
res+=DF(M/i-1);
res%=mod;
if(M!=i*i){
res+=DF(i-1);
}
}
}
DP[M]=res;
return res;
}
int main() {
ll N,P;
cin>>N>>P;
vll A(N+3);
A[0]=0,A[1]=1;
ll S=0,D=0;
rep(i,N-2){
A[i+2]=A[i+1]*P+A[i];
A[i+2]%=mod;
}
rep(i,N){
S+=A[i];
D+=A[i]*A[i];
S%=mod;
D%=mod;
}
ll AN=S*S;
AN%=mod;
AN+=D;
AN+=mod;
AN%=mod;
AN*=(mod+1)/2;
AN%=mod;
cout<<AN<<endl;
}