hdu 6712 sakura_av6712-程序员宅基地

技术标签: lucas与扩展lucas  

公式的话官方题解已经非常详细,这里就不再写公式了,大致推导为n步有x+y步是j,k两维移动,有n-x-y步是在i轴上移动。 在x+y的两维中,有y步是在y轴上移动,x步在x轴上移动。然后算上C(n,x+y)*C(x+y,x)*t1^(x/p)*t2(y/p)。就是每个点的贡献。这题卡常卡的太恶心了。

#include"bits/stdc++.h"
using namespace std;
typedef long long LL;
typedef long long ll;

const LL mod = 998244353;
LL pi = 2;
LL pe = 8388608, Pa[8388610];
LL per = pe-1;
LL fac[2][30],inv[2][30],mul0,mul1,mul2,qpe[1007];
inline LL qpow(LL a, LL b, LL m)
{
    LL ans = 1;
    a %= m;
    while(b > 0)
    {
        if(b&1) ans = ans*a%m;
        a = a*a%m;
        b >>= 1;
    }
    return ans;
}

inline LL fast_pow(LL a, LL b, LL m)
{
    LL ans = 1;
    LL mm = m-1;
    a &= mm;
    while(b > 0)
    {
        if(b&1) ans = ans*a&mm;
        a = a*a&mm;
        b >>= 1;
    }
    return ans;
}

inline LL exgcd(LL a,LL b,LL &x,LL &y) {
    if(b==0) {
        x=1,y=0;
        return a;
    }
    LL d=exgcd(b,a%b,y,x);
    y-=a/b*x;
    return d;
}
inline LL Inv(LL a,LL P){//求a在膜P下的逆
    if(a==0)return 1;
    LL x,y;
    exgcd(a,P,x,y);
    return (x%P+P)%P;
}


inline LL cal1(LL n) {
    LL ans=0;
    while(n > 0) {
        ans += (n>>1);
        n >>= 1;
    }
    return ans;
}

inline LL cal2(LL n) {
    if(n == 0) return 1;
    LL ans = 1;
    if(n/pe) ans = qpe[n/pe];
    ans = Pa[n&per];
    return cal2(n>>1)*ans&per;
}

inline LL cal(LL n,LL m, LL P) {
    LL cnt = cal1(n) - cal1(m) - cal1(n-m);
    LL a = cal2(n), b = cal2(m), c = cal2(n-m);
    LL ans =  ((a*Inv(b,pe)&per)*Inv(c,pe)&per)*fast_pow(pi,cnt,pe)&per;
    return ans;//范围可到达P*Pi,注意是否用快速乘
}


inline LL C(LL n, LL m, LL p, LL id)//组合数C(n, m) % p
{
    if(m > n) return 0;
    return fac[id][n] * inv[id][m]%p * inv[id][n - m] % p;
}
inline LL Lucas(LL n, LL m, LL p, LL id)
{
    if(m == 0) return 1;
    return C(n % p, m % p, p, id) * Lucas(n / p, m / p, p, id) % p;
}

inline LL solve(LL n,LL m,LL P) {
    if(m>n)return 0;
    LL ans = cal(n,m,P)*mul0%P;
    LL p = 7;
    LL t1 = Lucas(n,m,p,0)*mul1%P;
    p = 17;
    LL t2 = Lucas(n,m,p,1)*mul2%P;
    ans = (ans + t1 + t2)%P;
    return ans;
}

void init(LL P){
    mul0 = (P/pe)*Inv(P/pe,pe)%P;
    mul1 = (P/7)*Inv(P/7,7)%P;
    mul2 = (P/17)*Inv(P/17,17)%P;
    Pa[0] = Pa[1] = 1;
    for(LL j = 2; j <= pe; j++){
        if(j&1) Pa[j] = Pa[j-1]*j&per;
        else Pa[j] = Pa[j-1];
    }
    qpe[0] = Pa[pe];
    for(int i = 1; i <= 1000; i++)
        qpe[i] = qpe[i-1]*Pa[pe]%pe;
    LL p = 7;
    fac[0][0] = fac[0][1] = 1;
    inv[0][0] = inv[0][1] = 1;
    for(LL j = 2; j <= p; j++){
        fac[0][j] = fac[0][j-1]*j%p; 
        inv[0][j] = inv[0][j-1]*qpow(j,p-2,p)%p;
    }

    p = 17;
    fac[1][0] = fac[1][1] = 1;
    inv[1][0] = inv[1][1] = 1;
    for(LL j = 2; j <= p; j++){
        fac[1][j] = fac[1][j-1]*j%p; 
        inv[1][j] = inv[1][j-1]*qpow(j,p-2,p)%p;
    }
}

int main(){
#ifndef ONLINE_JUDGE
    freopen("in.txt","r",stdin);
#endif
    int t1,t2,p,q,n,m;
    init(mod-1);
    while(~scanf("%d%d%d%d%d%d",&t1,&t2,&p,&q,&n,&m)){
        LL ans = 1;
        for(LL i = 1; i <= m; i++){
            int ax,ay,av;
            scanf("%d%d%d",&ax,&ay,&av);

            if(ax%p != 0 || ay%q != 0){
                continue;
            }
            LL x = ax/p , y = ay/q;
            if(x+y > n) continue;
            LL tmp = solve(x+y, x, mod-1)*solve(n,x+y,mod-1)%(mod-1)*qpow(t1,x,mod-1)%(mod-1)*qpow(t2,y,mod-1)%(mod-1);
            ans = ans*qpow(av,tmp,mod)%mod;
        }
        printf("%lld\n",(ans%mod+mod)%mod);
    }
    return 0;
}

 

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/qq_18869763/article/details/100057852

智能推荐

软件测试流程包括哪些内容?测试方法有哪些?_测试过程管理中包含哪些过程-程序员宅基地

文章浏览阅读2.9k次,点赞8次,收藏14次。测试主要做什么?这完全都体现在测试流程中,同时测试流程是面试问题中出现频率最高的,这不仅是因为测试流程很重要,而是在面试过程中这短短的半小时到一个小时的时间,通过测试流程就可以判断出应聘者是否合适,故在测试流程中包含了测试工作的核心内容,例如需求分析,测试用例的设计,测试执行,缺陷等重要的过程。..._测试过程管理中包含哪些过程

政府数字化政务的人工智能与机器学习应用:如何提高政府工作效率-程序员宅基地

文章浏览阅读870次,点赞16次,收藏19次。1.背景介绍政府数字化政务是指政府利用数字技术、互联网、大数据、人工智能等新技术手段,对政府政务进行数字化改革,提高政府工作效率,提升政府服务质量的过程。随着人工智能(AI)和机器学习(ML)技术的快速发展,政府数字化政务中的人工智能与机器学习应用也逐渐成为政府改革的重要内容。政府数字化政务的人工智能与机器学习应用涉及多个领域,包括政策决策、政府服务、公共安全、社会治理等。在这些领域,人工...

ssm+mysql+微信小程序考研刷题平台_mysql刷题软件-程序员宅基地

文章浏览阅读219次,点赞2次,收藏4次。系统主要的用户为用户、管理员,他们的具体权限如下:用户:用户登录后可以对管理员上传的学习视频进行学习。用户可以选择题型进行练习。用户选择小程序提供的考研科目进行相关训练。用户可以进行水平测试,并且查看相关成绩用户可以进行错题集的整理管理员:管理员登录后可管理个人基本信息管理员登录后可管理个人基本信息管理员可以上传、发布考研的相关例题及其分析,并对题型进行管理管理员可以进行查看、搜索考研题目及错题情况。_mysql刷题软件

根据java代码描绘uml类图_Myeclipse8.5下JAVA代码导成UML类图-程序员宅基地

文章浏览阅读1.4k次。myelipse里有UML1和UML2两种方式,UML2功能更强大,但是两者生成过程差别不大1.建立Test工程,如下图,uml包存放uml类图package com.zz.domain;public class User {private int id;private String name;public int getId() {return id;}public void setId(int..._根据以下java代码画出类图

Flume自定义拦截器-程序员宅基地

文章浏览阅读174次。需求:一个topic包含很多个表信息,需要自动根据json字符串中的字段来写入到hive不同的表对应的路径中。发送到Kafka中的数据原本最外层原本没有pkDay和project,只有data和name。因为担心data里面会空值,所以根同事商量,让他们在最外层添加了project和pkDay字段。pkDay字段用于表的自动分区,proejct和name合起来用于自动拼接hive表的名称为 ..._flume拦截器自定义开发 kafka

java同时输入不同类型数据,Java Spring中同时访问多种不同数据库-程序员宅基地

文章浏览阅读380次。原标题:Java Spring中同时访问多种不同数据库 多样的工作要求,可以使用不同的工作方法,只要能获得结果,就不会徒劳。开发企业应用时我们常常遇到要同时访问多种不同数据库的问题,有时是必须把数据归档到某种数据仓库中,有时是要把数据变更推送到第三方数据库中。使用Spring框架时,使用单一数据库是非常容易的,但如果要同时访问多个数据库的话事件就变得复杂多了。本文以在Spring框架下开发一个Sp..._根据输入的不同连接不同的数据库

随便推点

EFT试验复位案例分析_eft电路图-程序员宅基地

文章浏览阅读3.6k次,点赞9次,收藏25次。本案例描述了晶振屏蔽以及开关电源变压器屏蔽对系统稳定工作的影响, 硬件设计时应考虑。_eft电路图

MR21更改价格_mr21 对于物料 zba89121 存在一个当前或未来标准价格-程序员宅基地

文章浏览阅读1.1k次。对于物料价格的更改,可以采取不同的手段:首先,我们来介绍MR21的方式。 需要说明的是,如果要对某一产品进行价格修改,必须满足的前提条件是: ■ 1、必须对价格生效的物料期间与对应会计期间进行开启; ■ 2、该产品在该物料期间未发生物料移动。执行MR21,例如更改物料1180051689的价格为20000元,系统提示“对于物料1180051689 存在一个当前或未来标准价格”,这是因为已经对该..._mr21 对于物料 zba89121 存在一个当前或未来标准价格

联想启天m420刷bios_联想启天M420台式机怎么装win7系统(完美解决usb)-程序员宅基地

文章浏览阅读7.4k次,点赞3次,收藏13次。[文章导读]联想启天M420是一款商用台式电脑,预装的是win10系统,用户还是喜欢win7系统,该台式机采用的intel 8代i5 8500CPU,在安装安装win7时有很多问题,在安装win7时要在BIOS中“关闭安全启动”和“开启兼容模式”,并且安装过程中usb不能使用,要采用联想win7新机型安装,且默认采用的uefi+gpt模式,要改成legacy+mbr引导,那么联想启天M420台式电..._启天m420刷bios

冗余数据一致性,到底如何保证?-程序员宅基地

文章浏览阅读2.7k次,点赞2次,收藏9次。一,为什么要冗余数据互联网数据量很大的业务场景,往往数据库需要进行水平切分来降低单库数据量。水平切分会有一个patition key,通过patition key的查询能..._保证冗余性

java 打包插件-程序员宅基地

文章浏览阅读88次。是时候闭环Java应用了 原创 2016-08-16 张开涛 你曾经因为部署/上线而痛苦吗?你曾经因为要去运维那改配置而烦恼吗?在我接触过的一些部署/上线方式中,曾碰到过以下一些问题:1、程序代码和依赖都是人工上传到服务器,不是通过工具进行部署和发布;2、目录结构没有规范,jar启动时通过-classpath任意指定;3、fat jar,把程序代码、配置文件和依赖jar都打包到一个jar中,改配置..._那么需要把上面的defaultjavatyperesolver类打包到插件中

VS2015,Microsoft Visual Studio 2005,SourceInsight4.0使用经验,Visual AssistX番茄助手的安装与基本使用9_番茄助手颜色-程序员宅基地

文章浏览阅读909次。1.得下载一个番茄插件,按alt+g才可以有函数跳转功能。2.不安装番茄插件,按F12也可以有跳转功能。3.进公司的VS工程是D:\sync\build\win路径,.sln才是打开工程的方式,一个是VS2005打开的,一个是VS2013打开的。4.公司库里的线程接口,在CmThreadManager.h 里,这个里面是我们的线程库,可以直接拿来用。CreateUserTaskThre..._番茄助手颜色

推荐文章

热门文章

相关标签