【洛谷2822】组合数问题(组合数的递推公式和杨辉三角)

题目地址:https://www.luogu.org/problemnew/show/P2822

解题思路:


组合数可能会很大,所以要对k取模,模为0则是k的倍数,组合数实际上是杨辉三角(c[i][j]=c[i-1][j]+c[i-1][j-1],两个的计算方法都是这个式子)

数组c[i][j]即C(i,j),ans[i][j]存答案,先打表

方法1:

ans[i][j]=ans[i-1][j]-ans[i-1][j-1]+ans[i][j-1];
if(!c[i][j]&&j<=i) ans[i][j]++;

注意:这个时候i,j都是从1~maxn遍历的,因为在求ans[i][j]时,上一行对应的ans[i-1][j]不能为0,否则-a[i-1][j-1]就可能是负数,结果错误

方法2:

只需要根据输入的n的最大值打表,不要建maxn*maxn大小的表
if(!c[i][j])
    ans[i][j]=ans[i-1][j]+1;
else ans[i][j]=ans[i-1][j];

对于组合数数组c[][]的存储有两种方法:只是控制输出的条件不太一样。。搞了好久才知道区别在哪。。太弱了

<figcaption> i和j都从1开始存,实际上就是存了个杨辉三角,第一列都是C(i,0) </figcaption>
<figcaption> i和j都从1开始存,不存杨辉三角的第一列,此时的第一列都是C(i,1) </figcaption>

ac代码:


方法1

#include <iostream>
#include <cmath>
#include <cstring>
#include <cctype>
#include <algorithm>
#include <stdlib.h>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <sstream>
#define maxn 2005
using namespace std;
int t,k,n,m;//n>=m
long long  c[maxn][maxn]={0},ans[maxn][maxn]={0};//求组合数的数组和答案数组
void pre_solve()
{
    c[0][0]=1;
    for(int i=1;i<=maxn;i++)
        c[i][i]=c[i][0]=1;
    for(int i=1;i<maxn;i++)
    {
        for(int j=1;j<=i/2;j++)
        {
            c[i][j]=(c[i-1][j]+c[i-1][j-1])%k;
            c[i][i-j]=c[i][j];
        }
    }
    for(int i=1;i<maxn;i++)
    {
        for(int j=1;j<maxn;j++)
        {
            ans[i][j]=ans[i-1][j]-ans[i-1][j-1]+ans[i][j-1];
            if(!c[i][j]&&j<=i) ans[i][j]++;
            //cout<<i<<" "<<j<<" "<<ans[i][j]<<endl;
        }
    }

}
int main()
{
    //freopen("/Users/zhangkanqi/Desktop/11.txt","r",stdin);
    scanf("%d %d",&t,&k);
    pre_solve();
    while(t--)
    {
        scanf("%d %d",&n,&m);
        printf("%lld\n",ans[n][m]);
    }
    return 0;
}

方法2:

i和j从1开始存,实际上存了个杨辉三角

#include <iostream>
#include <cmath>
#include <cstring>
#include <cctype>
#include <algorithm>
#include <stdlib.h>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <sstream>
#define maxn 2005
using namespace std;
int t,k,n,m;//n>=m
long long  c[maxn][maxn]={0},ans[maxn][maxn]={0};//求组合数的数组和答案数组
int maxnn=0,maxm=0;
struct node{
    int n,m;
}a[10001];
int main()
{
    //freopen("/Users/zhangkanqi/Desktop/11.txt","r",stdin);
    scanf("%d %d",&t,&k);
    for(int i=0;i<t;i++)
    {
        scanf("%d %d",&a[i].n,&a[i].m);
        maxnn=max(maxnn,a[i].n);
        maxm=max(maxm,a[i].m);
    }
    maxnn++;//!!
    for(int i=1;i<=maxnn;i++)
    {
        c[i][i]=c[i][1]=1;
    }
    for(int i=2;i<=maxnn;i++)
        for(int j=2;j<=i;j++)
         c[i][j]=(c[i-1][j]+c[i-1][j-1])%k;
    /*for(int i=1;i<=maxnn;i++) {
        for (int j = 1; j <= maxnn; j++) {
            cout << c[i][j] << " ";
        }
        cout << endl;
    }//输出杨辉三角*/
    for(int i=1;i<=maxnn;i++)
        for(int j=1;j<=i;j++)
        {
            if(!c[i][j])
                ans[i][j]=ans[i-1][j]+1;
            else ans[i][j]=ans[i-1][j];
            //cout<<i<<" "<<j<<" "<<ans[i][j]<<endl;
        }
    for(int i=0;i<t;i++)
    {
        long long res=0;
        for(int j=0;j<=a[i].m+1;j++)//a[i].m+1
        {
            res+=ans[a[i].n+1][j];//a[i].n+1
        }
        printf("%lld\n",res);
    }
    return 0;
}

i和j从1开始存,不存杨辉三角的第一列

#include <iostream>
#include <cmath>
#include <cstring>
#include <cctype>
#include <algorithm>
#include <stdlib.h>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <sstream>
#define maxn 2005
using namespace std;
int t,k,n,m;//n>=m
long long  c[maxn][maxn]={0},ans[maxn][maxn]={0};//求组合数的数组和答案数组
int maxnn=0,maxm=0;
struct node{
    int n,m;
}a[10001];
int main()
{
    freopen("/Users/zhangkanqi/Desktop/11.txt","r",stdin);
    scanf("%d %d",&t,&k);
    for(int i=0;i<t;i++)
    {
        scanf("%d %d",&a[i].n,&a[i].m);
        maxnn=max(maxnn,a[i].n);
        maxm=max(maxm,a[i].m);
    }
    for(int i=1;i<=maxnn;i++)
    {
        c[i][1]=i%k;//注意这句!
        c[i][i]=1;
    }
    for(int i=2;i<=maxnn;i++)
        for(int j=2;j<i;j++)
            c[i][j]=(c[i-1][j]+c[i-1][j-1])%k;
    for(int i=1;i<=maxnn;i++)
        for(int j=1;j<=i;j++)
        {
            if(!c[i][j])
                ans[i][j]=ans[i-1][j]+1;
            else ans[i][j]=ans[i-1][j];
        }
    for(int i=0;i<t;i++)
    {
        long long res=0;
        for(int j=0;j<=a[i].m;j++)//a[i].m
        {
            res+=ans[a[i].n][j];//a[i].n
        }
        printf("%lld\n",res);
    }
    return 0;
}

 

全部评论

相关推荐

09-26 19:45
门头沟学院 Java
点赞 评论 收藏
分享
emmm别问我为啥上一条帖子隔了两个月我才开始投简历和拿offer,因为我懒😰简单流程如下:周一凌晨改好的简历,然后到处乱投简历;周二接到了三维家的一面通知,临时抱佛脚的背了一些八股;周三上午一面下午通知第二天hr面;周四上午hr面下午拿offer,遂收手支线:在BOSS上顺手投了几个大厂,投字节的时候不小心投城客户端了,结果过了一天HR突然把我简历要走了,还问我能不能整客户端,我直接一口答应(脏面评警告😢)结果在周三下午的时候给我打电话,说前端有空缺实习岗,问我有没有兴趣,然后就跟我约了周四下午一面😰我都没咋准备啊,咩都不会啊😭结果周四下午面完,晚上打电话通知过一面了,赶紧把二面约在下周一下午,留点缓冲时间。逆大天了,我一半的问题都不会,他居然给我过了?运气未免有点好了😥现在正在恶补计网、网安、性能优化的东西(这三大板块我是几乎一点不会,一面几乎一点答不出来,加上我又没怎么背八股,这块被干烂了😵)心得体会与经验:1.&nbsp;我giao怎么这么快就结束了,我还以为要找好久😨2.&nbsp;大厂的面试问题真的和中厂小厂很大不同,比如在三维家我能自己吹水到vue的数据劫持、Proxy代理响应式之类的他们就觉得很不错了,但是在字节你但凡敢提到一下就会追问你细节了,一追问马脚就全漏出来了3.&nbsp;有信心真的很重要,我感觉我能拿中厂offer最重要的就是吹水吹出自信来了,以至于三维家面试反问面试官有哪里还需要改进的时候,他就说很不错了解的很多😦4.&nbsp;理解很重要,我从头到尾真没背过很多八股,不过有一些知识确实是敲过代码验证过,所以面试的时候能吹水吹得出来😇想了解面经啥的可以直接评论区问我,但我可能也说不全,因为我没有记录,而且今天摆了一天感觉记忆快清空了😵下面是故事时间:我暑假刚开始的时候才开始准备八股,印象很深那个时候连什么原型、事件循环、闭包这些名词都没听过,资料也不知道怎么找,就一直零零散散的准备,感觉也只有js稍微背了一下八股,其他很多时候都是靠完全理解和手写熟悉一些机制的,但这样做效率很低,反正准备了一个多星期半个月就开摆了😭结果一摆就摆到了开学,笔记是乱七八糟的,八股是忘光光的,简历是一直没改的,实习也是一直没投过的。直到上周日晚上偶然和师兄聊天,他突然问我“你怎么还不找实习”,那天晚上才幡然醒悟,是时候做点事情了😡然后就按照上面描述的来走了。其实我感觉我从头到尾都没背特别多八股,也没怎么找刷题资料啥的,早期就是翻尚硅谷或者黑马的入门视频从头学起,中期用面试鸭看了一点点题,主要是在学js机制和敲js代码,后期才发现了w3c的面经网站,然后在那里看着学(那个时候已经懒得敲了,因为有些问题与代码感觉不像是给找实习的看的,忒细了点😂)接下来继续准备字节二面吧,虽然几乎没啥可能可以通过,但是万一有奇迹呢?😍😍😍也祝大家能够早日拿到心仪的offer
我的offer呢😡:我已经预见10天后你会发,节孝子启动了
投递三维家等公司10个岗位
点赞 评论 收藏
分享
评论
点赞
收藏
分享

创作者周榜

更多
牛客网
牛客网在线编程
牛客网题解
牛客企业服务