BZOJ 2339: [HNOI2011]卡农动态规划+容斥原理

最新推荐文章于 2023-09-03 20:03:20 发布

EM-LGH

最新推荐文章于 2023-09-03 20:03:20 发布

阅读量137

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.csdn.net/liyong1009s/article/details/103867498

本文深入探讨了使用C++进行编程竞赛的高级技巧，包括快速幂运算、动态规划及组合数学的应用。通过具体代码示例，讲解了如何优化算法效率，解决复杂问题。

code:

#include <bits/stdc++.h>  
#define N 1000007 
#define ll long long 
#define mod 100000007 
#define setIO(s) freopen(s".in","r",stdin) 
using namespace std;
int qpow(int x,int y) 
{
    int tmp=1; 
    for(;y;y>>=1,x=(ll)x*x%mod) 
        if(y&1) tmp=(ll)tmp*x%mod; 
    return tmp; 
}
int A[N],f[N];  
int main()
{ 
    // setIO("input");  
    int i,j,n,m,inv; 
    scanf("%d%d",&n,&m);   
    inv=A[0]=1;   
    for(i=2;i<=m;++i) inv=(ll)i*inv%mod;   
    inv=qpow(inv,mod-2);   
    int mx=(qpow(2,n)-1+mod)%mod;     
    for(i=1;i<=m;++i) 
        A[i]=(ll)A[i-1]*(mx-i+1+mod)%mod;    
    f[0]=1,f[1]=0;   
    for(i=2;i<=m;++i) 
    {
        f[i]=(ll)(A[i-1]-f[i-1]+mod)%mod;  
        f[i]=(ll)(f[i]-(ll)f[i-2]*(i-1)%mod*(qpow(2,n)-1-(i-2)+mod)%mod+mod)%mod;  
    }
    printf("%d\n",(ll)f[m]*inv%mod);   
    return 0;
}