令 f ( p , k ) f(p,k) f(p,k)表示从 [ 0 , n − 1 ] [0,n-1] [0,n−1]中选择 k k k个数，使得和为 p p p的倍数的方案数。
接下来用一种很诡异的姿势容斥。
直接算并不好算，考虑枚举最后有 i i i个数固定且相同，前面选择 k − i k-i k−i个不同的数和为 s u m sum sum，这两部分分开计算，互不影响。
发现会出现 i + 1 i+1 i+1个数重复的情况，需要减掉，然后又多减掉了 i + 2 i+2 i+2个数重复的情况，就这样容斥下去。
考虑算固定 i i i个数的情况。
显然我们需要下面方程有解：
p ∣ i x + s u m p\mid ix+sum p∣ix+sum
则 g c d ( i , p ) ∣ s u m gcd(i,p)\mid sum gcd(i,p)∣sum，设 s u m = t ⋅ g c d ( i , p ) sum=t\cdot gcd(i,p) sum=t⋅gcd(i,p)，则：
p ∣ i x + t ⋅ g c d ( i , p ) p g c d ( i , p ) ∣ i g c d ( i , p ) x + t p\mid ix+t\cdot gcd(i,p)\\\frac{p}{gcd(i,p)}\mid \frac{i}{gcd(i,p)}x+t p∣ix+t⋅gcd(i,p)gcd(i,p)p∣gcd(i,p)ix+t
显然由于 p ∣ n p\mid n p∣n， x x x在 [ 0 , n − 1 ] [0,n-1] [0,n−1]中的合法取值有 n p g c d ( i , p ) = n p ⋅ g c d ( i , p ) \frac{n}{\frac{p}{gcd(i,p)}}=\frac{n}p\cdot gcd(i,p) gcd(i,p)pn=pn⋅gcd(i,p)，且与 t t t无关，也就是说子问题为 f ( g c d ( i , p ) , k − i ) f(gcd(i,p),k-i) f(gcd(i,p),k−i)。
容斥完了之后，需要考虑一个东西，在算没有任何数重复的合法方案的时候，我们实际上还是固定了一个数，则每一个方案被算了 k k k次，需要除掉。
那么，很容易地，我们得到容斥递归求解的式子：
f ( p , k ) = 1 k ∑ i = 1 k ( − 1 ) i + 1 ⋅ n ⋅ g c d ( i , p ) p ⋅ f ( g c d ( i , p ) , k − i ) f(p,k)=\frac{1}{k}\sum_{i=1}^{k}(-1)^{i+1}\cdot \frac{n\cdot gcd(i,p)}{p}\cdot f(gcd(i,p),k-i) f(p,k)=k1i=1∑k(−1)i+1⋅pn⋅gcd(i,p)⋅f(gcd(i,p),k−i)
总复杂度是一个调和级数 O ( k log k ) O(k\log k) O(klogk)？常数较大。
代码：
#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const
using std::cerr;
using std::cout;
inline int gcd(int a,int b){
if(!a||!b)return a|b;
int shift=__builtin_ctz(a|b);
for(b>>=__builtin_ctz(b);a;a-=b)if((a>>=__builtin_ctz(a))<b)std::swap(a,b);
return b<<shift;
}
cs int mod=1e9+7;
inline int add(int a,int b){
return (a+=b)>=mod?a-mod:a;}
inline int dec(int a,int b){
return (a-=b)<0?a+mod:a;}
inline int mul(int a,int b){
ll r=(ll)a*b;return r>=mod?r%mod:r;}
inline int power(int a,int b,int res=1){
for(;b;b>>=1,a=mul(a,a))(b&1)&&(res=mul(res,a));
return res;
}
inline void Inc(int &a,int b){
(a+=b)>=mod&&(a-=mod);}
inline void Dec(int &a,int b){
(a-=b)<0&&(a+=mod);}
cs int K=5e3+3;
int n,k;
int Cn[K],inv[K];
inline void init(){
inv[0]=inv[1]=1;
Cn[0]=1,Cn[1]=n;
for(int re i=2;i<=k;++i){
inv[i]=mul(inv[mod%i],mod-mod/i);
Cn[i]=mul(Cn[i-1],mul(dec(n+1,i),inv[i]));
}
}
int dp[K][K];
inline int f(int p,int k){
if(!k)return 1;
if(p==1)return Cn[k];
if(~dp[p][k])return dp[p][k];
int ans=0,tmp,val,d=n/p;
for(int re i=1;i<=k;++i){
tmp=gcd(p,i);
val=mul(f(tmp,k-i),tmp);
(i&1)?Inc(ans,val):Dec(ans,val);
}
return dp[p][k]=mul(ans,mul(d,inv[k]));
}
signed main(){
#ifdef zxyoi
freopen("slay.in","r",stdin);
#endif
std::cin>>n>>k;
init();
memset(dp,-1,sizeof dp);
int ans=0,tmp,val;
for(int re i=1;i<=k;++i){
tmp=gcd(n,i);
val=mul(f(tmp,k-i),tmp);
(i&1)?Inc(ans,val):Dec(ans,val);
}
cout<<mul(ans,inv[k])<<"\n";
return 0;
}
文章浏览阅读1k次。通过使用ajax方法跨域请求是浏览器所不允许的，浏览器出于安全考虑是禁止的。警告信息如下：不过jQuery对跨域问题也有解决方案，使用jsonp的方式解决，方法如下：$.ajax({ async:false, url: 'http://www.mysite.com/demo.do', // 跨域URL ty..._nginx不停的xhr
文章浏览阅读2k次。关于在 Oracle 中配置 extproc 以访问 ST_Geometry，也就是我们所说的 使用空间SQL 的方法，官方文档链接如下。http://desktop.arcgis.com/zh-cn/arcmap/latest/manage-data/gdbs-in-oracle/configure-oracle-extproc.htm其实简单总结一下，主要就分为以下几个步骤。..._extproc
文章浏览阅读1.5w次。linux下没有上面的两个函数，需要使用函数 mbstowcs和wcstombsmbstowcs将多字节编码转换为宽字节编码wcstombs将宽字节编码转换为多字节编码这两个函数，转换过程中受到系统编码类型的影响，需要通过设置来设定转换前和转换后的编码类型。通过函数setlocale进行系统编码的设置。linux下输入命名locale -a查看系统支持的编码_linux c++ gbk->utf8
文章浏览阅读750次。今天准备从生产库向测试库进行数据导入，结果在imp导入的时候遇到“ IMP-00009:导出文件异常结束” 错误，google一下，发现可能有如下原因导致imp的数据太大，没有写buffer和commit两个数据库字符集不同从低版本exp的dmp文件，向高版本imp导出的dmp文件出错传输dmp文件时，文件损坏解决办法：imp时指定..._imp-00009导出文件异常结束
文章浏览阅读143次。当下是一个大数据的时代，各个行业都离不开数据的支持。因此，网络爬虫就应运而生。网络爬虫当下最为火热的是Python，Python开发爬虫相对简单，而且功能库相当完善，力压众多开发语言。本次教程我们爬取前程无忧的招聘信息来分析Python程序员需要掌握那些编程技术。首先在谷歌浏览器打开前程无忧的首页，按F12打开浏览器的开发者工具。浏览器开发者工具是用于捕捉网站的请求信息，通过分析请求信息可以了解请..._初级python程序员能力要求
文章浏览阅读7.6k次，点赞2次，收藏6次。@Service标注的bean，类名：ABDemoService查看源码后发现，原来是经过一个特殊处理：当类的名字是以两个或以上的大写字母开头的话，bean的名字会与类名保持一致public class AnnotationBeanNameGenerator implements BeanNameGenerator { private static final String C..._@service beanname
文章浏览阅读6.9w次，点赞73次，收藏463次。1.前序创建#include<stdio.h>#include<string.h>#include<stdlib.h>#include<malloc.h>#include<iostream>#include<stack>#include<queue>using namespace std;typed_二叉树的建立
文章浏览阅读7.1k次。在Asp.net上使用Excel导出功能，如果文件名出现中文，便会以乱码视之。 解决方法： fileName = HttpUtility.UrlEncode(fileName, System.Text.Encoding.UTF8);_asp.net utf8 导出中文字符乱码
文章浏览阅读2.1k次，点赞4次，收藏23次。第一次实验 词法分析实验报告设计思想词法分析的主要任务是根据文法的词汇表以及对应约定的编码进行一定的识别，找出文件中所有的合法的单词，并给出一定的信息作为最后的结果，用于后续语法分析程序的使用；本实验针对 PL/0 语言 的文法、词汇表编写一个词法分析程序，对于每个单词根据词汇表输出： (单词种类, 单词的值) 二元对。词汇表：种别编码单词符号助记符0beginb..._对pl/0作以下修改扩充。增加单词
文章浏览阅读773次。我在使用adb.exe时遇到了麻烦.我想使用与bash相同的adb.exe shell提示符,所以我决定更改默认的bash二进制文件(当然二进制文件是交叉编译的,一切都很完美)更改bash二进制文件遵循以下顺序> adb remount> adb push bash / system / bin /> adb shell> cd / system / bin> chm..._adb shell mv 权限
文章浏览阅读6.8k次，点赞12次，收藏125次。1. 单目相机标定引言相机标定已经研究多年，标定的算法可以分为基于摄影测量的标定和自标定。其中，应用最为广泛的还是张正友标定法。这是一种简单灵活、高鲁棒性、低成本的相机标定算法。仅需要一台相机和一块平面标定板构建相机标定系统，在标定过程中，相机拍摄多个角度下（至少两个角度，推荐10~20个角度）的标定板图像（相机和标定板都可以移动），即可对相机的内外参数进行标定。下面介绍张氏标定法（以下也这么称呼）的原理。原理相机模型和单应矩阵相机标定，就是对相机的内外参数进行计算的过程，从而得到物体到图像的投影_相机-投影仪标定
文章浏览阅读2.2k次。文章目录Wayland 架构Wayland 渲染Wayland的 硬件支持简 述： 翻译一篇关于和 wayland 有关的技术文章, 其英文标题为Wayland Architecture .Wayland 架构若是想要更好的理解 Wayland 架构及其与 X (X11 or X Window System) 结构；一种很好的方法是将事件从输入设备就开始跟踪, 查看期间所有的屏幕上出现的变化。这就是我们现在对 X 的理解。 内核是从一个输入设备中获取一个事件，并通过 evdev 输入_wayland