技术标签: 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;
}
文章浏览阅读128次。西班牙已经发布了基于Firefox OS的手机,但是不是任何人能体验到,很都人都在苦苦等待,现在我们有另一种方法,不需要任何命令,没有纷繁复杂的安装步骤,在该文中,会教大家在ubuntu系统中一种轻松体验 Firefox OS。只需在火狐浏览器中就可是实现:firefoxosvs4火狐OS模拟器是为开发者使用的,用来测试他们的程序的的火狐浏览..._基于ubuntu的firefox os
文章浏览阅读755次。本文出自 “至简李云” 博客,http://yunli.blog.51cto.com/831344/168865,作者:李云,转载时请务必以超链接形式标明文章 原始出处 、作者信息和本声明 单元测试(Unit Test, UT)是一个老生常谈的话题,在对这篇文章进行博客归类时,我还是将其归类为开发技术,尽管其带有测试两个字。如何做单元测试不是我这里想说的,而是业界对其认识的认识及重视是我想指出的。_单元测试有利于重构
文章浏览阅读949次,点赞2次,收藏7次。聚合函数Aggregations一、简单聚合 1.1 数据准备 1.2 count 1.3 countDistinct 1._sparksql 聚合函数
文章浏览阅读308次。Mysql 想要远程客户端链接,必须给root修改可以远程访问的权限一、在连接服务器后,操作mysql系统数据库mysql -u root -puse mysql;查询用户表命令:select User,authentication_string,Host from user这里也可以看出host默认都是localhost访问权限二、接下来就是最重要的部分了:..._ubuntu 连接远程mysql
文章浏览阅读1.3w次,点赞38次,收藏70次。graylog使用 graylog报警配置 graylog接springboot graylog使用总结 graylog邮件通知_graylog
文章浏览阅读1.6k次。使用IP摄像头,需要在手机上下载一个网络摄像头app,之后在app上开启云服务,就可以把手机摄像头当做电脑的另外一个摄像头。下载支持网络摄像头的app首先下载任意款网络摄像头,我的手机是华为mate20,我下载的软件是IP摄像头。打开app后,点击开启云服务就会出现以下内容了这里出现192.168..:8080等字样,把手机作为摄像头。点击开启服务器运行如下代码:拿到192..._opencv打开ip摄像头 [tcp @ 0x5586bdf2cbc0] connection to
文章浏览阅读466次。DOCTYPE html>富文本编辑器实现原理title>#edit{height:260px;width:100%;overflow:scroll;border:solid 1px black}style>head>div>居中button>左对齐button>右对齐button>添加缩进button>去掉缩进button>宋体button>大字体button>红色字体button>..._html标签在线转富文本
文章浏览阅读6.7w次,点赞18次,收藏159次。GIS派文章地址:13个免费的开源GIS软件QGISQGIS是一个开放源码的地理信息系统。该项目诞生于2002年5月,并于同年6月作为SourceForge上的一个项目建立。我们一直在努力使GIS软件(传统上是昂贵的专有软件)成为任何人都可以使用个人电脑的可行前景。QGIS目前运行在大多数Unix平台、Windows和macOS上。QGIS是使用Qt工具包(https://www.qt.io)..._开源 gis
文章浏览阅读5.3k次,点赞4次,收藏8次。报错:Tomcat Manager App--403 Access Denied You are not authorized to view this page解决办法:1.关闭Tomcat(安装目录下bin文件下shutdown批处理文件)2.找到安装目录下conf文件夹tomcat-users.xml文件3.在xml文件<tomcat-users>&l..._you are not authorized to view this page.
文章浏览阅读8w次,点赞438次,收藏1.8k次。一.Javascript组成JavaScript的实现包括以下3个部分:ECMAScript(核心)描述了JS的语法和基本对象。文档对象模型 (DOM)处理网页内容的方法和接口浏览器对象模型(BOM)与浏览器交互的方法和接口javascript 有三部分构成,ECMAScript,DOM和BOM,根据宿主(浏览器)的不同,具体的表现形式也不尽相同,ie和其他的浏览器风格迥异,IE 扩展了 BOM,..._element.lastchild返回元素的最后一个子元素。element.namespaceuri返回元素的nam
文章浏览阅读1.1k次。第一种:native为id自动生成策略,生成的是数字id,添加数据到MySQL数据库时不需要设置id的值,hibernate框架会帮你生成,但是会给框架执行时间造成压力。第二种:assigned(手动委派)主键策略需要在添加数据时自己设置id,因为它不能自动生成id,不麻烦别人,自己动手。第三种:uuid.hex程序会自动生成16进制uuid主键,添加数据到MySQL数据库时不..._hibernate id生成器配置uuid和assign有什么区别
文章浏览阅读594次。第二周_决策项 评分分 神经网络