NOIP2018模板总结【数学】】的更多相关文章

质因数分解 //质因数分解 int prime[MAXN], tim[MAXN], cnt; void Divide(int N) { printf("%d = ", N); for(int i = 2; i * i <= N; i++) if(N % i == 0) { prime[++cnt] = i; while(N % i == 0) N /= i, tim[cnt]++; } if(N > 1) prime[++cnt] = N, tim[cnt] = 1; pr…
D. Mahmoud and Ehab and another array construction task 因子分解模板 题意 给出一个原序列a 找出一个字典序大于a的序列b,使得任意 \(i!=j\),\(gcd(a[i],a[j])==1\),现在要你找出这样的序列b,并且满足所有合法序列中输出字典序最小的那个 思路 维护一个set,set里面装所有当前可以取的合法元素,先把所有的数字放进set里面,因为要求字典序最小的序列b,并且b的字典序要大于a,当构造的b到当前位置截止时和a相同时…
例子1 乘法除法运算: 1格式:{$number|number_format} 保留小数点后2位小数. {$v/$total*100|string_format:"%0.2f"|cat:'%'} 截取字符串 {$pcomment|truncate:176:"...":true}…
/*   smarty 引入对象 */ require_once 'libs/Smarty.class.php';  require 'MySmarty.class.php';  $msma = new MySmarty(); class stu{     public $name = '蒋金佚';     public $sex = '男';     function sayAge(){       echo '我今年21岁了<hr/>';     }  }  $stu1 = new stu…
暑假写的一些博客复习一遍.顺便再写一遍或者以现在的角度补充一点东西. 盛暑七月 初涉基环外向树dp&&bzoj1040: [ZJOI2008]骑士 比较经典的基环外向树dp.可以借鉴的技巧在于将每一个环拆出一条边,使剩下部分成为树.再然后就是max(f[u][0],f[v][0])思考中可能会出现的纰漏. ; i<=n; i++) { v[i] = read(), tt = read(); if (get(tt)!=get(i)){ addedge(i, tt); fa[fa[tt]…
手动博客搬家: 本文发表于20181125 13:25:03, 原地址https://blog.csdn.net/suncongbo/article/details/84487306 题目链接: https://www.luogu.org/problemnew/show/P4725 题目大意: 给定一个\(n\)次多项式\(A(x)\), 求一个\(n\)次多项式\(B(x)\)满足\(B(x)\equiv \ln A(x) (\mod x^n)\) 题解: 神数学模板题-- 数学真奇妙! 前驱…
%该模板用于数学答题 \documentclass[UTF8]{ctexart}%[中文编码 UTF8] \usepackage{fancyhdr}%{页眉页脚页码} \pagestyle{fancy}%{样式 fancy} \usepackage[fleqn]{amsmath}%[公式左对齐]{数学支持} \usepackage{amssymb}%{数学符号} \begin{document} \title{硕士研究生入学考试复试试卷答案}%标题 \author{DengChaohai}%作者…
hdu 4111  Alice and Bob 博弈:http://www.cnblogs.com/XDJjy/p/3350014.html hdu 4112 Break the Chocolate 数学题: http://www.cnblogs.com/XDJjy/p/3353386.html 努力半月10/1 ~10/17 结合2011上海交通大学命题的成都现场赛练习: 并: 基础数论题 简易的概率题 简单的几何题 简单的数据结构 易懂的建模搜索题 不是特复杂的模拟题 低中档的博弈 最后得整…
在使用文档写数学类的文章时候,常常会涉及到一些数学公式,由于数学公式中包含了很多的数学符号,如果使用文档自带的公式编辑器往往会发现很多的符号都不全或者不符合自己的要求.这个时候就需要一款专业的数学公式编辑器,MathType就是来解决这个问题的.利用MathType可以编辑出各种各样的数学符号与公式,很多常用函数也是可以编辑出来.下面来给大家介绍介绍用MathType编辑反三角函数的方法? 原文:http://www.mathtype.cn/jiqiao/fan-sanjiao.html 具体操…
第一段代码 \documentclass{article} \usepackage{ctex} \begin{document} \section{文字} 特可爱模板 \section{数学} \[ a^=b^+c^ \] \end{document} 第二段代码 \documentclass[UTF8]{ctexart} \title{杂谈勾股定理} \author{张三} \date{\today} \bibliographystyle{plain} \begin{document} \ma…
模板—数学—Exgcd Code: #include <cstdio> #include <algorithm> using namespace std; int ex_gcd(int a,int b,int &x,int &y) { if(!b) {x=1,y=0;return a;} int gcd=ex_gcd(b,a%b,x,y),tmp=x; x=y,y=tmp-a/b*y; return gcd; } int a,b,c,gcd,x1,y1; int m…
模板—数学—Lucas Code: #include <cstdio> #include <algorithm> using namespace std; #define N 100010 int n,m,p,inv[N],powq[N]; int lucas(int n,int m) { if(n<m) return 0; if(n<=p&&m<=p) return 1ll*powq[n]*inv[m]%p*inv[n-m]%p; return…
点此看题面 大致题意: 让你完成三种操作:求\(Y^Z\%P\)的值,求满足\(XY\equiv Z(mod\ P)\)的最小非负整数\(X\),求满足\(Y^X\equiv Z(mod\ P)\)的最小非负整数\(X\). 关于三个模板 只要你熟悉各类数学模板,就应该不难看出,其实这就是一道数学模板三合一. 第一个任务,显然是快速幂. 第二个任务,可是徐xgcd\(exgcd\)最经典的运用啊. 第三个任务,则是\(BSGS\)算法. 这样一来,这题就成了一道练模板的水题. 关于此题的数据 很…
首先让我们来谈一谈数学意义上的矩阵(在座各位也可以简单地将它理解为一个二维数组) 这样可以帮助我们理解矩阵加速及其运用的原理(矩阵加速是一个及其玄学的东西,所以请重点理解矩阵乘法)  这里给出一段严格的数学定义来帮助理解矩阵的概念,简单地看一下就可以了: 在数学中,矩阵(Matrix)是一个按照长方阵列排列的实数或复数集合,最早来自于方程组的系数及常数所构成的方阵.       矩阵是高等代数学中的常见工具,也常见于统计分析等应用数学学科中,矩阵的运算是数值分析领域的重要问题.在物理学中,矩阵在…
代码 请使用XeLatex编译 main.tex \documentclass{article} \usepackage{ctex} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 页面设置 \usepackage{geometry} %\geometry{a4paper,scale=0.8} %\usepackage{fontspec} %\setmonofont{Consolas} \usepa…
[SinGuLaRiTy-1047] Copyright (c) SinGuLaRiTy 2017. All Rights Reserved. 质因数分解 void solve(int n) { ==) { printf(); n/=; } ;i<=sqrt(n);i+=) { ) { n/=i; printf("%d*",i); i-=; } } printf("%d\n",n); } 欧拉线性筛素数 #define MAXN 100005 #define…
1.求逆元 int inv(int a) { ) ; return (MOD - MOD / a) * inv(MOD % a); } 2.线性筛法 bool isPrime[MAXN]; int label[MAXN], prime[MAXN]; int n, total; void makePrime() { n = ; ; i <= n; ++i) { if(!label[i]) { prime[total++] = i; label[i] = total; } ; j < label[…
一.素数 欧拉筛 void prime(){ check[]=; ;i<=n;i++){ if(!check[i])prim[++cnt]=i;//这个if语句后面没有大括号!! ;j<=cnt&&prim[j]*i<=n;j++){ check[i*prim[j]]=true; )break; } } } 简单的素数判定 bool check(int x){ )return false; ;i*i<=x;i++) )return false; return tru…
灵感来源是因为调试了一下午dij,就想把错误记下来 dij:结构体里的重载运算符不会写      结构体声明后要加引号     用于排序的结构体按照边长度排序 匈牙利:在dfs中的if语句中,要继续搜索的是fa[v]而不是v 倍增求lca:预处理出祖先时循环边界是dep[x]而不是n 割点:如果low[v]>=dfn[u]并且u!=fa,u才是一个割点,而不是low[v]>=low[u] 线段树求k小值:k_th_lazy[]只有在有值的时候才会被更新/下放(建议另开一个函数记录,这样既简洁又…
普通的扩展欧几里得算法,通过了洛谷的扩展欧几里得算法找乘法逆元.修复了容易溢出的bug,虽然新版本仍有可能会溢出longlong,假如参与运算的数字都是longlong,假如可以的话直接使用__int128或者去抄一个RoundGod的BigInt的模板(这里的C题).事不宜迟明天就抄这个大数模板. 求解模n意义下a的逆元,即求方程LCE2(a,1,n,x),结果放入x中,返回值指示是否有解. ll gcd(ll a, ll b) { if(b == 0) return a; while(ll…
ACM常用模板合集 #include <bits/stdc++.h> using namespace std; typedef long long ll; ll pr; ll pmod(ll a, ll b, ll p) { return (a * b - (ll)((long double)a / p * b) * p + p) % p; } //普通的快速乘会T ll gmod(ll a, ll b, ll p) { ll res = 1; while (b) { if (b &…
筛选法+求一个整数的分解+快速模幂运算+递归求计算1+p+p^2+````+p^nPOJ 1845 Sumdiv求A^B的所有约数之和%9901 */#include<stdio.h>#include<math.h>#include<iostream>#include<algorithm>#include<string.h>using namespace std;#define MOD 9901const int MAXN=10000;int p…
传送门 Description 在网友的国度中共有 \(n\) 种不同面额的货币,第 \(i\) 种货币的面额为 \(a[i]\),你可以假设每一种货币都有无穷多张.为了方便,我们把货币种数为 \(n\).面额数组为 \(a[1..n]\) 的货币系统记作 \((n,a)\). 在一个完善的货币系统中,每一个非负整数的金额 \(x\) 都应该可以被表示出,即对每一个非负整数 \(x\),都存在 \(n\) 个非负整数 \(t[i]\) 满足 \(a[i] \times t[i]\) 的和为 \(…
题目描述 给出一个有理数 c=a/b ,求 c mod 19260817的值. 说明 对于所有数据, 0≤a,b≤10^10001 分析: 一看题 这么短 哇简单!况且19260817还是个素数!(美滋滋 再一看数据 我天 可怕 10^10001 一看 完了 要打高精 (但我打高精肯定GG啊 一想 根据同余好像可做(前面的大佬讲过了我就不赘述了 哦 对了 还有费马小定理: a^phi(p)≡1(mod p) (只对于p是质数的情况哦 然后对于快读 略做修改就可以了 ---------------…
$O(n)$递推求逆元 #include<cstdio> #include<cstring> #include<algorithm> using namespace std; typedef long long ll; ]; int main(){ int n,p; scanf("%d%d",&n,&p); inv[]=; printf("1\n"); ;i<=n;i++){ inv[i]=(ll)(p-p/…
#1068 : RMQ-ST算法 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 小Hi和小Ho在美国旅行了相当长的一段时间之后,终于准备要回国啦!而在回国之前,他们准备去超市采购一些当地特产——比如汉堡(大雾)之类的回国. 但等到了超市之后,小Hi和小Ho发现者超市拥有的商品种类实在太多了——他们实在看不过来了!于是小Hi决定向小Ho委派一个任务:假设整个货架上从左到右拜访了N种商品,并且依次标号为1到N,每次小Hi都给出一段区间[L, R],小Ho要做的是选出标…
先看看. 通常模数常见的有998244353,1004535809,469762049,这几个的原根都是3.所求的项数还不能超过2的23次方(因为998244353的分解). 感觉没啥用. #include <cstdio> #include <cstring> template <class T> inline void swap(T &a, T &b) { T c; c = a; a = b; b = c; } ; , G = ; inline in…
费马(Fermat)小定理 当 \(p\) 为质数,则 \(a^{p-1}\equiv 1 \mod p\) 反之,费马小定理的逆定理不成立,这样的数叫做伪质数,最小的伪质数是341. 欧拉(Euler)定理 扩展欧拉(Euler)定理 根据扩展欧拉定理,不管a和p是不是互质,都可以缩小到 \([\varphi(p),2\varphi(p)]\) 之间,然后暴力用快速幂求解.…
使用Fermat小定理(Fermat's little theorem)的原理进行测试,不满足 \(2^{n-1}\;\mod\;n\;=\;1\) 的n一定不是质数:如果满足的话则多半是质数,满足上式(通过2为底的Fermat小定理测试)且是合数的,被称为"伪质数"(pseudoprime number),一个简单的伪质数是341.一个合数可能在a=2时通过了测试,但a=3时的计算结果却排除了素数的可能.于是,人们扩展了伪素数的定义,称满足 \(a^{n-1}\;\mod\;n\;=…
终于知道发明者的正确的名字了,是Min_25,这个筛法速度为亚线性的\(O(\frac{n^{\frac{3}{4}}}{\log x})\),用于求解具有下面性质的积性函数的前缀和: 在 \(p\) 处是简单的低次多项式 在 \(p^c\) 处可以快速求值 貌似积性函数是指取一个积性函数 \(f(x)\) ,其在质数的位置上取值与所求函数相同.所以可以用来求n以内的质数的个数(取常函数 \(f(x)=1\) )以及质数的和(取恒等函数 \(f(x)=x\) ). 参考资料: loj#6235.…