結果

問題 No.444 旨味の相乗効果
コンテスト
ユーザー vjudge1
提出日時 2024-12-24 21:09:43
言語 C++23
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 132 ms / 2,500 ms
コード長 1,400 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 3,475 ms
コンパイル使用メモリ 253,760 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-12-24 21:09:48
合計ジャッジ時間 4,897 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 5
other AC * 23
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>
#define int long long
#define all(v) v.begin(),v.end()
#define rall(v) v.rbegin(),v.rend()
using namespace std;
template<typename T>istream&operator>>(istream&I,vector<T>&v){for(auto&i:v)I>>i;return I;}
template<typename T>ostream&operator<<(ostream&O,vector<T>&v){for(auto&i:v)O<<i<<' ';return O;}
namespace AC{
const int mod=1e9+7;
struct matrix{
	vector<vector<int>>m;
	matrix(int N,int M):m(N,vector<int>(M,0)){}
	vector<int>&operator[](int x){return m[x];}
	matrix operator*(matrix&a){
		matrix res(m.size(),a.m[0].size());
		for(int i=0;i<m.size();i++)for(int j=0;j<a.m[0].size();j++)for(int k=0;k<m[0].size();k++)(res[i][j]+=m[i][k]*a[k][j])%=mod;
		return res;
	}
	matrix operator*=(matrix&a){return *this=*this*a;}
};
void pow(matrix&res,matrix a,int n){
	while(n){
		if(n&1)res*=a;
		a*=a;
		n>>=1;
	}
}
int pow(int a,int n){
	int res=1;
	while(n){
		if(n&1)res=res*a%mod;
		a=a*a%mod;
		n>>=1;
	}
	return res;
}
void solve(){
	int n,c;
	cin>>n>>c;
	vector<int>v(n);
	cin>>v;
	matrix m1(1,n),m2(n,n);
	for(int i=0;i<n;i++)m1[0][i]=v[i];
	for(int i=0;i<n;i++)for(int j=i;j<n;j++)m2[i][j]=v[j];
	pow(m1,m2,c-1);
	int ans=0;
	for(int i=0;i<n;i++)ans=(ans+m1[0][i])%mod;
	for(int i=0;i<n;i++)ans=(ans+mod-pow(v[i],c))%mod;
	cout<<ans;
}
}
signed main(){
	/*freopen("spice.in","r",stdin);
	freopen("spice.out","w",stdout);*/
	int t=1;
	//cin>>t;
	while(t--)AC::solve();
}
0