TopCoder SRM502 Div1 1000 动态规划
原文链接https://www.cnblogs.com/zhouzhendong/p/SRM502-1000.html
SRM502 Div1 1000
题意
从 [0,n-1] 中选择 k 个不同的整数,使得他们的和是 n 的倍数,求方案数。对 \(10^9+7\) 取模。
\(n\leq 10^9,k\leq 1000\)
题解
首先我们考虑从 n 个里面选择 k 个并进行排列的方案数,最终只需要除以 k! 就好了。
设 \(M=n\) ;
设 \(f(n,m,t)\) 表示 在 \(0,1,\cdots M-1\) 中任选 \(n\) 个,并使得 \((\sum_{1\leq i < n } x_i ) + t x_n \equiv 0 \pmod m\cdots (1)\) (设取的第 \(i\) 个为 \(x_i\) ) 。
则答案显然是 \(f(k,n,1)\) 。
于是我们考虑如何求解这个函数。
\(f(n,m,t)=\)
如果 $n=0 $ ,那么返回 \(1\) 。
如果 \(m = 1\) ,那么 (1) 式恒成立,答案就是 \(n!\binom{M}{n} = M^{\underline{n}}\)
否则,我们通过容斥,分类讨论(结果就是下面 (1) 的贡献减掉 (2) 的贡献)。
(1) \(x_i(i<n)\) 中可能存在与 \(x_n\) 相等的数(也可能不存在):相当于选择 \(n-1\) 个数存在 \(x_n\) 使得 \((\sum_{1\leq i < n } x_i ) + t x_n \equiv 0 \pmod m\) ,设 \(g = \gcd(t,m)\),则方案数 \(=f(n-1,g,1)\) ,对于 \(x_n\) ,我们也可以确定其取值个数,即 \(Mg/m\) 。
(2) \(x_i(i<n)\) 中存至少一个与 \(x_n\) 相等的数。则在 \(x_1\cdots x_{n-1}\) 个中任选一个和 \(x_n\) 相等,有 \(n-1\) 种可能,乘上对应的方案数 \(f(n-1,m,t+1)\) 即可。
由于我们经常要用到 \(f(a,b,1)\) ,所以我们对于 \(f(a,b,1)\) 记忆化一下,然后搜索即可。注意 \(f(a,n,1)\) 这种要特殊处理。
代码
static const int N=1005,mod=1e9+7;
int M;
int dp[N][N],Fac[N];
int Pow(int x,int y){
int ans=1;
for (;y;y>>=1,x=1LL*x*x%mod)
if (y&1)
ans=1LL*ans*x%mod;
return ans;
}
int gcd(int x,int y){
return y?gcd(y,x%y):x;
}
int f(int n,int m,int t){
if (n==0)
return 1;
if (m==1)
return Fac[n];
int g=gcd(m,t);
// sum + t * x = 0
// c(sum = 0) - c(sum = 0 && (t+1))
if (t==1){
int _m=min(m,1000+1);
if (~dp[n][_m])
return dp[n][_m];
dp[n][_m]=((1LL*f(n-1,g,1)*(M/m*g)
-1LL*f(n-1,m,t+1)*(n-1))%mod+mod)%mod;
return dp[n][_m];
}
int res=((1LL*f(n-1,g,1)*(M/m*g)
-1LL*f(n-1,m,t+1)*(n-1))%mod+mod)%mod;
return res;
}
int find(int N, int K){
M=N;
memset(dp,-1,sizeof dp);
int t=1;
Fac[0]=1;
for (int i=1;i<=K;i++){
t=1LL*t*i%mod;
Fac[i]=1LL*Fac[i-1]*(N-i+1)%mod;
}
t=Pow(t,mod-2);
int ans=1LL*t*f(K,N,1)%mod;
return ans;
}
TopCoder SRM502 Div1 1000 动态规划的更多相关文章
- TopCoder SRM500 Div1 1000 其他
原文链接https://www.cnblogs.com/zhouzhendong/p/SRM500-1000.html SRM500 Div1 1000 设 \(v_1,v_2,\cdots ,v_9 ...
- TopCoder SRM502 Div1 500 贪心 01背包
原文链接https://www.cnblogs.com/zhouzhendong/p/SRM502-500.html SRM502 Div1 500 好题. 首先,如果已经确定了解决所有问题的优先级, ...
- TopCoder 649 div1 & div2
最近一场TC,做得是在是烂,不过最后challenge阶段用一个随机数据cha了一个明显错误的代码,最后免于暴跌rating,还涨了一点.TC题目质量还是很高的,非常锻炼思维,拓展做题的视野,老老实实 ...
- TopCoder SRM500 Div1 250 其他
原文链接https://www.cnblogs.com/zhouzhendong/p/SRM500-250.html SRM500 Div1 250 题意 (看题用了半个小时--) 有 n 个人(编号 ...
- TopCoder SRM500 Div1 500 分治
原文链接https://www.cnblogs.com/zhouzhendong/p/SRM500-500.html SRM500 Div1 500 没想到 double 的精度居然没有爆-- 考虑以 ...
- topcoder SRM712 Div1 LR
题目: Problem Statement We have a cyclic array A of length n. For each valid i, element i-1 the l ...
- TopCoder 603 div1 & div2
div2 250pts MiddleCode 题意:s串长度为奇数时,将中间字符取掉并添加到t末尾:长度为偶数时,将中间两个较小的字符取掉并添加到末尾. 分析:直接做,学习了一下substr(s, p ...
- TopCoder SRM704 Div1 800 构造
原文链接https://www.cnblogs.com/zhouzhendong/p/SRM704-800.html 题解 考虑构造一个 $n = 20$ 的图. 先把所有 $i$ 都连向 $i-1$ ...
- topcoder SRM642 div1 hard WheelofFortune
题目链接:vjudge 大意:有两个人参加一场游戏,这个游戏在一个编号为\(0\text~n-1\)的轮盘上进行,一开始轮盘上的数字均为0:一共有\(m\)轮,每一轮都有一个操作参数\(s_i\),主 ...
随机推荐
- Ex 4_10 给定一个有向图G=(V,E),其中边...(bellman-ford算法的应用).._第十二次作业
在bellman-ford算法中,循环n-1(n为顶点个数)次可以找出从源点到其他顶点的最多n-1条边的最短路径,若循环k次则可以找出从源点到其他顶点的最多k条边的最短路径. package org. ...
- Spring加载加密的配置文件
一.继承并实现自己的属性文件配置器类 /** * 带加密的Spring属性配置文件扩展类 * 加密方式:AES * @author simon * */ public class EncryptPro ...
- 1)django-建立步骤和目录说明
一:前言 django是python最流行的WEB框架. 二:django安装 pip install django 三:django项目建立步骤 1.创建django工程 django-admin ...
- Net 4.5 WebSocket 在 Windows 7, Windows 8 and Server 2012上的比较
.Net 4.5 WebSocket Server Running on Windows 7? Net 4.5 WebSocket Server 可以运行在 Windows 7,但是Net 4.5的 ...
- 整理oracle 树形查询
注:本文参考了<整理oracle 树形查询> sql树形递归查询是数据库查询的一种特殊情形,也是组织结构.行政区划查询的一种最常用的的情形之一.下面对该种查询进行一些总结: create ...
- CentOS7图形界面与命令行界面切换(转载)
在图形界面使用 ctrl+alt+F2切换到dos界面 dos界面 ctrl+alt+F2切换回图形界面 在命令上 输入 init 3 命令 切换到dos界面 输入 init 5命令 切换到图形界面 ...
- ERROR 1045 (28000): Access denied for user 'mysql'@'localhost' (using password: YES
一.有时可以直接输入命令: mysql进入数据库 启动数据库:# mysqld_safe & 二.查看用户命令: mysql> use mysql; Reading table info ...
- Kali安装问题解决方案
一.对今天安装Kalilinux 失败的原因做一个简单的总结, 1.安装之前的电脑状况,电脑配置Windows7 64位旗舰版配置,虚拟版本是VMware Workstation Pro14 ver ...
- bzoj 4007
非常好的树形dp 首先,有个很显然的状态:记状态f[i][j]表示以i为根节点的子树中选了j个叶节点作战,那么很显然有转移:f[i][j1+j2]=f[i<<1][j1]+f[i<& ...
- Python字符串编码转换
使用encode()方法编码 str.encode([encoding="utf-8"][,errors="strict"]) str:表示需要转换的字符串 e ...