結果
問題 | No.766 金魚すくい |
ユーザー |
![]() |
提出日時 | 2018-12-14 00:41:01 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 72 ms / 1,500 ms |
コード長 | 1,303 bytes |
コンパイル時間 | 948 ms |
コンパイル使用メモリ | 100,180 KB |
実行使用メモリ | 8,192 KB |
最終ジャッジ日時 | 2024-09-25 04:50:31 |
合計ジャッジ時間 | 3,471 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 42 |
ソースコード
#include <cstdio>#include <cstring>#include <iostream>#include <string>#include <cmath>#include <bitset>#include <vector>#include <map>#include <set>#include <queue>#include <deque>#include <algorithm>#include <complex>#include <unordered_map>#include <unordered_set>#include <random>using namespace std;typedef long long int ll;typedef pair<int, int> P;const ll MOD=1e9+7;ll powmod(ll a, ll k){ll ap=a, ans=1;while(k>0){if(k%2==1){ans*=ap;ans%=MOD;}ap=ap*ap;ap%=MOD;k/=2;}return ans;}ll inv(ll a){return powmod(a, MOD-2);}int main(){int n, m; ll p;cin>>n>>m>>p;p=p*inv(100ll)%MOD;ll f[200001]; f[0]=1;for(ll i=1; i<=m+n; i++) f[i]=f[i-1]*i%MOD;ll invf[200001];invf[m+n]=inv(f[m+n]);for(ll i=m+n-1; i>=0; i--) invf[i]=invf[i+1]*(i+1)%MOD;ll s[100001]; s[0]=0;ll v[100001];for(int i=0; i<n; i++){cin>>v[i];}sort(v, v+n, greater<ll>());for(int i=0; i<n; i++) s[i+1]=(s[i]+v[i])%MOD;ll ans=0;ll ps=0;for(int i=0; i<=n-1; i++){ll q=f[i+m-1]*invf[i]%MOD*invf[m-1]%MOD*powmod((1-p+MOD)%MOD, (ll)i)%MOD*powmod(p, (ll)m)%MOD;ps=(ps+q)%MOD;ans+=(q*s[i]%MOD);ans%=MOD;}ans+=((1-ps+MOD)%MOD*s[n]%MOD);ans%=MOD;cout<<ans<<endl;return 0;}