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

智能推荐

FX3/CX3 JLINK 调试_ezusbsuite_qsg.pdf-程序员宅基地

文章浏览阅读2.1k次。FX3 JLINK调试是一个有些麻烦的事情,经常有些莫名其妙的问题。 设置参见 c:\Program Files (x86)\Cypress\EZ-USB FX3 SDK\1.3\doc\firmware 下的 EzUsbSuite_UG.pdf 文档。 常见问题: 1.装了多个版本的jlink,使用了未注册或不适当的版本 选择一个正确的版本。JLinkARM_V408l,JLinkA_ezusbsuite_qsg.pdf

用openGL+QT简单实现二进制stl文件读取显示并通过鼠标旋转缩放_qopengl如何鼠标控制旋转-程序员宅基地

文章浏览阅读2.6k次。** 本文仅通过用openGL+QT简单实现二进制stl文件读取显示并通过鼠标旋转缩放, 是比较入门的级别,由于个人能力有限,新手级别,所以未能施加光影灯光等操作, 未能让显示的stl文件更加真实。****效果图:**1. main.cpp```cpp#include "widget.h"#include <QApplication>int main(int argc, char *argv[]){ QApplication a(argc, argv); _qopengl如何鼠标控制旋转

刘焕勇&王昊奋|ChatGPT对知识图谱的影响讨论实录-程序员宅基地

文章浏览阅读943次,点赞22次,收藏19次。以大规模预训练语言模型为基础的chatgpt成功出圈,在近几日已经给人工智能板块带来了多次涨停,这足够说明这一风口的到来。而作为曾经的风口“知识图谱”而言,如何找到其与chatgpt之间的区别,找好自身的定位显得尤为重要。形式化知识和参数化知识在表现形式上一直都是大家考虑的问题,两种技术都应该有自己的定位与价值所在。知识图谱构建往往是抽取式的,而且往往包含一系列知识冲突检测、消解过程,整个过程都能溯源。以这样的知识作为输入,能在相当程度上解决当前ChatGPT的事实谬误问题,并具有可解释性。

如何实现tomcat的热部署_tomcat热部署-程序员宅基地

文章浏览阅读1.3k次。最重要的一点,一定是degbug的方式启动,不然热部署不会生效,注意,注意!_tomcat热部署

用HTML5做一个个人网站,此文仅展示个人主页界面。内附源代码下载地址_个人主页源码-程序员宅基地

文章浏览阅读10w+次,点赞56次,收藏482次。html5 ,用css去修饰自己的个人主页代码如下:&lt;!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd"&gt;&lt;html xmlns="http://www.w3.org/1999/xh..._个人主页源码

程序员公开上班摸鱼神器!有了它,老板都不好意思打扰你!-程序员宅基地

文章浏览阅读201次。开发者(KaiFaX)面向全栈工程师的开发者专注于前端、Java/Python/Go/PHP的技术社区来源:开源最前线链接:https://github.com/svenstaro/gen..._程序员怎么上班摸鱼

随便推点

UG\NX二次开发 改变Block UI界面的尺寸_ug二次开发 调整 对话框大小-程序员宅基地

文章浏览阅读1.3k次。改变Block UI界面的尺寸_ug二次开发 调整 对话框大小

基于深度学习的股票预测(完整版,有代码)_基于深度学习的股票操纵识别研究python代码-程序员宅基地

文章浏览阅读1.3w次,点赞18次,收藏291次。基于深度学习的股票预测数据获取数据转换LSTM模型搭建训练模型预测结果数据获取采用tushare的数据接口(不知道tushare的筒子们自行百度一下,简而言之其免费提供各类金融数据 , 助力智能投资与创新型投资。)python可以直接使用pip安装tushare!pip install tushareCollecting tushare Downloading https://files.pythonhosted.org/packages/17/76/dc6784a1c07ec040e74_基于深度学习的股票操纵识别研究python代码

中科网威工业级防火墙通过电力行业测评_电力行业防火墙有哪些-程序员宅基地

文章浏览阅读2k次。【IT168 厂商动态】 近日,北京中科网威(NETPOWER)工业级防火墙通过了中国电力工业电力设备及仪表质量检验测试中心(厂站自动化及远动)测试,并成为中国首家通过电力协议访问控制专业测评的工业级防火墙生产厂商。   北京中科网威(NETPOWER)工业级防火墙专为工业及恶劣环境下的网络安全需求而设计,它采用了非X86的高可靠嵌入式处理器并采用无风扇设计,整机功耗不到22W,具备极_电力行业防火墙有哪些

第十三周 ——项目二 “二叉树排序树中查找的路径”-程序员宅基地

文章浏览阅读206次。/*烟台大学计算机学院 作者:董玉祥 完成日期: 2017 12 3 问题描述:二叉树排序树中查找的路径 */#include #include #define MaxSize 100typedef int KeyType; //定义关键字类型typedef char InfoType;typedef struct node

C语言基础 -- scanf函数的返回值及其应用_c语言ignoring return value-程序员宅基地

文章浏览阅读775次。当时老师一定会告诉你,这个一个"warning"的报警,可以不用管它,也确实如此。不过,这条报警信息我们至少可以知道一点,就是scanf函数调用完之后是有一个返回值的,下面我们就要对scanf返回值进行详细的讨论。并给出在编程时利用scanf的返回值可以实现的一些功能。_c语言ignoring return value

数字医疗时代的数据安全如何保障?_数字医疗服务保障方案-程序员宅基地

文章浏览阅读9.6k次。十四五规划下,数据安全成为国家、社会发展面临的重要议题,《数据安全法》《个人信息保护法》《关键信息基础设施安全保护条例》已陆续施行。如何做好“数据安全建设”是数字时代的必答题。_数字医疗服务保障方案

推荐文章

热门文章

相关标签