1. 程式人生 > >2018牛客多校第九場E(動態規劃,思維,取模)

2018牛客多校第九場E(動態規劃,思維,取模)

pac namespace for ons mod 思維 space scan 動態規劃

#include<bits/stdc++.h>
using namespace std;
const long long mod=1000000007,inv=570000004;
long long i,j,n,m,a[1010]/*成功幾率*/,sum=0,dp[1010][1010]/*動態規劃*/,def[1010]//預處理;
long long qpow(long long x,long long y,long long mod)//快速冪
{
long long ans=1,tmp=x;
while(y)
{
if(y&1)
ans=ans*tmp%mod;
tmp=tmp*tmp%mod;
y>>=1;
}
return ans;
}
int main()
{
scanf("%lld%lld",&n,&m);
a[0]=0;
a[n+1]=0;
sum=0;
for(i=1;i<=n;i++)
{
scanf("%lld",&a[i]);
def[i]=qpow(i,m,mod);//預處理每段長度正確的分數
}
for(i=1;i<=n;i++)
{
dp[i][i]=a[i]*inv%mod;//單題成功率
for(j=i+1;j<=n;j++)
{
dp[i][j]=dp[i][j-1]*a[j]%mod*inv%mod;//動態規劃,i~j-1正確的基礎上j題正確
}
}
for(i=0;i<=n;i++)
{
for(j=i+2;j<=n+1;j++)
{
sum+=dp[i+1][j-1]*def[j-i-1]%mod*(100-a[i])%mod*inv%mod*(100-a[j])%mod*inv%mod;//i和j都不對而i+1~j-1正確
sum%=mod;
}
}
printf("%lld",sum);
return 0;
}
//動態規劃,思維,每進行一次運算都要進行取模

2018牛客多校第九場E(動態規劃,思維,取模)