(为了方便,以下记$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的更多相关文章

  1. The 13th Chinese Northeast Collegiate Programming Contest

    题解: solution Code: A. Apple Business #include<cstdio> #include<algorithm> #include<ve ...

  2. [LeetCode] Rising Temperature 上升温度

    Given a Weather table, write a SQL query to find all dates' Ids with higher temperature compared to ...

  3. SharePoint 2010 Survey的Export to Spreadsheet功能怎么不见了?

    背景信息: 最近用户报了一个问题,说他创建的Survey里将结果导出成Excel文件(Export to spreadsheet)的按钮不见了. 原因排查: 正常情况下,这个功能只存在于SharePo ...

  4. SharePoint Tricks - Survey

    1. SharePoint 2010中,在Survey的问题框中输入HTML代码可以用于插入图片或者链接,具体方法为: 1.1 在问题框中输入html, 1.2 在New Form和Edit Form ...

  5. SharePoint - 添加图片到Survey的某一问题之上

    Survey是SharePoint常用功能之一,而曾经被用户多次问到的问题是能否在Survey的某一问题上添加图片,经过查看,SharePoint Survey不提供此方法,只得谷歌之,得一比较懒但又 ...

  6. 转:关于数据库压缩技术的Survey

    原文来自于:http://outofmemory.cn/mysql/database-compression-tech 昨天给团队内的小伙伴做了一个关于数据库压缩技术的Survey,现将其中可以公开的 ...

  7. BZOJ2276: [Poi2011]Temperature

    2276: [Poi2011]Temperature Time Limit: 20 Sec  Memory Limit: 32 MBSubmit: 293  Solved: 117[Submit][S ...

  8. leetcode-Rising Temperature

    Given a Weather table, write a SQL query to find all dates' Ids with higher temperature compared to ...

  9. 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. ...

随机推荐

  1. 20 行代码:Serverless 架构下用 Python 轻松搞定图像分类和预测

    作者 | 江昱 前言 图像分类是人工智能领域的一个热门话题.通俗解释就是,根据各自在图像信息中所反映的不同特征,把不同类别的目标区分开来的图像处理方法. 它利用计算机对图像进行定量分析,把图像或图像中 ...

  2. 寻找最佳路径(ArcPy实现)

    一.背景 随着社会经济发展需求,公路的重要性日益提高.在一些交通欠发达的地区,公路建设迫在眉睫.如何根据实际地形情况设计出比较合理的公路规划,是一个值得研究的问题. 二.实验目的: (1)通过练习,熟 ...

  3. 洛谷2093 JZPFAR + KD-Tree学习笔记 (KD-Tree)

    KD-Tree这玩意还真的是有趣啊.... (基本完全不理解) 只能谈一点自己的对KD-Tree的了解了. 首先这个玩意就是个暴力... 他的结构有点类似二叉搜索树 每一层都是以一个维度作为划分标准. ...

  4. 纯前端H5小应用_localStorage存储

    开发缘由[需求发现和分析] 想要送朋友一个礼物,但是想了想,街上买的东西,em~,我们这样的猿确实不会选礼物啊,由此就想利用自己手中的工具和知识做点有用的东西吧,抱枕是礼物,钢笔是礼物,电子产品也是礼 ...

  5. ❤️【Python从入门到精通】(二十七)更进一步的了解Pillow吧!

    您好,我是码农飞哥,感谢您阅读本文,欢迎一键三连哦. 进一步介绍Pillow库的使用,详细了解 干货满满,建议收藏,需要用到时常看看. 小伙伴们如有问题及需要,欢迎踊跃留言哦~ ~ ~. 前言 本文是 ...

  6. 第十一章 Dockerfile安装Jenkins-2.249.3-1.1

    一.安装Docker Docker部署Jenkins前提已经安装Docker,这边脚本安装Docker. #1.编写Docker安装脚本 [root@ip-10-0-12-212 ~]# vim In ...

  7. python OptionParser的用法

    from optparse import OptionParser parser = OptionParser(usage = "usage: %prog [options] arg&quo ...

  8. 036—环境变量path

    day04 课堂笔记 1.开发第一个java程序:HelloWorld 1.1.程序写完以后,一定要ctrl+s进行保存 源代码若修改,需重新进行编译 1.2.编译阶段 怎么编译?使用什么命令?这个命 ...

  9. 原生js-无缝滚动

    <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...

  10. matlab添加永久路径

    addpath('D:\MATLAB6p5\toolbox\svm'); 临时添加路径,不能添加子目录 addpath(genpath('D:\MATLAB6p5\toolbox\svm'));临时添 ...