P4345 [SHOI2015]超能粒子炮·改

P4345 [SHOI2015]超能粒子炮·改

题意:

求解式子∑i=0kCni%p\sum_{i=0}^{k}C_{n}^{i} \% pi=0kCni%p
n,k<=1e18

题解:

f(n,k)=∑i=0kCnif(n,k)=\sum_{i=0}^{k}C_{n}^{i}f(n,k)=i=0kCni
开始化简:
由卢卡斯定理得:
f(n,k)=∑i=0kCni=∑i=0kCn/pi/p∗Cn%pi%pf(n,k)=\sum_{i=0}^{k}C_{n}^{i}=\sum_{i=0}^{k}C_{n/p}^{i/p}*C_{n\%p}^{i\%p}f(n,k)=i=0kCni=i=0kCn/pi/pCn%pi%p

我们将相同的Cn/pi/pC_{n/p}^{i/p}Cn/pi/p为一组将其展开:
Cn/p0∑i=0p−1Cn%pi=Cn/p1∑i=0p−1Cn%pi+...+Cn/pk/p∑i=0k%pCn%piC_{n/p}^{0}\sum_{i=0}^{p-1}C_{n\%p}^{i}=C_{n/p}^{1}\sum_{i=0}^{p-1}C_{n\%p}^{i}+...+C_{n/p}^{k/p}\sum_{i=0}^{k\%p}C_{n\%p}^{i}Cn/p0i=0p1Cn%pi=Cn/p1i=0p1Cn%pi+...+Cn/pk/pi=0k%pCn%pi

此时有0到k/p-1这些整块,还有k/p这个不完整的块
我们先考虑完整的块,
我们将∑i=0p−1Cn%pi\sum_{i=0}^{p-1}C_{n\%p}^{i}i=0p1Cn%pi提出来
变成:∑i=0p−1Cn%pi∗(Cn/p0+...+Cn/pk/p−1)\sum_{i=0}^{p-1}C_{n\%p}^{i}*(C_{n/p}^{0}+...+C_{n/p}^{k/p-1})i=0p1Cn%pi(Cn/p0+...+Cn/pk/p1)
我们一开始定义:
f(n,k)=∑i=0kCnif(n,k)=\sum_{i=0}^{k}C_{n}^{i}f(n,k)=i=0kCni
这就可以写成:
f(n%p,p−1)∗f(n/p,k/p−1)f(n\%p,p-1)*f(n/p,k/p-1)f(n%p,p1)f(n/p,k/p1)

再考虑不完整的第k/p块
∑i=0k%pCn%pi\sum_{i=0}^{k\%p}C_{n\%p}^{i}i=0k%pCn%pi也可以写成f(n%p,k%p)f(n\%p,k\%p)f(n%p,k%p)

总结:
f(n,k)=f(n%p,p−1)∗f(n/p,k/p−1)+Cn/pk/p∗f(n%p,k%p)f(n,k)=f(n\%p,p-1)*f(n/p,k/p-1)+C_{n/p}^{k/p}*f(n\%p,k\%p)f(n,k)=f(n%p,p1)f(n/p,k/p1)+Cn/pk/pf(n%p,k%p)
Cn/pk/pC_{n/p}^{k/p}Cn/pk/p可以用Lucas求
f(n%p,p−1)f(n\%p,p-1)f(n%p,p1)f(n%p,k%p)f(n\%p,k\%p)f(n%p,k%p)因为数值不大,直接预初理就行(杨辉三角再求和)
复杂度O(p2+Tlog23332n)O(p^{2}+Tlog^{2}_{2333}n)O(p2+Tlog23332n)

代码:

#include <bits/stdc++.h>
#include <unordered_map>
#define debug(a, b) printf("%s = %d\n", a, b);
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> PII;
clock_t startTime, endTime;
//Fe~Jozky
const ll INF_ll= 1e18;
const int INF_int= 0x3f3f3f3f;
void read(){};
template <typename _Tp, typename... _Tps> void read(_Tp& x, _Tps&... Ar)
{x= 0;char c= getchar();bool flag= 0;while (c < '0' || c > '9')flag|= (c == '-'), c= getchar();while (c >= '0' && c <= '9')x= (x << 3) + (x << 1) + (c ^ 48), c= getchar();if (flag)x= -x;read(Ar...);
}
template <typename T> inline void write(T x)
{if (x < 0) {x= ~(x - 1);putchar('-');}if (x > 9)write(x / 10);putchar(x % 10 + '0');
}
void rd_test()
{
#ifdef LOCALstartTime= clock();freopen("in.txt", "r", stdin);
#endif
}
void Time_test()
{
#ifdef LOCALendTime= clock();printf("\nRun Time:%lfs\n", (double)(endTime - startTime) / CLOCKS_PER_SEC);
#endif
}
const int maxn= 3000;
const int P= 2333;
ll c[maxn + 2][maxn + 2];
ll f[maxn + 2][maxn + 2];
inline ll Lucas(ll a, ll b)
{if (!b)return 1;if (a == b)return 1;if (a < b)return 0;return c[a % P][b % P] * Lucas(a / P, b / P) % P;
}
inline ll F(ll a, ll k)
{if (k < 0)return 0;if (!a)return 1;if (!k)return 1;if (a < P && k < P)return f[a][k];return (f[a % P][P - 1] * F(a / P, k / P - 1) % P + Lucas(a / P, k / P) * f[a % P][k % P] % P) % P;
}
int main()
{int T;scanf("%d", &T);c[0][0]= 1;for (int i= 1; i <= maxn; i++)c[i][i]= c[i][0]= 1;for (int i= 1; i <= maxn; i++)for (int j= 1; j < i; j++)c[i][j]= (c[i - 1][j] + c[i - 1][j - 1]) % P;f[0][0]= 1;for (int i= 1; i <= maxn; i++)f[i][0]= 1;for (int i= 0; i <= maxn; i++)for (int j= 1; j <= maxn; j++)f[i][j]= (c[i][j] + f[i][j - 1]) % P;ll a, k;while (T--) {scanf("%lld%lld", &a, &k);printf("%lld\n", F(a, k));}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/316218.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

.NET Framework VS .NET Core

本文对应的原文来至 c-sharpcorner 的一篇文章&#xff0c;文末有链接。如有错误&#xff0c;还请指正。前言你会为你的下一个应用程序选择哪一种开发平台 - .NET Framework 或者 .NET Core&#xff1f;在这篇文章中&#xff0c;让我们比较一下这两个开发平台的特点&#xff0c…

P6669 [清华集训2016] 组合数问题

P6669 [清华集训2016] 组合数问题 题意&#xff1a; 给你n&#xff0c;m&#xff0c;k&#xff0c;问有多少对(i,j)满足K∣CijK|C_{i}^{j}K∣Cij​ (Cij是k的倍数C_{i}^{j}是k的倍数Cij​是k的倍数) n,m<1e18 题解&#xff1a; n和m非常大&#xff0c;非常非常大&#x…

解决 VS2019 中.net core WPF 暂时无法使用 Designer 的临时方法

以下方法来自于微软github开源项目WPF:dotnet/samples - WPF Hello World sample with linked files&#xff0c;请放心使用。此篇文章是上篇文章解决vs2019中暂时无法为.net core WinForms使用 Designer 的临时方法的姊妹篇&#xff0c;但对WPF而言实现起来比WinForms简单很多…

Building Fire Stations

题目链接 首先找到树的直径&#xff0c;直径左端点是a&#xff0c;直径右端点是b&#xff0c;中间的点是mid&#xff08;偶数的情况下mid可以看做两个&#xff09;&#xff0c;两点因该是左右分布&#xff1b; 假设两点都不在直径上&#xff0c;那么移到直径上的话距离更短&a…

.NET Core 迁移躺坑记

最近将自己负责的一个核心接口系统从.Net Framework迁移到了.Net Core。整体过程&#xff0c;从业务层面说一般般吧(整体还好但还是搞的业务有感&#xff0c;没出严重故障&#xff09;但是技术层面上感觉其实并没有达到要求&#xff0c;不过预期也是应该不会那么顺利&#xff0…

P4720 【模板】扩展卢卡斯定理/exLucas(无讲解,纯记录模板)

P4720 【模板】扩展卢卡斯定理/exLucas 题意&#xff1a; CnmmodpC_{n}^{m}\bmod pCnm​modp 对于 100% 的数据&#xff0c;1≤m≤n≤1018&#xff0c;2≤p≤106&#xff0c;不保证 p 是质数。 题解&#xff1a; 模板题&#xff0c;单纯写本文章记录板子 代码&#xff1a; …

[Abp vNext 源码分析] - 2. 模块系统的变化

一、简要说明本篇文章主要分析 Abp vNext 当中的模块系统&#xff0c;从类型构造层面上来看&#xff0c;Abp vNext 当中不再只是单纯的通过 AbpModuleManager 来管理其他的模块&#xff0c;它现在则是 IModuleManager 和 IModuleLoader 来协同工作&#xff0c;其他的代码逻辑并…

P3301 [SDOI2013]方程

P3301 [SDOI2013]方程 题意&#xff1a; 题解&#xff1a; 插板法介绍 首先要先讲组合数学的一个方法&#xff1a;插板法 问题引出&#xff1a;把10个球放进三个盒子&#xff0c;每个箱子至少一个有多少种分法&#xff1f; 10个球就有9个空隙&#xff0c;我们可以考虑在这个…

201912-3 化学方程式

他这个好像之和大写字母有关系&#xff1b; 小写字母跟着前面的的大写字母&#xff1b; 和代表要处理了&#xff1b; &#xff08;&#xff09;代表要乘了&#xff1b; #include<iostream> #include<cstdio> #include<algorithm> #include<cstring>…

.NET Framework 4.8发布

原文地址&#xff1a;https://devblogs.microsoft.com/dotnet/announcing-the-net-framework-4-8/我们很高兴地宣布今天发布.NET Framework 4.8。它包含在Windows 10 May 2019更新中。.NET Framework 4.8也可在Windows 7和Windows Server 2008 R2 上使用。您可以从我们的 .NET下…

Deltix Round, Summer 2021 (open for everyone, rated, Div. 1 + Div. 2)

Deltix Round, Summer 2021 (open for everyone, rated, Div. 1 Div. 2) 题号题目知识点AA Variety of OperationsBTake Your Places!CCompressed Bracket SequenceDTake a GuessEEquilibriumFSports BettingGGates to Another WorldHDIY Tree

[NewLife.XCode]数据层缓存(网站性能翻10倍)

NewLife.XCode是一个有10多年历史的开源数据中间件&#xff0c;支持nfx/netcore&#xff0c;由新生命团队(2002~2019)开发完成并维护至今&#xff0c;以下简称XCode。整个系列教程会大量结合示例代码和运行日志来进行深入分析&#xff0c;蕴含多年开发经验于其中&#xff0c;代…

cf1556A. A Variety of Operations

cf1556A. A Variety of Operations 题意&#xff1a; 有两个数a&#xff0c;b一开始都是0&#xff0c;现在有三种操作&#xff1a; 给a和b都加ka加k&#xff0c;b减ka减k&#xff0c;b加k 问从a0&#xff0c;b0到ac&#xff0c;bd最少需要几步&#xff1f; 题解&#xff1a…

[NewLife.XCode]高级查询(化繁为简、分页提升性能)

NewLife.XCode是一个有10多年历史的开源数据中间件&#xff0c;支持nfx/netcore&#xff0c;由新生命团队(2002~2019)开发完成并维护至今&#xff0c;以下简称XCode。整个系列教程会大量结合示例代码和运行日志来进行深入分析&#xff0c;蕴含多年开发经验于其中&#xff0c;代…

cf1556B B. Take Your Places!

cf1556B B. Take Your Places! 题意&#xff1a; 有n个数&#xff0c;你可以将相邻两个数交换&#xff0c;使得奇偶性一样的数不相邻。问最少操作步数 题解&#xff1a; 最终排列无非是&#xff1a;奇&#xff0c;偶&#xff0c;奇…或者偶&#xff0c;奇&#xff0c;偶… …

微软云Azure训练营 | 八城联动,全球盛会

Global Azure Bootcamp是由微软发起、MVP参与组织的全球化学习交流活动。每年会挑选一个特定的时间&#xff0c;在同一天内&#xff0c;全球不同地区将同时开展。2019年全球Azure训练营&#xff08;Global Azure Bootcamp&#xff09;将于2019年4月27日在全球270多个城市同时举…

2021-09-211547G - How Many Paths?

#include <iostream> #include <cstdio> #include <cstring> #include <algorithm> #include <set> using namespace std; typedef long long ll;const int N4e510,M1e610,mod998244353;int h[N],hs[N],e[M],ne[M],idx0; // h 原图 hs新图 v…

cf1556Compressed Bracket Sequence

cf1556Compressed Bracket Sequencex 题意&#xff1a; 给你n个数&#xff0c;奇数位置上的数表示左括号的数量&#xff0c;偶数位置上的数表示右括号的数量。问有多少个[l,r]是满足括号匹配的 题解&#xff1a; 括号匹配也算是经典问题了 直接统计不好计算&#xff0c;我们…

长沙4月21日开发者大会暨.NET社区成立大会活动纪实

活动总结2019年4月21日是一个斜风细雨、微风和煦的美好日子&#xff0c;由长沙.NET技术社区、腾讯云云加社区、微软Azure云技术社区、中国.NET技术社区、长沙柳枝行动、长沙互联网活动基地&#xff08;唐胡子俱乐部&#xff09;等多家单位共同主办的长沙开发者技术大会暨长沙.N…

树上启发式合并 简单例题

U41492 树上数颜色 #include <iostream> #include <cstdio> #include <cstring> #include <algorithm> #include <set> #include <vector> #include <queue> using namespace std; typedef long long LL; typedef pair<int,int…