又上了节课。。。俩题

计算系数    组合数问题。。。

要不是大佬指点就只能阶乘暴力算了

(主要还是我忘了杨辉三角)

杨辉三角与组合数C有着千丝万缕的联系,在计算,使用方面相当方便。

先说计算系数

计算系数

给定一个多项式(by+ax)^k,请求出多项式展开后x^n*y^m项的系数。

输出结果%10007

#include<bits/stdc++.h>
using namespace std;
long long k,a,b;
long long n,m;
long long t=;
long long c[][];
int main(){
cin>>a>>b>>k>>m>>n;
c[][]=;
a=a%;
b=b%;
k=k%;
n=n%;
m=m%;//输入需要先取模,不然会很大
for(int i=;i<=k;i++){
for(int j=;j<=i;j++){
c[i][]=;
c[i][j]=c[i-][j]+c[i-][j-];//处理原理,杨辉三角
c[i][j]%=;//预处理,待会直接调用
}
}
for(int i=;i<=m;i++){
t*=a;
t%=;
}
for(int i=;i<=n;i++){
t*=b;
t%=;
}//忘记*a*b很致命...
cout<<(t%*c[k][n]%)%;
  return ;
}

组合数问题

给定 n,m 和 k,对于所有的  0≤i≤n,0≤j≤min(i,m) 有多少对 (i,j) 满足C​ 是 k的倍数。

#include<bits/stdc++.h>
using namespace std;
int k;
int n,m;
int t;
int c[][];
int ans[][];
int main(){
scanf("%d%d",&t,&k);
for(int i=;i<=;i++) c[i][]=;
for(int i=;i<=;i++){
for(int j=;j<=i;j++){
c[i][j]=(c[i-][j]+c[i-][j-])%k;//这里取模免去以后的判断
}
}
for(int i=;i<=;i++){
for(int j=;j<=i;j++){
ans[i][j]=ans[i-][j]+ans[i][j-]-ans[i-][j-];//前缀和,减去重复的
if(!c[i][j]) ans[i][j]++;
}
ans[i][i+]=ans[i][i];
}
for(int a=;a<=t;a++){
scanf("%d%d",&n,&m);
if(n<m) m=n;
printf("%d\n",ans[n][m]);
}
return ;
}

同样杨辉三角

最后%下子谦。

OTZ%%%子谦。大佬的更多相关文章

  1. 简单数据结构题(from 钟子谦——IOI2018集训队自选题)

    简单数据结构题(from 钟子谦--IOI2018集训队自选题) 试题描述 给一棵 \(n\) 个点的树,点权开始为 \(0\) ,有 \(q\) 次操作,每次操作是选择一个点,把周围一圈点点权 \( ...

  2. test20191210 钟子谦

    100+40+0=140.暴力没写满-- 简单模拟 很久很久以前,有一个 \(1\sim n\) 的排列 \(a\),还有一个长度为 \(q\) 的,每个元素在 \(1\) 到 \(n\) 之间的序列 ...

  3. 【洛谷P1313 计算系数】

    题目连接 #include<algorithm> #include<iostream> #include<cstring> #include<cstdio&g ...

  4. 【洛谷P2822 组合数问题】

    题目连接 #include<iostream> #include<cstring> #include<cstdio> #include<cctype> ...

  5. noip2018——题解&总结

    近期正在疯狂复习某些东西,这篇博客尽量年底更完……(Day2T2除外) 好了,所有的希望都破灭了,原来这就是出题人的素质.——一个被欺骗的可怜 $OIer$ 人生中倒数第三次 $noip$ (Mayb ...

  6. Eva 剧情解析

    Eva 剧情解析 来源 https://zhuanlan.zhihu.com/p/20864898 [0.写在前面的话] 相信和我年龄差不多的小伙伴们对<新世纪福音战士>( <Neo ...

  7. ZJOI2017 day1滚粗记

    这几天去温州作为外省选手参加了$ZJOI day1$.打了几天的酱油,考试也滚粗了.. $day -2$ 中午从学校出发,坐飞机去温州.到了温州以后吃完晚饭就回宾馆.把一直想做的糖果公园做完了以后就堕 ...

  8. 45.更新一下scrapy爬取工商信息爬虫代码

    这里是完整的工商信息采集代码,不过此程序需要配合代理ip软件使用.问题:1.网站对ip之前没做限制,但是采集了一段时间就被检测到设置了反爬,每个ip只能访问十多次左右就被限制访问.2.网站对请求头的检 ...

  9. SDOI2019Round1游记

    SDOI2019Round1游记 Day 0 报道日,早晨五点睡的觉,一觉醒来已经一点半了,收拾收拾东西报道去了.因为没吃饭,坐着出租车晕车了,我让师傅把我放到历下大润发,我去金拱门吃了点饭才去的23 ...

随机推荐

  1. Fundebug前端JavaScript插件更新至1.7.1,拆分录屏代码,还原部分Script error.

    摘要: BUG监控插件压缩至18K. 1.7.1拆分了录屏代码,BUG监控插件压缩至18K,另外我们还原了部分Script error,帮助用户更方便地Debug.请大家及时更新哈~ 拆分录屏代码 从 ...

  2. laravel5集成支付宝alipay扫码支付流程(Laravel 支付解决方案)

    首先我们来探讨如何在Laravel应用中使用支付宝进行支付,对此,GitHub上有很多相关的包,其中最流行的两个包:Omnipay For Laravel 5 & Lumen 和 Larave ...

  3. django入门与实践 - 关于升级到django 3.7,三种模板超链接配置(编辑中)

    第一种方法: 在myblog/urls.py模块中: from django.contrib import admin from django.urls import path, include ur ...

  4. C++ 11 创建和使用共享 weak_ptr

    1.为什么需要weak_ptr? 在正式介绍weak_ptr之前,我们先来回忆一下shared_ptr的一些知识.我们知道shared_ptr是采用引用计数的智能指针,多个shared_ptr实例可以 ...

  5. win10系统关闭自动更新

    win10关闭自动更新 步骤①右键“此电脑”选择“管理”选项 步骤②(如下图所示): 步骤③:     步骤④: 好啦!这样就大功告成了!

  6. java 中文乱码以及转码

    查看此文章需要对字符集编码有一定的认识:任意门:字符集编码基础 一.字符串的内部表示? 重点:字符串在java(指在JVM中.在内存中)中统一用unicode表示( 即utf-16 LE) , 下面解 ...

  7. 初级c++编码规范

        想了很久,第一篇文章还是应该写编码规范好一点.编码规范是一个仁者见仁的问题,为了避免复杂庞大,自己总结了一套简单版本的规范. 简介     本文介绍一份自己使用的C++编码规范.第一次正式进入 ...

  8. Yii2.0调用sql server存储过程并获取返回值

    1.首先展示创建sql server存储过程的语句,创建一个简单的存储过程,测试用. SET ANSI_NULLS ON GO SET QUOTED_IDENTIFIER ON GO CREATE P ...

  9. 【原】Java学习笔记005 - 流程控制

    package cn.temptation; public class Sample01 { public static void main(String[] args) { // 程序的流程控制(流 ...

  10. CentOS 6.2 中文

    在虚拟机里面安装好centos6.2之后,默认是英文!   对于命令行操作无所谓啦,但是如果想看界面,就不是很适应!   修改方法如下:   1.用root登录系统,密码为创建虚拟机时候的密码.创建虚 ...