摘要: 题意是要我们求$\tan (x)$,可以用$\tan (x)=\frac{\sin (x)}{\cos (x)}$,或者这个 $\tan (x)=\sum _{x=1}^{\infty } \frac{\left(( 1)^{n 1} 2^{2 n} \left(2^{2 n} 1\right) B 阅读全文
posted @ 2017-07-17 16:43 czp001 阅读(206) 评论(0) 推荐(0) 编辑
摘要: 二次无理数的连分数是循环的,循环节从$a_1$开始然后到某一个$a_i=2 a_1$的时候结束,后面以此循环,即$\left[a_0;a_1,a_2,a_3,\text{...},2 a_1\right]$,[wiki百科][xx]里有一个针对$\sqrt{n}$的连分数求法,不涉及浮点数求倒数。 阅读全文
posted @ 2017-02-22 01:03 czp001 阅读(213) 评论(0) 推荐(0) 编辑
摘要: 模板题,多项式求逆元,FFT/NTT啦。想看[FFT/NTT系列教程][1]请点击,感觉写得很好。 cpp include using namespace std; define LL long long define N 530000 define M 530000 //const LL P = 阅读全文
posted @ 2017-02-18 22:41 czp001 阅读(163) 评论(0) 推荐(0) 编辑
摘要: 这题分块搞一搞,算到最后发现结果就是算矩阵J的行列式,要取模m,那个sign消掉了。 参考文献在这,[欧几里德算法的应用,金斌][1] 代码还是比较好写的,python代码 [1]: https://github.com/ldfaiztt/algorithm 3/blob/master/resour 阅读全文
posted @ 2017-02-04 00:52 czp001 阅读(307) 评论(0) 推荐(0) 编辑
摘要: 公式推导偷个懒,直接用别人的了 如果117是p的二次剩余,矩阵乘法部分,M的循环节就是p-1,否则是p+1。判断方法是勒让德符号,很好写,快速幂判断117的(p-1)/2次方是不是1.当然你也可以不判循环节,循环节取p^2-1,但是会慢。 python写的,本来做好了用c++重写的准备的,没想到居然 阅读全文
posted @ 2017-02-02 18:29 czp001 阅读(317) 评论(0) 推荐(0) 编辑
摘要: 这两题原理是一样的,不过第二题数据量大一些。这个累加操作相当于一个矩阵乘法,然后用矩阵中的第一列数和输入的数组做卷积,比如这个样例处理2次矩阵就是这样的 然后取出第一列数和输入的数做卷积,也就是多项式乘法 $\left(1+2x+3x^2+4x^3\right) \left(1+3x+5x^2+6x 阅读全文
posted @ 2017-01-22 14:05 czp001 阅读(198) 评论(0) 推荐(0) 编辑
摘要: 这题大概有两种解法: 1.预先打表放到数组里,大概100万间隔的素数的值$p_n$,然后用区间筛法。可以分的再细一点。 2.先估计,再用区间筛法。你需要一个第n个素数的估计值,一个素数统计函数和区间筛法,根据wiki上的关于第n个素数的公式,通过打表,我发现他有时大于实际值,有时小于实际值,我觉得让 阅读全文
posted @ 2017-01-21 23:46 czp001 阅读(829) 评论(0) 推荐(0) 编辑
摘要: 这个博客用来记录一些算法相关的东西吧,可能还会有别的什么的。 代码主要是C/C++和Python 希望可以坚持下去。 哈哈哈!! 好像立了flag。。。 就这样开始吧!! 炮姐镇一下。。。 阅读全文
posted @ 2017-01-18 22:51 czp001 阅读(113) 评论(0) 推荐(0) 编辑