技术标签: 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;
}
1计算阶乘的和v2.0(4分)题目内容:假设有这样一个三位数m,其百位、十位和个位数字分别是a、b、c,如果m= a!+b!+c!,则这个三位数就称为三位阶乘和数(约定0!=1)。请编程计算并输出所有的三位阶乘和数。函数原型:?long Fact(int n);函数功能:计算n的阶乘输入格式:?无输出格式:"%d\n"为避免出现格式错误,请直接拷贝粘贴题目中给的格式字符串和提示信息到你的程序中。时..._c语言小孩吃梨问题
第1章 Java快速入门1.1 Java历史简介1.2 Java语言与面向对象的程序设计1.3 Java程序概述1.4 本章小结br>第2章 Java语言基础2.1 Java程序的结构2.2 Java程序的一些特殊语句2.3 变量、常量及数据类型2.4 Java标识符和关键字2.5 运算符和表达式2.6 本章小结br>第3章 流程控制3.1 算法简介3.2 选择结构第1章 Java快速..._java程序设计案例
#include<bits/stdc++.h>using namespace std;int w,h,ans=1;char mp[200][200];bool vis[200][200];int dir[4][2]={{1,0},{-1,0},{0,1},{0,-1}};struct node{ int x, y;};void bfs(node s){ vis[s.x][s.y]=1; queue<node> q; q.push(s); while(!._1216:红与黑
1、查看第一个例子:viewer.entities.add({position : Cesium.Cartesian3.fromDegrees(-75.59777, 40.03883),point : {pixelSize : 10,color : Cesium.Color.YELLOW}});可以看到entity通过viewer中的entities加载到场景中,entities是entity的集...
1 创建守护进程1.让init进程成为新产生进程的父进程。调用fork函数创建子进程后,使父进程立即退出。这样,产生的子进程将变成孤儿进程,并被init进程接管,同时,所产生的新进程将变为在后台运行。2.调用setsid()使得新创建的进程脱离控制终端,同时创建新的进程组,并成为该进程组的首进程。进程组 & 会话 & 控制终端进程组是一个或多个进程的集合,进程组ID是由领头进程..._syslogd 记录用户态
作者:寒小阳&&龙心尘 时间:2016年2月。 出处:http://blog.csdn.net/han_xiaoyang/article/details/50759472 http://blog.csdn.net/longxinchen_ml/article/details/50749614 声明:版权所有,转载请联系作者并注明出处1. 引言也许你和这个叫『机器学习』的家伙一点也不..._机器学习路线
出色的Java工具有助于提高工作效率。Java IDE 工具提供了多种用户独特需求和个人偏好来创建编程环境的方法。今天,播妞给大家分享8个程序员最爱的Java开源IDE工具,没有用过的小伙...
16进制与10进制的对应关系是:0-9对应0-9;A-F对应10-15。十六进制(英文名称:Hexadecimal),是计算机中数据的一种表示方法。同我们日常生活中的表示法不一样。它由0-9,A-F组成,字母不区分大小写。扩展资料:进制转换的理论:1、 二进制数、十六进制数转换为十进制数:用按权展开法把一个任意R 进制数a n a n-1 ...a1a 0 . a -1 a -2...a -m转换..._进制代表的字母
### 问题描述编译安装make时错误```bashIn file included from /root/swoole-4.5.9/src/core/base.cc:17:/root/swoole-4.5.9/include/swoole.h:533:26: error: ‘clock_id_t’ was not declared in this scopeint swoole_clock_ge..._note: suggested alternative:
当类间继承关系很简单时,super()的使用很简单。class A(object):def __init__(self):print('a')class B(A):def __init__(self):super(B, self).__init__()print('b')b = B()输出结果:ab当一个类继承多个类时,问题就复杂起来了,请看下例:class A(object):def __ini..._python多继承同名函数
前面的文章有讲过前端和后端分开部署的文章前端分离部署https://blog.csdn.net/qiumen/article/details/111661326后端分离部署https://blog.csdn.net/qiumen/article/details/111363308当前用的jboss版本是wildfly-18.0.1.Final项目背景maven管理springboot项目,创建了frontend和backend后端两个模块。前端是vue,后端是spri.._eclipse vue和springboot整合打成war
05年的关于P2P技术的综述,还不错,比较全面。是中科院计算机技术研究所的文章。第 1 章 Peer-To-Peer 介绍罗杰文中科院计算技术研究所 最近几年,Peer-to-Peer (对等计算,简称P2P) 迅速成为计算机界关注的热门话题之一,财富杂志更将P2P列为影响Internet未来的四项科技之一。“Peer”在英语里有“对等者”和“伙