JZOJ.4709【NOIP2016提高A组模拟8.17】Matrix

最新推荐文章于 2016-08-17 22:14:23 发布

Luckfort

最新推荐文章于 2016-08-17 22:14:23 发布

阅读量356

点赞数 2

本文链接：https://blog.csdn.net/huangjingyuan107/article/details/52234288

版权

Problem

Description

这里写图片描述

Input

这里写图片描述

Output

这里写图片描述

Sample Input

4 3 5
4 1 7 3
4 7 4 8

Sample Output

59716

Data Constraint

这里写图片描述

Solution

不要看到 $N≤10^5$ 就害怕，其实这题很简单的。
对于a=0，答案就是 $T_n*b^{n-1}$ .
我们可以发现，只有L数组和T数组控制着答案。每次向右就乘以a，向下就乘以b，所以我们可以得到答案是

Σ n i = 2 C n - i 2 n - i - 2 * (a n - 1 b n - i * l [i] + a n - i b n - 1 * t [i])

$\Sigma_{i=2}^n C_{2n-i-2}^{n-i}*(a^{n-1}b^{n-i}*l[i]+a^{n-i}b^{n-1}*t[i])$
可是这样还是超时啊!!!!!!!!!!!!
那么我们可以将阶乘、逆元预处理啊……………………

Code

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#define N 100010
#define LL long long
#define mo 1000000007
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
int a,b,i,j,n,m;
LL l[N],t[N],ans,jc[N*2],jj[N*2];
LL ksm(LL x,LL y)
{
    LL s=1,t=x;
    while (y)
    {
        if (y%2==1) s=(s*t)%mo;
        t=(t*t)%mo;
        y/=2;
    }
    return s;
}
int main()
{
    scanf("%d%d%d",&n,&a,&b);
    fo(i,1,n) scanf("%lld",&l[i]);
    fo(i,1,n) scanf("%lld",&t[i]);
    jj[0]=jj[1]=jc[0]=jc[1]=1;
    fo(i,2,n*2)
    {
        jc[i]=(jc[i-1]*i)%mo;
        jj[i]=ksm(jc[i],mo-2);
    }
    ans=0;
    fo(i,2,n)
    {
        LL t1,t2,t3;
        t1=(((jc[2*n-i-2]*jj[n-i])%mo)*jj[n-2])%mo;
        t2=((ksm(a,n-1)*ksm(b,n-i)%mo)*l[i])%mo;
        t3=((ksm(a,n-i)*ksm(b,n-1)%mo)*t[i])%mo;
        ans=(ans+t1*((t2+t3)%mo))%mo;
    }
    printf("%lld",ans);
}