NOIP2016提高A组模拟8.17(雅礼联考day1)Matrix

Posted 无尽的蓝黄

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了NOIP2016提高A组模拟8.17(雅礼联考day1)Matrix相关的知识,希望对你有一定的参考价值。

题目

技术分享图片

分析

假设,我们从\(F_{i,2}\)出发,那么对\(F_{n,n}\)的贡献就是\(某个系数乘以a^{n-i}b^{n-1}r_i\)
同理,如果从\(F_{2,i}\)出发,那么对\(F_{n,n}\)的贡献就是\(某个系数乘以a^{n-1}b^{n-i}l_i\)
那么某个系数是什么呢?
感性理解一下,就是从出发点到(n,n)的方案数\(C_{2*n-i-2}^{n-i}\)
那么答案就是\[\sum_{i-2}^{n}C_{2*n-i-2}^{n-i}(a^{n-1}b^{n-i}l_i+a^{n-i}b^{n-1}r_i)\]

#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const long long maxlongint=2147483647;
const long long mo=1000000007;
const long long N=100005;
using namespace std;
long long l[N],r[N],a[N],b[N],jc[N*2],n,ny[N*2];
long long mi(long long x,long long y)
{
    long long sum=1;
    while(y)
    {
        if(y&1) sum=sum*x%mo;
        x=x*x%mo;
        y/=2;
    }
    return sum;
}
int main()
{
    scanf("%lld%lld%lld",&n,&a[1],&b[1]);
    for(long long i=1;i<=n;i++) scanf("%lld",&l[i]);
    for(long long i=1;i<=n;i++) scanf("%lld",&r[i]);
    jc[0]=1;
    for(long long i=1;i<=n*2;i++)
    {
        jc[i]=jc[i-1]*i%mo;
    }
    a[0]=b[0]=1;
    for(long long i=2;i<=n;i++)
    {
        a[i]=a[i-1]*a[1]%mo;
        b[i]=b[i-1]*b[1]%mo;
    }
    long long ans=0;
    for(int i=2;i<=n;i++)
    {
        int sum=jc[2*n-i-2]*mi(jc[n-i],mo-2)%mo*mi(jc[n-2],mo-2)%mo;
        sum=sum*(a[n-1]*b[n-i]%mo*l[i]%mo+a[n-i]*b[n-1]%mo*r[i]%mo)%mo;
        ans=(ans+sum)%mo;
    }
    printf("%lld\n",ans);
}

以上是关于NOIP2016提高A组模拟8.17(雅礼联考day1)Matrix的主要内容,如果未能解决你的问题,请参考以下文章

NOIP2016提高A组模拟8.17(雅礼联考day1)Value

2017.11.25NOIP提高组模拟赛A组

2017.12.09NOIP提高组模拟赛A组

2017.12.02NOIP提高组模拟赛A组

Jzoj4742NOIP2016提高A组模拟9.2快速幂单峰

NOIP2016提高A组五校联考4ksum