[bzoj2173]整数的lqp拆分

Description

lqp在为出题而烦恼,他完全没有头绪,好烦啊… 他首先想到了整数拆分。整数拆分是个很有趣的问题。给你一个正整数N,对于N的一个整数拆分就是满足任意m>0,a1 ,a2 ,a3…am>0,且a1+a2+a3+…+am=N的一个有序集合。通过长时间的研究我们发现了计算对于N的整数拆分的总数有一个很简单的递推式,但是因为这个递推式实在太简单了,如果出这样的题目,大家会对比赛毫无兴趣的。然后lqp又想到了斐波那契数。定义F0=0,F1=1,Fn=Fn-1+Fn-2 (n>1),Fn就是斐波那契数的第n项。但是求出第n项斐波那契数似乎也不怎么困难… lqp为了增加选手们比赛的欲望,于是绞尽脑汁,想出了一个有趣的整数拆分,我们暂且叫它:整数的lqp拆分。和一般的整数拆分一样,整数的lqp拆分是满足任意m>0,a1 ,a2 ,a3…am>0,且a1+a2+a3+…+am=N的一个有序集合。但是整数的lqp拆分要求的不是拆分总数,相对更加困难一些。对于每个拆分,lqp定义这个拆分的权值Fa1Fa2…Fam,他想知道对于所有的拆分,他们的权值之和是多少?简单来说,就是求 由于这个数会十分大,lqp稍稍简化了一下题目,只要输出对于N的整数lqp拆分的权值和mod 109(10的9次方)+7输出即可。

Input

输入的第一行包含一个整数N。

Output

输出一个整数,为对于N的整数lqp拆分的权值和mod 109(10的9次方)+7。

Sample Input

3

Sample Output

5

【样例说明】

F0=0,F1=1,F2=1,F3=2。

对于N=3,有这样几种lqp拆分:

3=1+1+1, 权值是1*1*1=1。

3=1+2,权值是1*2=2。

3=2+1,权值是2*1=2。

所以答案是1*1*1+1*2+2*1=5。

HINT

20%数据满足:1≤N≤25 50%数据满足:1≤N≤1000 100%数据满足:1≤N≤1000000

给出两种解法

1.

母函数还不知道,先当逆元和组合数练习题

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #include<cstring>
 5 using namespace std;
 6 typedef long long lol;
 7 int Mod=1000000007;
 8 lol ans,A[1000001],c[1000001];
 9 int n;
10 lol pow(lol x,int y)
11 {
12     lol res=1;
13     while (y)
14     {
15         if (y%2==1) res=(res*x)%Mod;
16         x=(x*x)%Mod;
17         y=y/2;
18     }
19     return res;
20 }
21 int main()
22 {int i;
23     cin>>n;
24     A[1]=1;
25      for (i=2;i<=n;i++)
26      A[i]=((Mod-Mod/i)*A[Mod%i])%Mod;
27      c[1]=n;
28      ans=(ans+n)%Mod;
29      for (i=1;i<=(n-1)/2;i++)
30      {
31          c[2*i+1]=(((c[2*i-1]*(n-2*i))%Mod)*(n-2*i+1)%Mod)%Mod;
32          c[2*i+1]=(((c[2*i+1]*A[2*i])%Mod)*A[2*i+1])%Mod;
33          ans=(ans+(c[2*i+1]*pow(2,i))%Mod)%Mod;
34      }
35     cout<<ans;
36 }

2.打表找规律:

G[n]=2G[n-1]+G[n-2],G[0]=0,G[1]=1

不过可以通过40分dp推导

暴力显然有G[n]=ΣG[i]*F[n-i]+F[n],O(n^2)

1 #include<cstdio>
2 int n;long long mod=1e9+7,g[1000010];
3 int main()
4 {
5      scanf("%d",&n);g[1]=1;
6      for(int i=2;i<=n;i++)g[i]=(2*g[i-1]+g[i-2])%mod;
7      printf("%lld\n",g[n]);
8 }
时间: 2024-10-14 17:37:18

[bzoj2173]整数的lqp拆分的相关文章

整数的lqp拆分

题目大意 lqp在为出题而烦恼,他完全没有头绪,好烦啊- 他首先想到了整数拆分.整数拆分是个很有趣的问题.给你一个正整数N,对于N的一个整数拆分就是满足任意m>0,a1 ,a2 ,a3-am>0,且a1+a2+a3+-+am=N的一个有序集合.通过长时间的研究我们发现了计算对于N的整数拆分的总数有一个很简单的递推式,但是因为这个递推式实在太简单了,如果出这样的题目,大家会对比赛毫无兴趣的. 然后lqp又想到了斐波那契数.定义F0=0,F1=1,Fn=Fn-1+Fn-2 (n>1),Fn就

[国家集训队] 整数的lqp拆分

题面:https://www.luogu.org/problemnew/show/P4451 题解链接:https://www.luogu.org/problemnew/solution/P4451 题解见代码注释 (或者我回来会放我在洛谷上写的题解) /* 设 g[i] 为i的lqp拆分的权值和,则 g[i] = ∑f[j] * g[i-j] + f[i],其中 g[0] = 0, g[1] = 1. 以前是推式子推出来结果的,那么今天就尝试用生成函数做一下: 设 A = ∑f[i] * x^

python写的大整数相乘的方法

输入 72106547548473106236 982161082972751393 两个大整数 输出结果 70820244829634538040848656466105986748 解题思路 首先根据 大整数相乘的原理的基础上,把大整数进行优化拆分,拆分的长度,要考虑语言中整形的长度.这里用的python,其实可以直接乘的,呵呵.我暂定设为了 4,也就是 说 7210 6547 5484 7310 6236 98 2161 0829 7275 1393 然后 逐份相乘,当然千万别忘了结果要补

UVa 10313 - Pay the Price

题目:求一个整数的重复拆分,限制拆分数的个数. 分析:dp,二维多重背包.整数拆分就用背包. 状态:设f(i,j)为j拆分成i个元素的拆法: 转移:f(i,j)= sum(f(i-1,j-k),f(i-1,j-2k),...,f(i-1,j-mk)){ 其中,1 ≤ k ≤ j }: 因为输入格式WA好多次,强大的sscanf( ⊙ o ⊙ )啊! 说明:注意数据范围,使用long long防止溢出. #include <iostream> #include <cstdlib> #

Lambda表达式的五个练习

Lambda表达式的五练习 语法: Lambda函数,是一个匿名函数,创建语法: lambda parameters:express parameters:可选,如果提供,通常是逗号分隔的变量表达式形式,即位置参数. expression:不能包含分支或循环(但允许条件表达式),也不能包含return(或yield)函数.如果为元组,则应用圆括号将其包含起来. 调用lambda函数,返回的结果是对表达式计算产生的结果. 1.题目:利用条件运算符的嵌套来完成此题:学习成绩>=90分的同学用A表示,

一些数论知识

1.本原勾股数: 概念:一个三元组(a,b,c),其中a,b,c没有公因数而且满足:a^2+b^2=c^2 首先,这种本原勾股数的个数是无限的,而且构造的条件满足: a=s*t,b=(s^2-t^2)/2,c=(s^2+t^2)/2 其中s>t>=1是任意没有公因数的奇数! 由以上概念就可以导出任意一个本原勾股数组. 2.素数计数(素数定理) 令π(x)为1到x中素数的个数 19世纪最高的数论成就就是以下这个玩意儿: lim(x->∞){π(x)/(x/ln(x))}=1 数论最高成就,

1028:Ignatius and the Princess III

本题应该有两种方法: 1.母函数法 2.递推法 母函数不了解,待充分了解之后,再进行补充! 这里为递推实现的方法: 思路: 定义:n为要拆分的整数: k为拆分的项数: f[n][k]代表 n的整数拆分中,最大项不超过k的方案数. 每一个整数n的拆分中,总有一项拆分为自己,即:n = n; 我们将其表示为f[n][1],而且f[n][1] = 1; 又,每一个整数n的拆分中,总有一项拆分为n个1,即:n = 1 + 1 + ...... + 1(n个1的加和); 我们将其表示为f[0][0],且f

mysql那些事(3)小数如何存储

创建mysql数据表的时候,经常会遇到存储小数(浮点数)的情况,如:价格,重量,身高等. 目前大的公司流行三种存储方案: 1.将数据扩大10的倍数达到使用整数类型存储目的. 比如价格,我们经常以分为单位进行存储,也就是将数据扩大100倍,这样元改成分存储. 重量可以用克为单位,如果克还是小数,就以毫克,微克等单位进行存储. 2.使用decimal类型的方式存储 比如价格,我们可以使用两位小数的精度进行存储.强烈不建议采用float或者double的类型存储,会存在精度损失的问题,在以后做值比较的

LCA 算法(二)倍增

 介绍一种解决最近公共祖先的在线算法,倍增,它是建立在任意整数的二进制拆分之上.   代码:   1 //LCA:Doubly 2 3 #include<cstdio> 4 #define swap(a,b) a^=b^=a^=b 5 #define maxn 500010 6 using namespace std; 7 8 int n,m,s,tot,head[maxn],deep[maxn],p[maxn][20]; 9 struct node 10 { 11 int nxt,to; 1