这玩意儿怎么看上去就很经典啊(

哦互不相同啊,那没事了(

考虑一个 \(\rm DP\)。由于限制了互不相同,那么我们考虑从值域开始想。

设 \(dp_{n,k}\) 为在 \([1,n]\) 中选了 \(k\) 个数的权值之和。

容易得到有:

\[dp_{n,k}=\sum_{i=1}^{n-1}dp_{i,k-1} \times n \times (N-k+1)
\]

很明显考虑 \(0\) 肯定也没有问题。

接下来开始生成函数:

\[F_k(x)=\sum_{i=0}^{\infty}dp_{i,k}x^i
\]

首先这个 \(N-k+1\) 肯定能放到最后去乘上一个阶乘,于是就变成了计算

\[dp_{n,k}=\sum_{i=0}^{n-1}dp_{i,k-1} \times n
\]

用 GF 来写的话,因为 \(n\) 的存在,考虑求导:

\[F_k(x)=x(\frac {xF_{k-1}(x)}{1-x})'
\]

然后我们有边界条件 \(F_0(x)=1\),答案是 \([x^n]\frac {F_k(x)}{1-x}\)。

注意到 \(F_k(x)\) 一定是类似 \(\frac {G(x)}{(1-x)^y}\) 之类的,所以考虑递推 \(F_k(x)\)。

根据 \(F_k(x)=x(\frac {xF_k(x)}{1-x})'\),我们直接保存 \(F_k(x)=\frac {H(x)}{(1-x)^y}\),然后暴力计算多项式乘法,除法求导利用 \((\frac {F(x)}{G(x)})'=\frac {F'(x)G(x)-F(x)G'(x)}{G^2(x)}\) 即可。

具体来说,计算这个东西:

\[(\frac {F(x)}{(1-x)^y})'
\]

我们知道 \((F^k(x))'=kF^{k-1}(x)F'(x)\):

\[\frac {F'(x)(1-x)^y+yF(x)(1-x)^{y-1}}{(1-x)^{2y}}
\]
\[\frac {F'(x)(1-x)+yF(x)}{(1-x)^{y+1}}
\]

答案是类似 \(\frac {F(x)}{(1-x)^{k+1}}\) 的形式,我们又知道 \(\frac 1 {(1-x)^{k+1}}=\sum_{i=0}^{\infty}\binom {i+k}ix^i\),直接卷就好了。

瓶颈在递推,复杂度 \(O(k^2)\)。

别忘记最后要乘上一个阶乘。

#include<cstdio>
typedef unsigned ui;
const ui M=505;
ui n,k,mod,f[M<<1],C[M<<1];
inline ui pow(ui a,ui b){
ui ans(1);for(a%=mod;b;b>>=1,a=1ull*a*a%mod)if(b&1)ans=1ull*ans*a%mod;return ans;
}
signed main(){
ui i,j,inx,ans(0);scanf("%u%u%u",&k,&n,&mod);C[inx=2*n+1>k?k:2*n+1]=f[0]=1;
for(i=1;i<=n;f[1]=f[0],f[0]=0,++i)for(j=i*2;j>=2;--j)f[j]=(1ull*f[j-1]*j+f[j-2]*(2ull*i-j))%mod;
if(inx==2*n+1)for(i=1;i<=2*n;++i)C[inx]=1ull*C[inx]*(k-i)%mod*pow(i,mod-2)%mod;
for(i=inx-1;i<=2*n;--i)C[i]=C[i+1]*(k+2ull*n-i)%mod*pow(k-i,mod-2)%mod;
for(i=0;i<=n*2;++i)ans=(ans+1ull*C[i]*f[i])%mod;for(i=1;i<=n;++i)ans=1ull*ans*i%mod;
printf("%u",ans);
}

LGP4463题解的更多相关文章

  1. 2016 华南师大ACM校赛 SCNUCPC 非官方题解

    我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...

  2. noip2016十连测题解

    以下代码为了阅读方便,省去以下头文件: #include <iostream> #include <stdio.h> #include <math.h> #incl ...

  3. BZOJ-2561-最小生成树 题解(最小割)

    2561: 最小生成树(题解) Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1628  Solved: 786 传送门:http://www.lyd ...

  4. Codeforces Round #353 (Div. 2) ABCDE 题解 python

    Problems     # Name     A Infinite Sequence standard input/output 1 s, 256 MB    x3509 B Restoring P ...

  5. 哈尔滨理工大学ACM全国邀请赛(网络同步赛)题解

    题目链接 提交连接:http://acm-software.hrbust.edu.cn/problemset.php?page=5 1470-1482 只做出来四道比较水的题目,还需要加强中等题的训练 ...

  6. 2016ACM青岛区域赛题解

    A.Relic Discovery_hdu5982 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Jav ...

  7. poj1399 hoj1037 Direct Visibility 题解 (宽搜)

    http://poj.org/problem?id=1399 http://acm.hit.edu.cn/hoj/problem/view?id=1037 题意: 在一个最多200*200的minec ...

  8. 网络流n题 题解

    学会了网络流,就经常闲的没事儿刷网络流--于是乎来一发题解. 1. COGS2093 花园的守护之神 题意:给定一个带权无向图,问至少删除多少条边才能使得s-t最短路的长度变长. 用Dijkstra或 ...

  9. CF100965C题解..

    求方程 \[ \begin{array}\\ \sum_{i=1}^n x_i & \equiv & a_1 \pmod{p} \\ \sum_{i=1}^n x_i^2 & ...

随机推荐

  1. 开发时Blocks跟Delegates如何选择----董鑫

    1.大多数delegate protocols 都拥有几个消息源. 以GKMatch为例(A GKMatch object provides a peer-to-peer network betwee ...

  2. Centos8安装virtualbox

    一.执行以下命令并启用 VirtualBox 和 EPEL 包仓库 dnf config-manager --add-repo=https://download.virtualbox.org/virt ...

  3. 安装Linux8.3.2011

    镜像地址:http://mirrors.aliyun.com/centos/8.3.2011/isos/x86_64/ 非DVD镜像安装时的安装源地址:http://mirrors.aliyun.co ...

  4. 面试题之java缓存总结,从单机缓存到分布式缓存架构

    1.缓存定义 高速数据存储层,提高程序性能 2.为什么要用缓存(读多写少,高并发) 1.提高读取吞吐量 2.提升应用程序性能 3.降低数据库成本 4.减少后端负载 5.消除数据库热点 6.可预测的性能 ...

  5. Vue 子组件更新父组件的值

    今天在使用Vue中遇到了一个新的需求:子组件需要修改由父组件传递过来的值,由于子组件的值是由父组件传递过来的,不能直接修改属性的值, 我们想改变传递过来的值只能通过自定义事件的形式修改父组件的值达到修 ...

  6. hashlib模块&日志模块

    内容概要 hashlib模块 logging模块 第三方模块下载 内容详细 hashlib模块 hashlib 是一个提供了一些流行的hash(摘要)算法的Python标准库.其中所包括的算法有 md ...

  7. Visual Studio Code 配置C、C++ 文件debug调试环境

    目录 vscode C/C++ Extension Pack 插件安装 vscode windows 端 debug 配置 window MinGW 环境安装 windows 端 C.CPP 单文件 ...

  8. 基于nodejs的阿里云DDNS服务,支持多网卡绑定

    前段几天写的,自用. 想着也许有和我一样需求的同学,所以开源发出来了. 提供阿里云的DDNS设置,支持多网卡多域名绑定. 特性 支持阿里云DDNS设置 支持针对多个线路多个网卡单独设置域名 支持IPv ...

  9. 从零开始,开发一个 Web Office 套件(5):Mouse hover over text

    <从零开始, 开发一个 Web Office 套件>系列博客目录 这是一个系列博客, 最终目的是要做一个基于HTML Canvas 的, 类似于微软 Office 的 Web Office ...

  10. 【转】k8s集群自定义clusterRole样例

    对pod资源可以删除,进入终端执行命令,其他资源只读权限 apiVersion: rbac.authorization.k8s.io/v1 kind: ClusterRole metadata: an ...