转眼到了收获的季节,由于有TT的专业指导,Lele获得了大丰收。特别是水果,Lele一共种了N种水果,有苹果,梨子,香蕉,西瓜……不但味道好吃,样子更是好看。
于是,很多人们慕名而来,找Lele买水果。
甚至连大名鼎鼎的HDU ACM总教头 lcy 也来了。lcy抛出一打百元大钞,"我要买由M个水果组成的水果拼盘,不过我有个小小的要求,对于每种水果,个数上我有限制,既不能少于某个特定值,也不能大于某个特定值。而且我不要两份一样的拼盘。你随意搭配,你能组出多少种不同的方案,我就买多少份!"
现在就请你帮帮Lele,帮他算一算到底能够卖出多少份水果拼盘给lcy了。
注意,水果是以个为基本单位,不能够再分。对于两种方案,如果各种水果的数目都相同,则认为这两种方案是相同的。
最终Lele拿了这笔钱,又可以继续他的学业了~
Input
本题目包含多组测试,请处理到文件结束(EOF)。
每组测试第一行包括两个正整数N和M(含义见题目描述,0
Output
对于每组测试,在一行里输出总共能够卖的方案数。
题目数据保证这个答案小于10^9
Sample Input
2 3
1 2
1 2
3 5
0 3
0 3
0 3
Sample Output
2
12
这是我写的第一道母函数的题看了很多别人的blog才搞明白什么回事;
这里简单讲下要点,就是下面的这个公式理解好:
(x^(v[0]*n1[0])+x^(v[0]*(n1[0]+1))+x^(v[0]*(n1[0]+2))+...+x^(v[0]*n2[0]))
(x^(v[1]*n1[1])+x^(v[1]*(n1[1]+1))+x^(v[1]*(n1[1]+2))+...+x^(v[1]*n2[1]))
...
(x^(v[K]*n1[K])+x^(v[K]*(n1[K]+1))+x^(v[K]*(n1[K]+2))+...+x^(v[K]*n2[K]))
K对应具体问题中物品的种类数。
v[i]表示该乘积表达式第i个因子的权重,对应于具体问题的每个物品的价值或者权重。
n1[i]表示该乘积表达式第i个因子的起始系数,对应于具体问题中的每个物品的最少个数,即最少要取多少个。
n2[i]表示该乘积表达式第i个因子的终止系数,对应于具体问题中的每个物品的最多个数,即最多要取多少个。
已经明白这个多项式展开后m阶的系数,那么只要把它求出来就行了。然后就是套板子了……逃
数组t是存一个中间结果的,数组v是逐项展开后的没项系数。
#include
using namespace std;
const int maxn=1e3+5;
int n,m;
int a[maxn];
int b[maxn];
int t[maxn];
int v[maxn];
int main()
{
while(~scanf("%d%d",&n,&m))
{
memset(v,0,sizeof v);//水果的权重
memset(t,0,sizeof t);
for(int i=1;i<=n;++i)
scanf("%d%d",&a[i],&b[i]);///每种水果的下界,上界
for(int i=a[1];i<=b[1];++i)
v[i]=1;///第1个多项式的系数
for(int i=2;i<=n;++i)///n个多项式的乘积
{
for(int j=0;j<=m;++j)///一个多项式的每一项,最多取m个,所以当指数超过了m就不取
for(int k=a[i];k+j<=m&&k<=b[i];++k)///a[i]是第一项的指数(取东西的下限)
///b[i]是最后一项的指数(取东西的上限) ///这个上限不能超过题目的上限m
t[j+k]+=v[j];
///原来的系数 加上 当前项的系数 [j+k] 表示取出j+k个方案数 类比a^(i+j) = a^i * a^j;
for(int j=0;j<=m;++j)
{
v[j]=t[j];
t[j]=0;///每个多项式算完之后清空
}
}
cout<