组合数学--母函数

最新推荐文章于 2021-05-19 09:34:59 发布
最新推荐文章于 2021-05-19 09:34:59 发布
阅读量116 收藏
点赞数
原文链接：http://www.cnblogs.com/lin1874/p/11382708.html
版权

模板：（转载：https://blog.csdn.net/xiaofei_it/article/details/17042651）

//为计算结果，b为中间结果。
int a[MAX],b[MAX];
//初始化a
memset(a,0,sizeof(a));
a[0]=1;
for (int i=1;i<=17;i++)//循环每个因子
{
    memset(b,0,sizeof(b));
    for (int j=n1[i];j<=n2[i]&&j*v[i]<=P;j++)//循环每个因子的每一项
        for (int k=0;k+j*v[i]<=P;k++)//循环a的每个项
            b[k+j*v[i]]+=a[k];//把结果加到对应位
    memcpy(a,b,sizeof(b));//b赋值给a
}
//模板二
//初始化a，因为有last，所以这里无需初始化其他位
a[0]=1;
int last=0;
for (int i=0;i<K;i++)
{
    int last2=min(last+n[i]*v[i],P);//计算下一次的last
    memset(b,0,sizeof(int)*(last2+1));//只清空b[0..last2]
    for (int j=n1[i];j<=n2[i]&&j*v[i]<=last2;j++)//这里是last2
        for (int k=0;k<=last&&k+j*v[i]<=last2;k++)//这里一个是last，一个是last2
            b[k+j*v[i]]+=a[k];
    memcpy(a,b,sizeof(int)*(last2+1));//b赋值给a，只赋值0..last2
    last=last2;//更新last
}

HDU：1085

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e5+10;
const int inf=0x3f3f3f3f;
int a[maxn],b[maxn];
int n2[maxn];
int v[maxn]={0,1,2,5};
int n1[]={0,0,0,0};
int p,K;
void muhs2(){
    a[0]=1;
    int last=0;
    for(int i=1;i<=3;i++){
        int last2=last+n2[i]*v[i];
        memset(b,0,sizeof(int)*(last2+1));
        for(int j=n1[i];j<=n2[i];j++){
            for(int k=0;k<=last;k++){
                b[k+j*v[i]]+=a[k];
            }
        }
        memcpy(a,b,sizeof(int)*(last2+1));
        last=last2;
    }
}
int main(){
    while(cin>>n2[1]>>n2[2]>>n2[3]){
        if(n2[1]+n2[2]+n2[3]==0)break;
        memset(a,0,sizeof(a));
        muhs2();
        int i=1;
        while(a[i]){
            i++;
        }
        cout<<i<<endl;
    }

    return 0;
}
View Code

 HDU:1171

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=2e6+10;
const int inf=0x3f3f3f3f;
int a[maxn],b[maxn];
int n2[maxn];
int v[maxn];
int n1[maxn]={0};
int n,p;
void muhs2(){
    a[0]=1;
    int last=0;
    for(int i=1;i<=n;i++){
        int last2=last+n2[i]*v[i];
        memset(b,0,sizeof(int)*(last2+1));
        for(int j=n1[i];j<=n2[i];j++){
            for(int k=0;k<=last;k++){
                b[k+j*v[i]]+=a[k];
            }
        }
        memcpy(a,b,sizeof(int)*(last2+1));
        last=last2;
    }
    int i=last/2;
    while(a[i]==0&&i>=0)i--;
    cout<<last-i<<' '<<i<<endl;
}
int main(){
    while(cin>>n){
        if(n<=0)break;
        memset(a,0,sizeof(a));
        for(int i=1;i<=n;i++){
            cin>>v[i]>>n2[i];
        }
        muhs2();
    }

    return 0;
}
View Code

HDU:1398

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e4+10;
const int inf=0x3f3f3f3f;
int a[maxn],b[maxn];
int n2[maxn];
int v[maxn];
int n1[maxn]={0};
int n,p;
void muhs1(){
    memset(a,0,sizeof(a));
    a[0]=1;
    for(int i=1;i<=17;i++){
        memset(b,0,sizeof(b));
        for(int j=0;j*v[i]<=n;j++){
            for(int k=0;k+j*v[i]<=n;k++){
                b[k+j*v[i]]+=a[k];
            }
        }
        memcpy(a,b,sizeof(b));
    }
    cout<<a[n]<<endl;
}
void init(){
    for(int i=1;i<=17;i++)v[i]=i*i;
}
int main(){
    init();
    while(cin>>n){
        if(n==0)break;
        memset(a,0,sizeof(a));
        muhs1();
    }

    return 0;
}
View Code

HDU:2079

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e4+10;
const int inf=0x3f3f3f3f;
int a[maxn],b[maxn];
int n2[maxn];
int v[maxn];
int n1[maxn]={0};
int n,m,p;
void muhs1(){
    memset(a,0,sizeof(a));
    a[0]=1;
    for(int i=1;i<=m;i++){
        memset(b,0,sizeof(b));
        for(int j=0;j*v[i]<=n&&j<=n2[i];j++){
            for(int k=0;k+j*v[i]<=n;k++){
                b[k+j*v[i]]+=a[k];
            }
        }
        memcpy(a,b,sizeof(b));
    }
    cout<<a[n]<<endl;
}
int main(){
    int T;
    cin>>T;
    while(T--){
        cin>>n>>m;
        for(int i=1;i<=m;i++){
            cin>>v[i]>>n2[i];
        }
        muhs1();
    }

    return 0;
}
View Code

HDU:2082

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e4+10;
const int inf=0x3f3f3f3f;
int a[maxn],b[maxn];
int n2[maxn];
int v[maxn];
int n1[maxn]={0};
int n,m,p=50;
void muhs1(){
    memset(a,0,sizeof(a));
    a[0]=1;
    for(int i=1;i<=26;i++){
        memset(b,0,sizeof(b));
        for(int j=0;j*v[i]<=p&&j<=n2[i];j++){
            for(int k=0;k+j*v[i]<=p;k++){
                b[k+j*v[i]]+=a[k];
            }
        }
        memcpy(a,b,sizeof(b));
    }
}
int main(){
    int T;
    cin>>T;
    while(T--){
        for(int i=1;i<=26;i++){
            v[i]=i;
            cin>>n2[i];
        }
        muhs1();
        ll ans=0;
        for(int i=1;i<=50;i++){
            ans+=a[i]*1LL;
        }
        cout<<ans<<endl;
    }

    return 0;
}
View Code

HDU:2110

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e5+10;
const int inf=0x3f3f3f3f;
int a[maxn],b[maxn];
int n2[maxn];
int v[maxn];
int n1[maxn]={0};
int n,m,p=50,sum;
void muhs1(){
    a[0]=1;
    int last=0;
    for(int i=1;i<=n;i++){
        int last2=min(last+n2[i]*v[i],sum);
        memset(b,0,sizeof(int)*(last2+1));
        for(int j=0;j<=n2[i]&&j*v[i]<=last2;j++){
            for(int k=0;k<=last&&k+j*v[i]<=last2;k++){
                b[k+j*v[i]]+=a[k];
                b[k+j*v[i]]%=10000;
            }
        }
        memcpy(a,b,sizeof(int)*(last2+1));
        last=last2;
    }
    if(a[sum]>0)cout<<a[sum]%10000<<endl;
    else cout<<"sorry"<<endl;
}
int main(){
    while(cin>>n){
        if(n==0)break;
        memset(a,0,sizeof(a));
        sum=0;
        for(int i=1;i<=n;i++){
            cin>>v[i]>>n2[i];
            sum+=(v[i]*n2[i]);
        }
        if(sum%3){
            cout<<"sorry"<<endl;
            continue;
        }
        sum/=3;
        muhs1();
    }

    return 0;
}
View Code

HDU:2152

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e6+10;
const int inf=0x3f3f3f3f;
int a[maxn],b[maxn];
int n2[maxn];
int v[maxn];
int n1[maxn];
int n,m,p=50;
void muhs1(){
    a[0]=1;
    int last=0;
    for(int i=1;i<=n;i++){
        int last2=min(last+n2[i]*v[i],m);
        memset(b,0,sizeof(int)*(last2+1));
        for(int j=n1[i];j<=n2[i]&&j*v[i]<=last2;j++){
            for(int k=0;k<=last&&k+j*v[i]<=last2;k++){
                b[k+j*v[i]]+=a[k];
            }
        }
        memcpy(a,b,sizeof(int)*(last2+1));
        last=last2;
    }
    if(last<m)cout<<0<<endl;
    else cout<<a[m]<<endl;
}
int main(){
    while(cin>>n>>m){
        for(int i=1;i<=n;i++){
            v[i]=1;
            cin>>n1[i]>>n2[i];
        }
        muhs1();

    }

    return 0;
}
View Code

 HDU：1284

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e6+10;
const int inf=0x3f3f3f3f;
int a[maxn],b[maxn];
int n2[maxn];
int v[maxn]={0,1,2,3};
int n1[maxn];
int n,m,p=50;
void muhs1(){
    a[0]=1;
    n=32768;
    int last=0;
    for(int i=1;i<=3;i++){
        int last2=n;
        memset(b,0,sizeof(int)*(last2+1));
        for(int j=0;j*v[i]<=last2;j++){
            for(int k=0;k<=last&&k+j*v[i]<=last2;k++){
                b[k+j*v[i]]+=a[k];
            }
        }
        memcpy(a,b,sizeof(int)*(last2+1));
        last=last2;
    }
}
int main(){
    muhs1();
    while(cin>>n){
        cout<<a[n]<<endl;
    }

    return 0;
}
View Code

 HDU：2189

#include <bits/stdc++.h>
using namespace std;
const int maxn=160;
int n,m;
int a[maxn],b[maxn];
int prime[maxn];
void getprime(){
    memset(prime,0,sizeof(prime));
    for(int i=2;i<=150;i++){
        if(!prime[i])prime[++prime[0]]=i;
        for(int j=1;j<=prime[0]&&prime[j]<=150/i;j++){
            prime[prime[j]*i]=1;
            if(i%prime[j]==0)break;
        }
    }
}
void muhs(){
    memset(a,0,sizeof(a));
    for(int i=0;i<=n;i+=2)a[i]=1;
    for(int i=2;i<=prime[0];i++){
        memset(b,0,sizeof(b));
        for(int j=0;j<=150;j++){
            for(int k=0;k+j<=n;k+=prime[i]){
                b[k+j]+=a[j];
            }
        }
        memcpy(a,b,sizeof(b));
    }

}
int main(){
    getprime();
    int T;
    cin>>T;
    while(T--){
        cin>>n;
        muhs();
        cout<<a[n]<<endl;
    }
}
View Code

HDU:2065

指数型母函数，根据泰勒展开（虽然还不懂……）最终答案为2^(n-1)+4^(n-1)

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=160;
const int mod=100;
ll n;
int a[maxn],b[maxn];
ll ksm(ll a,ll b){
    ll ans=1;
    while(b){
        if(b&1){
            ans=(ans*(a%mod))%mod;
        }
        b>>=1;
        a=((a%mod)*(a%mod))%mod;
        ans%=mod;
    }
    return ans%mod;
}

int main(){
    int T;
    while(cin>>T){
        if(T==0)break;
        int flag=1;
        while(T--){
            cin>>n;
            cout<<"Case "<<flag++<<": ";
            cout<<(ksm(2,n-1)+ksm(4,n-1))%mod<<endl;
        }
        cout<<endl;
    }

}
View Code

 POJ：3734

跟上题一样，泰勒展开

#include <iostream>
using namespace std;
typedef long long ll;
const int maxn=160;
const int mod=10007;
ll n;
int a[maxn],b[maxn];
ll ksm(ll a,ll b){
    ll ans=1;
    while(b){
        if(b&1){
            ans=(ans*(a%mod))%mod;
        }
        b>>=1;
        a=((a%mod)*(a%mod))%mod;
        ans%=mod;
    }
    return ans%mod;
}

int main(){
    int T;
    cin>>T;
    while(T--){
        cin>>n;
        cout<<(ksm(2,n-1)+ksm(4,n-1))%mod<<endl;
    }

}
View Code

HDU:1709

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e4+10;
int n;
int a[maxn],b[maxn],v[maxn];
int ans[maxn];
int sum=0;
void muhs(){
    memset(a,0,sizeof(a));
    a[v[1]]=1;
    a[0]=1;
    for(int i=2;i<=n;i++){
        memset(b,0,sizeof(b));
        for(int j=0;j<=sum;j++){
            if(b[j]==1||a[j]==1)b[j]=1;
            if(a[j]==1)b[j+v[i]]=b[abs(j-v[i])]=1;
        }
        memcpy(a,b,sizeof(b));
    }
    /*for(int i=1;i<=sum;i++)cout<<a[i]<<' ';
    cout<<endl;*/
}
int main(){
    ios::sync_with_stdio(false);
    while(cin>>n){
        sum=0;
        for(int i=1;i<=n;i++){
            cin>>v[i];
            sum+=v[i];
        }
        muhs();
        int tot=0;
        for(int i=1;i<=sum;i++){
            if(a[i]==0){
                ans[tot++]=i;
            }
        }
        cout<<tot<<endl;
        for(int i=0;i<tot;i++){
            if(i==tot-1)cout<<ans[i]<<endl;
            else cout<<ans[i]<<' ';
        }
    }

}
View Code

HDU:1261(JAVA)多重集的排列

import java.math.BigInteger;
import java.util.ArrayList;
import java.util.Scanner;

public class Main {

    public static void main(String[] args) {
        Scanner scn = new Scanner(System.in);
        int n;
        int a[]=new int [27];
        while(scn.hasNextInt()){
            n=scn.nextInt();
            if(n==0)break;
            int sum=0;
            for(int i=1;i<=n;i++){
                a[i]=scn.nextInt();
                sum+=a[i];
            }
            BigInteger ans=BigInteger.ONE;
            for(int i=1;i<=sum;i++){
                ans=ans.multiply(BigInteger.valueOf(i));
            }
            BigInteger ansm=BigInteger.ONE;
            for(int i=1;i<=n;i++){
                int x=1;
                for(int j=1;j<=a[i];j++){
                    x*=j;
                }
                ansm=ansm.multiply(BigInteger.valueOf(x));
            }
            System.out.println(ans.divide(ansm));
        }

    }
}
View Code

HDU：1028

#include <bits/stdc++.h>
using namespace std;
const int maxn=1e4+10;
int n;
int a[maxn],b[maxn];
void muhs(){
    memset(a,0,sizeof(a));
    for(int i=0;i<=n;i++)a[i]=1;
    for(int i=2;i<=n;i++){
        memset(b,0,sizeof(b));
        for(int j=0;j<=n;j++){
            for(int k=0;k<=n&&j+k<=n;k+=i){
                b[j+k]+=a[j];
            }
        }
        memcpy(a,b,sizeof(b));
    }
}
int main(){
    while(cin>>n){
        muhs();
        cout<<a[n]<<endl;
    }


    return 0;
}
View Code

 HDU：2069（二维母函数）控制了硬币的数量

#include <bits/stdc++.h>
using namespace std;
const int maxn=310;
int n;
int a[maxn][105],b[maxn][105];
int v[10]={1,5,10,25,50};
void muhs(){
    memset(a,0,sizeof(a));
    for(int i=0;i<=100;i++)a[i][i]=1;
    for(int i=1;i<5;i++){
        memset(b,0,sizeof(b));
        for(int j=0;j<=n;j++){
            for(int k=0;k+j<=n;k+=v[i]){
                for(int kk=0;kk+k/v[i]<=100;kk++){
                    b[k+j][kk+k/v[i]]+=a[j][kk];
                }
            }
        }
        for(int j=0;j<=n;j++){
            for(int k=0;k<=100;k++){
                a[j][k]=b[j][k];
                b[j][k]=0;
            }
        }
    }
}


int main(){
    while(cin>>n){
        muhs();
        int ans=0;
        for(int i=0;i<=100;i++)ans+=a[n][i];
        cout<<ans<<endl;
    }

    return 0;
}
View Code

 

转载于:https://www.cnblogs.com/lin1874/p/11382708.html

weixin_30746117
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
组合数学之母函数
qq_36921652的博客
01-27 241
//老版母函数 //首先用c1，c2两个数组，c1[j]表示j块钱有c1[j]中表达方式 //c1[0]=1,c2 0 //关键for循环 for (int i = 0; i < n; ++i)//该循环是指每*种*货币一个个的放进去，共n种 { for (int j = 0; j <= p; ++j)//这里是c1数组的循环,c1的数组界p { for (int k = 0; k <=
组合数学之母函数问题
weixin_34342905的博客
08-20 191
母函数问题是组合数学中非常经典的问题，大概是本科二年级的课程，非常有意思的一门课，当然也是非常精深的一门课。 定义 对于序列a0，a1，a2。…构造函数G(x)： 则称函数G(x)是序列a0，a1，a2。…的母函数。 非常明显，依据二项展开式。非常easy知道(1+x)^n是序列C(n,0),C(n,1),…,C(n,n)的母函数。假设已知序列a0，a1。a2，...
组合数学(母函数)
林苏泽
03-17 649
母函数生成函数定义一个小栗子3个糖栗子[A 找单词](https://vjudge.net/contest/361385#problem/A)母函数有一点特别注意一下就是初始值的问题，指数为零的初始值设为1，其他都为零，因为你要满是指数为零的系数为1 要不你的多项式相乘，不就没什么意思了吗？，当然最后，0的系数，不加进去，零不算单词另一解法背包方案数[E - Ignatius and the Pr...
组合数学-母函数
weixin_30802171的博客
04-02 88
母函数的定义 母函数计数法则 为什么要(1+x^k)，因为有可能不选质量为k的砝码，不选的情况就是x^0=1 整数拆分！ 转载于:https://www.cnblogs.com/zsben991126/p/10644246.html...
组合数学之 母函数
youngkl博客
11-14 1072
母函数当中经典的问题  正是整数拆分  所谓整数拆分即把整数分解成若干整数的和（相当于把n个无区别的球放到n个无标志的盒子，盒子允许空，也允许放多于一个球）。 整数拆分成若干整数的和，办法不一，不同拆分法的总数叫做拆分数。
组合数学-母函数-与递推关系
09-22
### 组合数学中的母函数与递推关系 #### 一、引言 组合数学是一门研究离散对象的科学，其关注的是如何选择、排列这些对象的方式。母函数是组合数学中的一个重要工具，它可以帮助我们解决许多组合问题，尤其是与...
组合数学- 母函数.rar
09-16
母函数是组合数学中的一种工具，用于求解序列的递推关系，特别是在解决计数问题时非常有用。 母函数（Generating Function）是一个形式幂级数，其系数对应于一个数列中的项。通常，我们将数列的项a_n用母函数A(x)...
组合数学 - 卢开澄
11-04
本书包括：排列与组合，母函数与递推关系，容斥原理与鸽巢原理，贝恩塞特引理与波利亚定理，区组设计与编码，组合算法与复杂性分析。
组合数学－acmer必备
05-24
卢开澄，卢华明，组合数学。ACM必备教程，包含母函数，排列组合，递推，容斥原理，鸽巢原理，组合算法等必备的程序设计知识
组合数学课件-教案
11-27
组合数学是数学的一个重要分支，尤其在计算机科学中扮演着不可或缺的角色。这门学科主要研究在确定模式下事物的配置，关注配置的存在性、计数、分类以及配置的性质。组合数学的应用范围广泛，包括算法设计、数据结构...
组合数学 —— 母函数
Alex_McAvoy的博客
08-04 4527
【概述】 某个序列的母函数是一种形式幂级数，其每一项的系数可以提供关于这个序列的信息。 给定数列，构造一个函数，称 F(x) 为数列的母函数，其中，序列只作为标志用，称为标志函数。 标志函数最重要的形式是，这种情况下的母函数一般形式为： 例如：就是序列 的母函数 也就是说，可以利用来讨论序列的性质，此外还可以引入适当的函数，将问题简化，把复杂的问题变成形式上的初等...
组合数学---母函数（1）
Zsss的博客
09-11 1345
母函数的一部分总结
组合数学5--母函数
weixin_39491014的博客
07-10 943
母函数 一 母函数的定义 1.定义 对于序列C0,C1,C2,…C_0,C_1,C_2,\dotsC0​,C1​,C2​,…构造一函数 G(x)=C0+C1x+C2x2+…G(x)=C_0+C_1x+C_2x^2+\dotsG(x)=C0​+C1​x+C2​x2+…称G(x)G(x)G(x)为序列C0,C1,C2,…C_0,C_1,C_2,\dotsC0​,C1​,C2​,…的母函数。 例如：(1+x)n(1+x)^n(1+x)n为序列C0,C1,C2,…CnC_0,C_1,C_2,\dots C_nC0​
母函数（排列组合问题）
qq_52384627的博客
05-19 304
母函数 一.母函数的定义 二.普通型母函数 例一 例二 例三 例四 例五 例六 三.指数型母函数 定义 例七 例八
组合数学之母函数一（卡特兰数）
Light2Chasers的博客
08-27 1007
卡特兰数公式： 2.一个栈（无穷大）的进栈序列为1，2，3，… ，n,有多少个不同出栈序列？ 分析： （1）对于每个数来说，必须进栈一次，出栈一次。我们把进栈设为状态“1”，出栈设为状态“0”。n个数的所有状态对应n个1和n个0组成的2n位二进制数。由于等待入栈的操作数按照1，…，n的顺讯排列，入栈的操作数b大于等于出栈的操作数a(a&amp;amp;lt;=b)，因此输出序列的总数目等于由左到右扫...
数理方法习题 前六章.pdf
最新发布
08-03
数理方法习题 前六章
ACM-母函数及其应用
这个概念主要源自于组合数学，它将一个序列的系数转换成一个多项式，从而通过多项式的运算来分析序列的性质。 母函数的基本思想是将序列的每一项与多项式的一个项相对应。例如，序列a0, a1, a2, ... 可以对应于...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值