这道题虽然不是计数问题,但同样可以用母函数解决,好在这道题的时间要求较宽。
代码如下:
#include
#include
#include
#include
typedef long long ll;
using namespace std;
ll c1[250005],c2[250005];
struct fac
{
ll v;
ll n;
}facs[55];
bool cmp(fac x,fac y)
{
return x.v>N)&&(N>0))
{
ma=0;
for(i=1;i<=N;i++)
{
cin>>facs[i].v>>facs[i].n;
ma+=facs[i].v*facs[i].n; //ma记录所有facilities的总值
}
sort(facs+1,facs+N+1,cmp); //从小到大排序,按v
for(i=0;i<=ma;i++)
c1[i]=c2[i]=0;
for(i=0;i<=facs[1].n;i++) //母函数第一个多项式系数初始化为1
c1[i*facs[1].v]=1;
for(i=2;i<=N;i++)
{
for(j=0;j<=ma;j++)
for(k=0;(k+j<=ma)&&(k<=facs[i].v*facs[i].n);k+=facs[i].v)
c2[j+k]=c2[j+k]+c1[j];
for(j=0;j<=ma;j++)
{
c1[j]=c2[j];
c2[j]=0;
}
}
/*for(i=1;i<=N;i++)
{
for(j=1;j<=facs[i].n;j++)
cout<=B) //因为i是从小到大,所以如果发现A>=B,就可以输出A,B,直接结束循环了
{
cout<