
ACM-数学
「已注销」
这个作者很懒,什么都没留下…
展开
-
数学问题合集
求n!首位数使用斯特林公式得到近似结果: limn→+∞n!2πn−−−√(ne)n=1 \lim_{n\to+\infty} {n!\over{\sqrt {2 \pi n}} ({n\over e})^n} = 1之后用10的幂计算取首位即可,代码如下:#include <stdio.h>#include <math.h>const double PI =3.14159265358979;原创 2017-11-08 11:59:16 · 302 阅读 · 0 评论 -
欧拉筛法求素数
C代码#include <stdio.h>#include <stdbool.h>bool flag[10001] = {0};int p[10001] = {0};int main(){ int i, n, cnt = 0, j; scanf("%d", &n); for (i = 2; i <= n; ++i...原创 2017-12-04 21:37:04 · 359 阅读 · 0 评论 -
EOJ 1748 Box Walking
题目简介给定长方体的长宽高以及长方体表面的一个点(x,y,z),求从点(0,0,0)出发到(x,y,z)的最短距离的平方。说明这是个看似简单实际上也简单的问题。首先如果该点在(0,0,0)所在的三个平面之一之内,那么显然连线长度就是答案。对于其他情况,考虑平铺长方体的六个面,这时连线的长度也几乎是答案,但是画画图不难发现有时候连线会跑到平铺后的长方体外面去……这意味着我们需...原创 2018-02-23 21:09:40 · 180 阅读 · 0 评论 -
EOJ 2018.1.29新生训练Week2
总体来说不算好做……但是确实也不难。A给定正整数n,k,求f(n,k)=∑i=1nik" role="presentation" style="position: relative;">f(n,k)=∑ni=1ikf(n,k)=∑i=1nikf(n,k)=\sum^n_{i=1}i^k,结果对19260817取模。 1 时限0.698s(???)上来用拉格朗日原创 2018-02-05 15:21:35 · 678 阅读 · 0 评论 -
EOJ 2018.2.5新生训练Week3 D.Game of Chairs
n个椅子,c种颜色排成一圈,间隔1m。随机选一种颜色,你要马上移动到这种颜色的椅子上(原本颜色相同则不动)。求走动距离的最小期望(输出最简分数)。 1 ≤ c ≤ n ≤ 1e6据说暴力模拟+优化(n^2—>nlogn?)2.5s内可过……不过这里用了一些数学知识,复杂度降到O(n)。(其实是2n,不过eoj评测姬太快了,可以忽略)首先肯定是常规的环拆链操作:复制一...原创 2018-02-10 17:55:17 · 353 阅读 · 0 评论 -
高斯消元法求行列式及逆矩阵
题目简介给定一个 n*n 的矩阵,输出它的行列式值和逆矩阵(保证存在)。说明期末复习线性代数时,发现自己求逆矩阵总是求错,于是干脆写了个程序来实现。。至于是列主元还是全主元还是约当什么的……以后来填这个坑。 高斯消元法求解线性方程组只要稍微修改下代码就可以,判断无解/无穷多解也不难,至于求自由未知量然后输出任意一解……还没有想好,也是以后再来填(逃#include us原创 2018-01-28 12:57:02 · 5492 阅读 · 0 评论 -
轻量级的巧妙解法
EOJ 1076 染气球 一共有 N 只气球,小强将 N 只气球从左到右依次编号为 1、2、3……N,每次给 2 个整数 a,b (a<=b),小强便骑上他的“小飞鸽 ” 牌电动车从气球 a 开始到气球 b 依次给每个气球涂一次颜色。但是 N 次以后小强已经忘记了第 I 个气球已经涂过几次颜色了,你能帮他算出每个气球被涂过几次颜色吗? N <= 100000da...原创 2018-02-04 12:52:43 · 304 阅读 · 0 评论 -
Miller-Rabin素性测试与二次探测
算法简介首先是一些概念:费马小定理:对于素数p和任意整数a,有ap≡a(modp)a^p ≡ a(mod p).反之,对于一个数p,如果满足ap≡a(modp)a^p ≡ a(mod p),则 p 很可能是素数。伪素性测试:瞎猜若干个x,只要不满足上式,那么p就不是素数。看起来没毛病了。Carmichael数:对于合数n,如果对所有正整数b(b和n互素)都有bn−1≡1(mod原创 2018-01-19 16:34:53 · 1583 阅读 · 2 评论 -
EOJ 1277/Ulm Local 1996 Humble Numbers
题目简介A number whose only prime factors are 2,3,5 or 7 is called a humble number. The sequence 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 15, 16, 18, 20, 21, 24, 25, 27, … shows the first 20 humble numb原创 2018-01-19 16:46:49 · 215 阅读 · 0 评论 -
EOJ 1499 矩阵快速幂求斐波那契数列
题目简介给定n,求斐波那契数列前n项和。 0 < n < 1e9说明快速幂实现方式和整数差不多,没什么好讲的。而对于斐波那契数列,不难发现: (Fn+2Fn+1)=(1110)(Fn+1Fn)(1)(1)(Fn+2Fn+1)=(1110)(Fn+1Fn)\begin{equation}\begin{pmatrix} F_{n+2}\\F_{n+1}\\...原创 2018-01-10 20:11:58 · 232 阅读 · 0 评论 -
计算多项式的系数
题目简介给定一个多项式 (ax+by)k" role="presentation" style="position: relative;">(ax+by)k(ax+by)k(ax+by)^k,计算多项式展开后 xnym" role="presentation" style="position: relative;">xnymxnymx^ny^m 项的系数。Input 第 1 行:一个原创 2018-01-25 15:32:00 · 9318 阅读 · 0 评论 -
EOJ 3322 万年历(历史题)
为什么说是历史题?来自某百科: 罗马教皇格里高利十三世在1582年组织了一批天文学家,根据哥白尼日心说计算出来的数据,对儒略历作了修改。将1582年10月5日到14日之间的10天宣布撤销,继10月4日之后为10月15日。后来人们将这一新的历法称为“格里高利历”,也就是今天世界上所通用的历法,简称格里历或公历。所以怎么算星期几?蔡勒公式(Zeller formul原创 2017-11-08 12:05:54 · 413 阅读 · 0 评论 -
EOJ 2784/Project Euler 120改 remainder
题目简介Let r be the remainder when (a-1)n + (a+1)n is divided by a2 .For example, if a = 7 and n = 3, then r = 42: ( 63 + 83 = 728 )%49= 42. And as n varies, so too will r, but for a = 7 it turns o原创 2018-01-04 13:54:23 · 220 阅读 · 0 评论 -
EOJ 1159/POJ 3286 How many 0's?
http://acm.ecnu.edu.cn/problem/1159/ 给定m,n问从m到n之间有多少0。又到了令人头大的数位dp时间…… 不过这里没有用数位dp,直接前缀和+递推解决的。比较费草稿纸而不是键盘,属于思维题吧。#include <bits/stdc++.h>using namespace std;typedef long long ll;...原创 2018-03-04 20:50:33 · 224 阅读 · 0 评论