[gym102220I]Temperature Survey
(为了方便,以下记$a_{0}=0,a_{n+1}=n$,并将$n$加上1)
构造一个$n$行的网格图,从上到下第$i$行有$a_{i}$个格子,格子左对齐
记第$i$行第$j$个格子为$(i,j)$,格子集合$\{(i,j)\mid i_{1}\le i\le i_{2}$且$j_{1}\le j\le j_{2}\}$为$([i_{1},i_{2}],[j_{1},j_{2}])$
此时,考虑一条从$(1,1)$到$(n,a_{n})$的路径(只能向下或向右),令$b_{i}$为路径中从第$i$行到第$i+1$行时的格子编号,不难发现这样的路径与合法的$b_{i}$一一对应,那么不妨统计路径数
关于路径数,显然可以dp解决,即令$f_{i,j}$表示从$(1,1)$到$(i,j)$的路经数,则$f_{i,j}=f_{i,j-1}+f_{i-1,j}$
使用分治优化dp转移,当分治到区间$[l,r]$时,需要根据$f_{[l,r],a_{l-1}}$的值求出$f_{r,(a_{l-1},a_{r}]}$的值
具体的,分治过程如下——
1.令$mid=\lfloor\frac{l+r}{2}\rfloor$,递归左区间$[l,mid]$,根据$f_{[l,mid],a_{l-1}}$的值求出$f_{mid,(a_{l-1},a_{mid}]}$的值
2.根据$f_{(mid,r],a_{l-1}}$和$f_{mid,(a_{l-1},a_{mid}]}$的值,快速求出$f_{(mid,r],a_{mid}}$和$f_{r,(a_{l-1},a_{mid}]}$的值
3.递归右区间$(mid,r]$,根据$f_{(mid,r],a_{mid}}$的值求出$f_{r,(a_{mid},a_{r}])}$的值
(其中,第二步显然可以写成卷积的形式,直接ntt即可)
边界条件:若$l>r$直接退出,若$l=r$令$f_{l,(a_{l-1},a_{r}]}=f_{l,a_{l-1}}$并退出
另外,有一些细节问题:
1.为了避免位置不合法,需要保证$a_{l-1}<a_{l}$,若不满足则不断增加$l$即可
2.关于$f$的存储,只需要用两个数组,分别存储当前每一行/列递归到最右边的一列/行的$f$值即可
总复杂度为$o(n\log^{2}n)$,可以通过


1 #include<bits/stdc++.h>
2 using namespace std;
3 #define N 200005
4 #define mod 998244353
5 #define ll long long
6 #define vi vector<int>
7 vi v,vl,vu;
8 int t,n,fac[N<<1],inv[N<<1],rev[N<<3],a[N],f[N],ans[N];
9 int c(int n,int m){
10 return (ll)fac[n]*inv[m]%mod*inv[n-m]%mod;
11 }
12 int qpow(int n,int m){
13 int s=n,ans=1;
14 while (m){
15 if (m&1)ans=(ll)ans*s%mod;
16 s=(ll)s*s%mod;
17 m>>=1;
18 }
19 return ans;
20 }
21 void ntt(vi &a,int n,int p){
22 for(int i=0;i<(1<<n);i++)
23 if (i<rev[i])swap(a[i],a[rev[i]]);
24 for(int i=2;i<=(1<<n);i<<=1){
25 int s=qpow(3,(mod-1)/i);
26 if (p)s=qpow(s,mod-2);
27 for(int j=0;j<(1<<n);j+=i)
28 for(int k=0,ss=1;k<(i>>1);k++,ss=(ll)ss*s%mod){
29 int x=a[j+k],y=(ll)a[j+k+(i>>1)]*ss%mod;
30 a[j+k]=(x+y)%mod;
31 a[j+k+(i>>1)]=(x+mod-y)%mod;
32 }
33 }
34 if (p){
35 int s=qpow((1<<n),mod-2);
36 for(int i=0;i<(1<<n);i++)a[i]=(ll)a[i]*s%mod;
37 }
38 }
39 vi mul(vi a,vi b){
40 int n=0,m=a.size()+b.size()-1;
41 while ((1<<n)<m)n++;
42 for(int i=0;i<(1<<n);i++)rev[i]=(rev[i>>1]>>1)+((i&1)*(1<<n)/2);
43 while (a.size()<(1<<n))a.push_back(0);
44 while (b.size()<(1<<n))b.push_back(0);
45 ntt(a,n,0);
46 ntt(b,n,0);
47 for(int i=0;i<(1<<n);i++)a[i]=(ll)a[i]*b[i]%mod;
48 ntt(a,n,1);
49 return a;
50 }
51 void calc(int l,int r){
52 while ((l<=r)&&(a[l-1]==a[l]))l++;
53 if (l>r)return;
54 if (l==r){
55 for(int i=a[l-1]+1;i<=a[r];i++)ans[i]=f[l];
56 return;
57 }
58 int mid=(l+r>>1);
59 calc(l,mid);
60 vl.clear(),vu.clear();
61 for(int i=mid+1;i<=r;i++)vl.push_back(f[i]);
62 for(int i=a[l-1]+1;i<=a[mid];i++)vu.push_back(ans[i]);
63
64 v.clear();
65 for(int i=0;i<r-mid;i++)v.push_back(c(i+a[mid]-a[l-1]-1,i));
66 v=mul(v,vl);
67 for(int i=0;i<r-mid;i++)f[i+mid+1]=v[i];
68
69 v.clear();
70 for(int i=0;i<r-mid+a[mid]-a[l-1]-1;i++)v.push_back(fac[i]);
71 for(int i=0;i<r-mid;i++)vl[i]=(ll)vl[i]*inv[r-mid-1-i]%mod;
72 v=mul(v,vl);
73 for(int i=0;i<a[mid]-a[l-1];i++)ans[i+a[l-1]+1]=(ll)v[i+r-mid-1]*inv[i]%mod;
74
75 v.clear();
76 for(int i=0;i<a[mid]-a[l-1];i++)v.push_back(c(i+r-mid-1,i));
77 v=mul(v,vu);
78 for(int i=0;i<a[mid]-a[l-1];i++)ans[i+a[l-1]+1]=(ans[i+a[l-1]+1]+v[i])%mod;
79
80 v.clear();
81 for(int i=0;i<r-mid+a[mid]-a[l-1]-1;i++)v.push_back(fac[i]);
82 for(int i=0;i<a[mid]-a[l-1];i++)vu[i]=(ll)vu[i]*inv[a[mid]-a[l-1]-1-i]%mod;
83 v=mul(v,vu);
84 for(int i=0;i<r-mid;i++)f[i+mid+1]=(f[i+mid+1]+(ll)v[i+a[mid]-a[l-1]-1]*inv[i])%mod;
85 calc(mid+1,r);
86 }
87 int main(){
88 fac[0]=inv[0]=inv[1]=1;
89 for(int i=1;i<(N<<1);i++)fac[i]=(ll)fac[i-1]*i%mod;
90 for(int i=2;i<(N<<1);i++)inv[i]=(ll)(mod-mod/i)*inv[mod%i]%mod;
91 for(int i=1;i<(N<<1);i++)inv[i]=(ll)inv[i-1]*inv[i]%mod;
92 scanf("%d",&t);
93 while (t--){
94 scanf("%d",&n);
95 for(int i=1;i<=n;i++)scanf("%d",&a[i]);
96 n++,a[n]=n;
97 for(int i=1;i<=n;i++)f[i]=ans[i]=0;
98 f[1]=1;
99 calc(1,n);
100 printf("%d\n",ans[n]);
101 }
102 return 0;
103 }
[gym102220I]Temperature Survey的更多相关文章
- The 13th Chinese Northeast Collegiate Programming Contest
题解: solution Code: A. Apple Business #include<cstdio> #include<algorithm> #include<ve ...
- [LeetCode] Rising Temperature 上升温度
Given a Weather table, write a SQL query to find all dates' Ids with higher temperature compared to ...
- SharePoint 2010 Survey的Export to Spreadsheet功能怎么不见了?
背景信息: 最近用户报了一个问题,说他创建的Survey里将结果导出成Excel文件(Export to spreadsheet)的按钮不见了. 原因排查: 正常情况下,这个功能只存在于SharePo ...
- SharePoint Tricks - Survey
1. SharePoint 2010中,在Survey的问题框中输入HTML代码可以用于插入图片或者链接,具体方法为: 1.1 在问题框中输入html, 1.2 在New Form和Edit Form ...
- SharePoint - 添加图片到Survey的某一问题之上
Survey是SharePoint常用功能之一,而曾经被用户多次问到的问题是能否在Survey的某一问题上添加图片,经过查看,SharePoint Survey不提供此方法,只得谷歌之,得一比较懒但又 ...
- 转:关于数据库压缩技术的Survey
原文来自于:http://outofmemory.cn/mysql/database-compression-tech 昨天给团队内的小伙伴做了一个关于数据库压缩技术的Survey,现将其中可以公开的 ...
- BZOJ2276: [Poi2011]Temperature
2276: [Poi2011]Temperature Time Limit: 20 Sec Memory Limit: 32 MBSubmit: 293 Solved: 117[Submit][S ...
- leetcode-Rising Temperature
Given a Weather table, write a SQL query to find all dates' Ids with higher temperature compared to ...
- Application to find the maximum temperature in the weather dataset
import org.apache.hadoop.fs.Path; import org.apache.hadoop.io.IntWritable; import org.apache.hadoop. ...
随机推荐
- 20 行代码:Serverless 架构下用 Python 轻松搞定图像分类和预测
作者 | 江昱 前言 图像分类是人工智能领域的一个热门话题.通俗解释就是,根据各自在图像信息中所反映的不同特征,把不同类别的目标区分开来的图像处理方法. 它利用计算机对图像进行定量分析,把图像或图像中 ...
- 寻找最佳路径(ArcPy实现)
一.背景 随着社会经济发展需求,公路的重要性日益提高.在一些交通欠发达的地区,公路建设迫在眉睫.如何根据实际地形情况设计出比较合理的公路规划,是一个值得研究的问题. 二.实验目的: (1)通过练习,熟 ...
- 洛谷2093 JZPFAR + KD-Tree学习笔记 (KD-Tree)
KD-Tree这玩意还真的是有趣啊.... (基本完全不理解) 只能谈一点自己的对KD-Tree的了解了. 首先这个玩意就是个暴力... 他的结构有点类似二叉搜索树 每一层都是以一个维度作为划分标准. ...
- 纯前端H5小应用_localStorage存储
开发缘由[需求发现和分析] 想要送朋友一个礼物,但是想了想,街上买的东西,em~,我们这样的猿确实不会选礼物啊,由此就想利用自己手中的工具和知识做点有用的东西吧,抱枕是礼物,钢笔是礼物,电子产品也是礼 ...
- ❤️【Python从入门到精通】(二十七)更进一步的了解Pillow吧!
您好,我是码农飞哥,感谢您阅读本文,欢迎一键三连哦. 进一步介绍Pillow库的使用,详细了解 干货满满,建议收藏,需要用到时常看看. 小伙伴们如有问题及需要,欢迎踊跃留言哦~ ~ ~. 前言 本文是 ...
- 第十一章 Dockerfile安装Jenkins-2.249.3-1.1
一.安装Docker Docker部署Jenkins前提已经安装Docker,这边脚本安装Docker. #1.编写Docker安装脚本 [root@ip-10-0-12-212 ~]# vim In ...
- python OptionParser的用法
from optparse import OptionParser parser = OptionParser(usage = "usage: %prog [options] arg&quo ...
- 036—环境变量path
day04 课堂笔记 1.开发第一个java程序:HelloWorld 1.1.程序写完以后,一定要ctrl+s进行保存 源代码若修改,需重新进行编译 1.2.编译阶段 怎么编译?使用什么命令?这个命 ...
- 原生js-无缝滚动
<!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...
- matlab添加永久路径
addpath('D:\MATLAB6p5\toolbox\svm'); 临时添加路径,不能添加子目录 addpath(genpath('D:\MATLAB6p5\toolbox\svm'));临时添 ...