結果
問題 | No.444 旨味の相乗効果 |
ユーザー | rickytheta |
提出日時 | 2016-11-11 23:38:35 |
言語 | C++11 (gcc 11.4.0) |
結果 |
AC
|
実行時間 | 181 ms / 2,500 ms |
コード長 | 1,583 bytes |
コンパイル時間 | 2,504 ms |
コンパイル使用メモリ | 151,308 KB |
実行使用メモリ | 4,380 KB |
最終ジャッジ日時 | 2023-09-07 02:10:50 |
合計ジャッジ時間 | 3,112 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge11 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
4,376 KB |
testcase_01 | AC | 2 ms
4,376 KB |
testcase_02 | AC | 1 ms
4,380 KB |
testcase_03 | AC | 1 ms
4,380 KB |
testcase_04 | AC | 181 ms
4,380 KB |
testcase_05 | AC | 2 ms
4,380 KB |
testcase_06 | AC | 1 ms
4,376 KB |
testcase_07 | AC | 8 ms
4,380 KB |
testcase_08 | AC | 3 ms
4,380 KB |
testcase_09 | AC | 3 ms
4,380 KB |
testcase_10 | AC | 1 ms
4,380 KB |
testcase_11 | AC | 2 ms
4,380 KB |
testcase_12 | AC | 92 ms
4,380 KB |
testcase_13 | AC | 47 ms
4,380 KB |
testcase_14 | AC | 2 ms
4,376 KB |
testcase_15 | AC | 2 ms
4,376 KB |
testcase_16 | AC | 180 ms
4,376 KB |
testcase_17 | AC | 2 ms
4,380 KB |
testcase_18 | AC | 2 ms
4,376 KB |
testcase_19 | AC | 19 ms
4,376 KB |
testcase_20 | AC | 141 ms
4,376 KB |
testcase_21 | AC | 2 ms
4,376 KB |
testcase_22 | AC | 2 ms
4,376 KB |
testcase_23 | AC | 1 ms
4,380 KB |
testcase_24 | AC | 2 ms
4,380 KB |
testcase_25 | AC | 2 ms
4,380 KB |
testcase_26 | AC | 1 ms
4,380 KB |
testcase_27 | AC | 2 ms
4,376 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; typedef long long ll; typedef vector<int> vi; typedef vector<ll> vl; typedef pair<int,int> pii; typedef pair<ll,ll> pll; typedef int _loop_int; #define REP(i,n) for(_loop_int i=0;i<(_loop_int)(n);++i) #define FOR(i,a,b) for(_loop_int i=(_loop_int)(a);i<(_loop_int)(b);++i) #define FORR(i,a,b) for(_loop_int i=(_loop_int)(b)-1;i>=(_loop_int)(a);--i) #define DEBUG(x) cout<<#x<<": "<<x<<endl #define DEBUG_VEC(v) cout<<#v<<":";REP(i,v.size())cout<<" "<<v[i];cout<<endl #define ALL(a) (a).begin(),(a).end() #define CHMIN(a,b) a=min((a),(b)) #define CHMAX(a,b) a=max((a),(b)) // mod const ll MOD = 1000000007ll; #define FIX(a) ((a)%MOD+MOD)%MOD // floating typedef double Real; const Real EPS = 1e-11; #define EQ0(x) (abs(x)<EPS) #define EQ(a,b) (abs(a-b)<EPS) typedef complex<Real> P; ll n,c; ll a[85]; typedef vector<vl> mat; mat mul(mat &a,mat &b){ int n = a.size(); mat ret(n,vl(n,0)); REP(k,n)REP(i,n)REP(j,n)ret[i][j] = (ret[i][j] + a[i][k]*b[k][j]%MOD)%MOD; return ret; } ll mypow(ll a,ll b){ b %= (MOD-1); ll r = 1; while(b){ if(b&1)r=r*a%MOD; a=a*a%MOD; b>>=1; } return r; } int main(){ scanf("%lld%lld",&n,&c); REP(i,n)scanf("%lld",a+i); mat A(n,vl(n,0)); REP(i,n)REP(j,i+1)A[i][j] = a[i]; mat I(n,vl(n,0)); REP(i,n)I[i][i]=1; mat B = I; ll k = c; --c; while(c){ if(c&1) B=mul(A,B); A=mul(A,A); c>>=1; } ll sum = 0; REP(i,n)REP(j,n)sum = (sum + B[i][j]*a[j]%MOD)%MOD; REP(i,n)sum = (sum-mypow(a[i],k)+MOD)%MOD; printf("%lld\n",sum); return 0; }