bzoj3992: [SDOI2015]序列统计

2024-04-02 10:48

本文主要是介绍bzoj3992: [SDOI2015]序列统计,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

传送门:http://www.lydsy.com:808/JudgeOnline/problem.php?id=3992

思路:M是一个质数,问题又是求乘积,于是我们就可以想到利用M的原根g把问题变成求和(我怎么想不到啊。。。)

根据原根的性质,我们可以把1到M-1中的数i表示为(g^b[i])%M,且指数互不相同

那么X就可以表示成(g^b[x])%M

问题就转化为:然后问题转化成了在序列b中,选出n个数(一个数可以取多次),且它们的和s满足s=b[x]

这个问题我们可以用母函数+NTT解决,答案就是多项式的b[x]次项的系数

因为n很大,所以再套一个快速幂即可

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int mod=(479<<21)+1,G=3,maxn=17000;
int n,k,sum,m,inv_G,inv_N,T[maxn],vis[maxn],tim,N,rev[maxn],pos[maxn],root;int qpow(int a,int b){int res=1;for (;b;b>>=1){if (b&1) res=1ll*res*a%mod;a=1ll*a*a%mod;}return res;
}bool check(int x){int now=1;++tim;for (int i=1;i<n;i++,now=now*x%n){if (vis[now]==tim) return 0;vis[now]=tim;}return 1;
}int rever(int x){int res=0,len=N;while (len--) res<<=1,res^=(x&1),x>>=1;return res;}
int findroot(){for (int i=2;i<=n;i++) if (check(i)) return i;}struct DFT{int a[maxn];void ntt(int op){for (int i=0;i<N;i++) if (rev[i]>i) swap(a[rev[i]],a[i]);int g=op==1?G:inv_G;for (int sz=2;sz<=N;sz<<=1){int t=qpow(g,(mod-1)/sz);for (int bg=0;bg<N;bg+=sz)for (int po=bg,w=1;po<bg+(sz>>1);po++){int x=a[po],y=1ll*a[po+(sz>>1)]*w%mod;a[po]=(x+y)%mod,a[po+(sz>>1)]=(x-y+mod)%mod;w=1ll*w*t%mod;}}if (op==-1) for (int i=0;i<N;i++) a[i]=1LL*a[i]*inv_N%mod;}
}a,b;void qpow(){b.a[0]=1;for (;k;k>>=1){a.ntt(1);if (k&1){b.ntt(1);for (int i=0;i<N;i++) b.a[i]=1ll*b.a[i]*a.a[i]%mod;b.ntt(-1);for (int i=N-1;i>=n-1;i--) b.a[i-n+1]=(b.a[i-n+1]+b.a[i])%mod,b.a[i]=0;//因为是在mod m(代码里的n)的条件下进行的,所以要把和超过m的后半部分的答案加到前半部分去,而不是简单的清空}for (int i=0;i<N;i++) a.a[i]=1ll*a.a[i]*a.a[i]%mod;a.ntt(-1);for (int i=N-1;i>=n-1;i--) a.a[i-n+1]=(a.a[i-n+1]+a.a[i])%mod,a.a[i]=0;}
}int main(){inv_G=qpow(G,mod-2);scanf("%d%d%d%d",&k,&n,&sum,&m);for (int i=1;i<=m;i++) scanf("%d",&T[i]);N=(int)ceil(log2(n))+1;for (int i=0;i<(1<<N);i++) rev[i]=rever(i);N=1<<N,inv_N=qpow(N,mod-2),root=findroot();for (int i=0,res=1;i<n-1;i++) pos[res]=i,res=res*root%n;for (int i=1;i<=m;i++) if (T[i]) a.a[pos[T[i]]]++;qpow(),printf("%d\n",b.a[pos[sum]]);return 0;
}


这篇关于bzoj3992: [SDOI2015]序列统计的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/869781

相关文章

C++从序列容器中删除元素的四种方法

《C++从序列容器中删除元素的四种方法》删除元素的方法在序列容器和关联容器之间是非常不同的,在序列容器中,vector和string是最常用的,但这里也会介绍deque和list以供全面了解,尽管在一... 目录一、简介二、移除给定位置的元素三、移除与某个值相等的元素3.1、序列容器vector、deque

一文详解SQL Server如何跟踪自动统计信息更新

《一文详解SQLServer如何跟踪自动统计信息更新》SQLServer数据库中,我们都清楚统计信息对于优化器来说非常重要,所以本文就来和大家简单聊一聊SQLServer如何跟踪自动统计信息更新吧... SQL Server数据库中,我们都清楚统计信息对于优化器来说非常重要。一般情况下,我们会开启"自动更新

最长公共子序列问题的深度分析与Java实现方式

《最长公共子序列问题的深度分析与Java实现方式》本文详细介绍了最长公共子序列(LCS)问题,包括其概念、暴力解法、动态规划解法,并提供了Java代码实现,暴力解法虽然简单,但在大数据处理中效率较低,... 目录最长公共子序列问题概述问题理解与示例分析暴力解法思路与示例代码动态规划解法DP 表的构建与意义动

关于最长递增子序列问题概述

《关于最长递增子序列问题概述》本文详细介绍了最长递增子序列问题的定义及两种优化解法:贪心+二分查找和动态规划+状态压缩,贪心+二分查找时间复杂度为O(nlogn),通过维护一个有序的“尾巴”数组来高效... 一、最长递增子序列问题概述1. 问题定义给定一个整数序列,例如 nums = [10, 9, 2

opencv实现像素统计的示例代码

《opencv实现像素统计的示例代码》本文介绍了OpenCV中统计图像像素信息的常用方法和函数,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一... 目录1. 统计像素值的基本信息2. 统计像素值的直方图3. 统计像素值的总和4. 统计非零像素的数量

如何使用 Bash 脚本中的time命令来统计命令执行时间(中英双语)

《如何使用Bash脚本中的time命令来统计命令执行时间(中英双语)》本文介绍了如何在Bash脚本中使用`time`命令来测量命令执行时间,包括`real`、`user`和`sys`三个时间指标,... 使用 Bash 脚本中的 time 命令来统计命令执行时间在日常的开发和运维过程中,性能监控和优化是不

hdu1496(用hash思想统计数目)

作为一个刚学hash的孩子,感觉这道题目很不错,灵活的运用的数组的下标。 解题步骤:如果用常规方法解,那么时间复杂度为O(n^4),肯定会超时,然后参考了网上的解题方法,将等式分成两个部分,a*x1^2+b*x2^2和c*x3^2+d*x4^2, 各自作为数组的下标,如果两部分相加为0,则满足等式; 代码如下: #include<iostream>#include<algorithm

uva 10131 最长子序列

题意: 给大象的体重和智商,求体重按从大到小,智商从高到低的最长子序列,并输出路径。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vect

flume系列之:查看flume系统日志、查看统计flume日志类型、查看flume日志

遍历指定目录下多个文件查找指定内容 服务器系统日志会记录flume相关日志 cat /var/log/messages |grep -i oom 查找系统日志中关于flume的指定日志 import osdef search_string_in_files(directory, search_string):count = 0

hdu4267区间统计

题意:给一些数,有两种操作,一种是在[a,b] 区间内,对(i - a)% k == 0 的加value,另一种操作是询问某个位置的值。 import java.io.BufferedInputStream;import java.io.BufferedReader;import java.io.IOException;import java.io.InputStream;import