题解 P5239 【回忆京都】
你们这些写题解的,就不能把话说清楚嘛!(吐槽1)
你们这些出题的,就不能多出点东方嘛!(吐槽2)
你们这些做题的,就不来写一篇详细一点的题解嘛!(吐槽3)
以上均是个人吐槽,纯属吐槽,不带任何针对性和感情色彩。
声明:
本题解适宜蒟蒻(比如我等)观看,若卡关,可以来此题解领提示。
小金羊写的题解致力于让刚刚学习二维数组的同学都能明白!!
把我顶上去让像我一样的juruo明白一下
回到正题。
首先还是看看这个题咋推出来的和杨辉三角&二位前缀和有关系?
我自己用python 3先手推了一下组合数,
先上python 3推组合数的代码:
import os
def jc(num):
if num is 0 or num is 1:
return 1
else :
return num*jc(num-1)
#阶乘递归版,适用于自己造的小型数据
#不是机惨w......
def zhs(i,j):
if i > j :
return 0
else :
return jc(j)/(jc(i)*jc(j-i))
#组合数配合阶乘递归,适用于小型数据推算
n=int(input())
while n is not 0:
n-=1
L=list(input().split())
print(zhs(int(L[0]),int(L[1])))
os.system("pause")
'''
print("Exit?",end=' ')
if input() is 'Yes':exit(1)
else :exit(1)
'''
一个个试一下,然后就发现这样一个鬼畜事件:
打眼一看:好像跟杨辉三角有那么点联系......
哦!缺了第一列的杨辉三角!
然后杨辉三角创造方法:\(O(1000^2\div 2)\)递推打表!
公式:
\]
别忘了取模(废话)
等等......
这个和此题有什么联系吗?
这个题让求组合数的和。
求和,先把区间用yellow色画出来。
然后发现......
(下图中填充黄色的是求和区域,紫色是和)
你发现了吗?这是一个二重的杨辉三角!
其实就是一个二维的前缀和预处理工作。
没有事情干的同学,树状数组&线段树都可解决这个问题。
根本就是两次递推打表!
然后这个题第一次我竟然只有10分......
但是要追求一个最优的方法。
这个时候我们要改一下原先的变量定义,是关于ans[i][j]
方面:
设ans[i][j]
表示的是前i
列前j
行的前缀和。
推出这个题前缀和公式的过程:
杨辉三角到底有什么好处?
其实杨辉三角给你预处理了单列上的前缀和。
然后单列上前缀和用汉语表示就是:
杨辉三角形第i
列上前j
行の前缀和就是杨辉三角形第i+1
行第j+1
列上的数据。
数学公式?(第一个中括号内暂定是1)
\]
根据以上推论,得出很多列的(就是二维的)前缀和推论。
汉语表达:
前i+1
列前j+1
行的前缀和就等于前i+1
列前j
行的前缀和加上前i+1
列第j+1
行的前缀和(即杨辉三角形第i+2
行第j+2
列那一项)。
公式?
\]
这样避免了许多不必要的记公式过程......
Upd4 2019/3/6:
有同学问,为啥是给你预处理了单列上的前缀和?
我们根据\(yh[i][j]=yh[i-1][j-1]+yh[i-1][j]\),
那么又\(\because yh[i-1][j]=yh[i-2][j-1]+yh[i-1][j]\)......
以此类推,得到yh[i+1][j+1]
,相当于我们得到了杨辉三角第j
列的前i
个数的和+\(yh[1][j]\)。
且根据我们杨辉三角靠左排列放置的方式,\(yh[1][j](j>0)\)必定为0。(观察可知\(yh[0][0]\)的右侧即\(yh[0][1]=0\),而右侧\(yh[0][1]\)和右侧数据都是0)
于是我们得出递推公式,ans[i][j]=ans[i][j-1]+yh[i+1][j+1]
。
到这里,我们找到了一个非常完美的没有过多数据+-的操作。
(qwq比cz dalao的算法的常数小)
坑点:
1.你以为取了mod就不会爆负数吗?太天真啦!
\(\therefore ans=(ans+mod)\)%\(mod\)
2.你以为我会\(O(2\cdot 1000^2)\)做吗?太天真啦!
复杂度\(O(1000^2\div 2+1000^2+query)\)
Upd1 2019/3/3:
关于代码和蒟蒻的二维前缀和求法补充完善
代码来辣!
求前缀和还是预处理吧......线段树什么的玩不来......
注意下面代码,求前缀和的时候意义和求杨辉三角的时候有所不同。
(原因见上面)
#include <iostream>
#include <cstdio>
#include <map>
using namespace std;
typedef long long int lli;
const int maxn=1008;
const lli mod=19260817;
lli yh[maxn+1][maxn+1],ans[maxn+1][maxn+1];
int n,m,q;
void Init()
{
for (register int i=0;i<=1004;i++)
{
yh[i][i]=1;
}
for (register int i=0;i<=1004;i++)
{
yh[i][0]=1;
}
for (register int i=2;i<=1004;i++)
{
for (register int j=1;j<=i;j++)
{
yh[i][j]=(yh[i-1][j-1]+yh[i-1][j]+mod)%mod;
}
}
//杨辉三角形的生成方式
for (register int i=1;i<=1004;i++)
{//前i列
for (register int j=1;j<=1004;j++)
{//前j行
ans[i][j]=(ans[i][j-1]+yh[i+1][j+1]+mod)%mod;
}
}
//二维前缀和的生成方式
}
int main()
{
Init();
scanf("%d",&q);
while (q--)
{
scanf("%d%d",&n,&m);
printf("%lld\n",ans[m][n]);
//注意ans[][]的定义!!
}
return 0;
}
Upd2 2019/3/3:
重新更正代码,实际4-WA
,现在AC。
原因在于这个算法的局限性:
实际上需要推到1000+,时间复杂度上虽然小了,但是容易边界数据卡没了......
实际上我就是这样4-WA:read 0
的......
题解 P5239 【回忆京都】的更多相关文章
- P5239 回忆京都
题目地址:P5239 回忆京都 杨辉三角即组合数的"打表"形式 再求一个二维前缀和 然后处理一下负数即可(因为在求前缀和的过程中有减法) #include <bits/std ...
- 洛谷 P5239 回忆京都 题解
题面 裸的杨辉三角前缀和,但----- 在求前缀和的时候有可能得到一个负数(由于取模的原因),所以一定要加上模数后再取模!!!! #include <bits/stdc++.h> #def ...
- P5239 回忆京都(洛谷3月月赛T2)
题目描述 射命丸文在取材中发现了一个好玩的东西,叫做组合数. 组合数的定义如下:从n个不同元素中,任取m(m≤n)个元素并成一组,叫做从n个不同元素中取出m个元素的一个组合.所有组合的数量,就是组合数 ...
- 洛谷P5239 回忆京都
和 NOIP2016TG 组合数问题 差不多是一样的-- 首先要知道杨辉三角和组合数之间的关系 看一下数据范围,很明显要避免重复计算,而且查询的复杂度要非常小 一看n, m <= 1000 这明 ...
- HBSX2019 3月训练
Day 1 3月有31天废话 今天先颓过了就只剩30天了 初步计划 每天一道字符串/数据结构题 图论学习 根据<若干图论模型探讨>(lyd)复习 二分图与网络流学习 <算法竞赛进阶指 ...
- 5239-回忆京都-洛谷3月赛gg祭
传送门 题目背景 第十五届东方人气投票 音乐部门 106名 第四次国内不知道东方的人对东方原曲的投票调查 51名 回忆京都副歌我tm吹爆,东方文花帖我tm吹爆! 题目描述 射命丸文在取材中发现了一个好 ...
- Codeforces Round #373 (Div. 2) E. Sasha and Array 线段树维护矩阵
E. Sasha and Array 题目连接: http://codeforces.com/contest/719/problem/E Description Sasha has an array ...
- 【BZOJ1264】[AHOI2006]基因匹配Match DP+树状数组
[BZOJ1264][AHOI2006]基因匹配Match Description 基因匹配(match) 卡卡昨天晚上做梦梦见他和可可来到了另外一个星球,这个星球上生物的DNA序列由无数种碱基排列而 ...
- 【BZOJ2795】[Poi2012]A Horrible Poem hash
[BZOJ2795][Poi2012]A Horrible Poem Description 给出一个由小写英文字母组成的字符串S,再给出q个询问,要求回答S某个子串的最短循环节.如果字符串B是字符串 ...
随机推荐
- java中extends和implements的区别
implements:接口 1.实现一个接口就是要实现该接口中的所有方法(抽象类除外) 2)接口中的方法都是抽象的 多个无关的类可以实现同一个接口,一个类可以实现多个无关的接口 extends:继承父 ...
- Convolutional Networks for Images,Speech,and Time-series
Convolutional Networks for Images,Speech,and Time-series Yann LeCun Yoshua Bengio 1995年的 1引言 多层BP网络 ...
- Mac软件推荐
其他: AP文档浏览器+代码片段管理工具:Dash 抓包工具:Charles 使用教程:http://www.cnblogs.com/dsxniubility/p/4621314.html 音乐播放 ...
- 20155222卢梓杰 实验一 逆向及Bof基础
实验一 逆向及Bof基础 1.实验对象为32位可执行文件pwn1,这个程序主要有main.foo.getshell这三个函数,其中foo函数功能为输出输入的字符串,getshell函数功能为打开一个s ...
- TPM及TSS协议栈的安装使用
TPM及TSS协议栈的安装 标签: 可信计算. 目录 安装环境介绍 TPM及TSS安装 软件包下载 TPM 安装 安装TSS 安装tpm-tools 交互过程 编写代码测试TPM是否可用 编写代码测试 ...
- 设计模式 笔记 命令模式 Command
//---------------------------15/04/25---------------------------- //Conmmand 命令模式----对象行为型模式 /* 1:意 ...
- win10+anaconda3+python3.6+opencv3.1.0
最近用windows系统比较多,就想在win10下搞一下深度学习这一方面的研究,那么就需要配置好环境巴拉巴拉的一堆东西.默默记个笔记,正所谓“好记性不如烂笔头”. 1.安装Anaconda 这个是一个 ...
- vector 简介
vector简介 vector是STL中最常见的容器,它是一种顺序容器,支持随机访问.vector是一块连续分配的内存,从数据安排的角度来讲,和数组极其相似, 不同的地方就是:数组是静态分配空间,一旦 ...
- 机器学习英雄访谈录之双料 Kaggle 大师:Dr. Jean-Francois Puget
目录 机器学习英雄访谈录之双料 Kaggle 大师:Dr. Jean-Francois Puget 正文 对我的启发 机器学习英雄访谈录之双料 Kaggle 大师:Dr. Jean-Francois ...
- 深度学习目标检测综述推荐之 Xiaogang Wang ISBA 2015
一.INTRODUCTION部分 (1)先根据时间轴讲了历史 (2)常见的基础模型 (3)讲了深度学习的优势 那就是feature learning,而不用人工划分的feature engineeri ...