- 博客(61)
- 资源 (2)
- 收藏
- 关注
原创 1951: [Sdoi2010]古代猪文
题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=1951题目大意:计算G^∑i|nC(n,i)%P题解:数论+中国剩余定理+费马小定理+逆元+lucas定理。根据 费马小定理a^(p-1)≡1(mod p)设∑i|nC(n,i)=M;设M=k*(p-1)+W:则G^∑i|nC(n,i)=G^M=G^(k*(p-1)...
2019-09-05 21:02:37
254
原创 1121: [POI2008]激光发射器SZK
题目:多边形相邻边垂直,边长为整数,边平行坐标轴。要在多边形的点上放一些激光发射器和接收器。满足下列要求: 1发射器和接收器不能放置在同一点; 2发射器发出激光可以沿壁反射,最终到达一个接收器; 3发射器只能沿角平分线发射激光。求:最多可放置多少对发射器和接收器?点数4<=n<=100000题解:因为根据光路的可逆性,从A射到B,必然有从B射到A,所以不可能有两个点同时射到同一个点...
2019-09-05 19:46:45
275
原创 1008: [HNOI2008]越狱
题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=1008题目大意:监狱有连续编号为1...N的N个房间,每个房间关押一个犯人,有M种宗教,每个犯人可能信仰其中一种。如果相邻房间的犯人的宗教相同,就可能发生越狱,求有多少种状态可能发生越狱思路:用所有的情况减去不可能越狱的情况。#include <bits/stdc++...
2019-09-05 19:33:55
251
原创 1041: [HAOI2008]圆上的整点
题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=1041题目大意:求一个给定的圆(x^2+y^2=r^2),在圆周上有多少个点的坐标是整数。题解:∵ x^2+y^2=r^2∴ x^2=r^2-y^2∴x^2=(r-y)(r+y) 设d=gcd(r-y,r+y); r-y=d*a^2; r+y=d*b^2; (a&l...
2019-09-04 21:05:33
272
原创 bzoj 3283: 运算器
题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=3283快速幂+拓展lucas+拓展GSBS#include <bits/stdc++.h>using namespace std;typedef long long ll;void exgcd( ll a, ll b, ll &d, ll &...
2019-09-01 19:47:27
247
原创 快速卢卡斯模板
ll fac[maxn];void init(){ fac[0] = 1; for(int i = 1; i <= maxn; i++) fac[i] = fac[i-1] * i % MOD; return;}ll _pow(ll x, ll y){ ll res = 1,tmp = x % MOD; while(y) ...
2019-08-16 16:26:03
223
原创 积性函数与线性筛
常见积性函数:μ(n):莫比乌斯函数φ(n):欧拉函数d(n):一个数nn的约数个数σ(n):一个数nn的约数和f(x)=xk(k∈N):这个玩意儿也是积性函数线性筛素数:int pri[N],tot,isprime[N];//isprime[i]为1的表示不是质数void sieve(){ isprime[1]=1; for (int i=2;i<...
2019-05-15 15:37:33
234
原创 HDU - TrickGCD
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6053题目大意:给一个长度为n的数组,构造出一个长度为n的数组满足1≤Bi≤Ai,且任意区间gcd>=2;思路:对于每一个gcd都有a[i]/gcd个然后把每个值乘起来,对于所有序列的重复性恰好为莫比乌斯函数的规律。(莫比乌斯函数+容斥+筛法)#include <bits/st...
2019-05-12 18:19:00
229
原创 第十届蓝桥杯题解
试题 A: 平方和本题总分:5 分【问题描述】小明对数位中含有 2、0、1、9 的数字很感兴趣,在 1 到 40 中这样的数包括 1、2、9、10 至 32、39 和 40,共 28 个,他们的和是 574,平方和是 14362。注意,平方和是指将每个数分别平方后求和。请问,在 1 到 2019 中,所有这样的数的平方和是多少?...
2019-03-29 17:32:35
742
原创 Edgy Trees CodeForces - 1139C
题目链接:http://codeforces.com/problemset/problem/1139/C题目大意:给了一棵树,n个点,m条边。让从中选k个点,使得从a1到a2,a2到a3,ak-1到ak的路径中至少经过一条黑色的边,问这样的集合有多少个。思路:并查集求红色边的所用联通块,然后用总的情况减去联通块的集合(用快速幂求)。代码:#include <bits/std...
2019-03-22 17:09:31
250
原创 Farey Sequence
题目链接:https://vjudge.net/problem/POJ-2478The Farey Sequence Fn for any integer n with n >= 2 is the set of irreducible rational numbers a/b with 0 < a < b <= n and gcd(a,b) = 1 arranged i...
2018-11-26 19:12:40
452
原创 2017ccpc杭州(数论)
题目链接:https://vjudge.net/contest/268674#problem/B题意:依次输入一个整数分解质因数的各项底数与指数,求算: ∑d|nφ(d)×n/d思路:迪利克雷卷积知识和欧拉函数的知识可以看出它是个积性函数推导如下图:代码:#include <bits/stdc++.h>using namespace std;...
2018-11-26 17:29:29
403
原创 “东信杯”广西大学第一届程序设计竞赛(同步赛)
链接:https://ac.nowcoder.com/acm/contest/283/B来源:牛客网 不吉利的数时间限制:C/C++ 2秒,其他语言4秒空间限制:C/C++ 131072K,其他语言262144K64bit IO Format: %lld题目描述在数学中,某个序列的子序列是从最初序列通过去除某些元素但不破坏余下元素的相对位置(在前或在后)而形成的新序列。...
2018-11-25 18:54:23
1939
原创 Pairs Forming LCM LightOJ - 1236
题目链接:https://cn.vjudge.net/problem/LightOJ-1236Find the result of the following code:long long pairsFormLCM( int n ) { long long res = 0; for( int i = 1; i <= n; i++ ) for( int ...
2018-11-20 20:42:55
354
原创 Help Hanzo LightOJ - 1197
题目链接:https://cn.vjudge.net/problem/LightOJ-1197Amakusa, the evil spiritual leader has captured the beautiful princess Nakururu. The reason behind this is he had a little problem with Hanzo Hattori, ...
2018-11-20 20:27:24
280
原创 Trailing Zeroes (III) LightOJ - 1138
题目链接:https://cn.vjudge.net/problem/LightOJ-1138You task is to find minimal natural number N, so that N! contains exactly Q zeroes on the trail in decimal notation. As you know N! = 1*2*...*N. For ex...
2018-11-20 20:17:11
421
原创 GCD - Extreme (II) UVA - 11426
题目链接:https://cn.vjudge.net/problem/UVA-11426题意:跟你一个n让你求(1,2....n-1)分别到n的所有gcd之和。思路:欧拉函数的运用。我们用b[n]表示(1....n-1)与ngcd之和。得到递推式sum[n]=sum[n-1]+b[n];设a[i]表示gcd(n,(1到n-1))=i;我们知道gcd(n/i,(1到n-1)/i)=1;我...
2018-11-20 20:13:53
244
原创 三角形
链接:https://ac.nowcoder.com/acm/contest/75/I来源:牛客网 题目描述 给你一个三角形的顶点A,B,C的坐标(坐标都为整数),请求出三角形的面积,三角形内的点的个数以及边AB、BC和AC边上的点的个数(不包括顶点ABC)输入描述:多组输入每组输入三行,每行两个整数第一行顶点A的坐标Xa,Ya.第二行顶点B的坐标Xb,Yb.第三...
2018-11-15 20:19:56
328
原创 小咪买东西
链接:https://ac.nowcoder.com/acm/contest/52/E来源:牛客网 题目描述小咪是一个土豪手办狂魔,这次他去了一家店,发现了好多好多(n个)手办,但他是一个很怪的人,每次只想买k个手办,而且他要让他花的每一分钱都物超所值,即:买下来的东西的总价值/总花费=max。请你来看看,他会买哪些东西吧。输入描述: 多组数据。第一行一个整数T,为...
2018-11-08 20:02:47
303
原创 Harmonic Number(LightOJ - 1234)
In mathematics, the nth harmonic number is the sum of the reciprocals of the first n natural numbers:In this problem, you are given n, you have to find Hn.InputInput starts with an integer T...
2018-11-08 19:31:44
212
原创 Leading and Trailing(LightOJ - 1282)
You are given two integers: n and k, your task is to find the most significant three digits, and least significant three digits of nk.InputInput starts with an integer T (≤ 1000), denoting the num...
2018-11-08 19:28:38
195
原创 Aladdin and the Flying Carpet(LightOJ - 1341)
It's said that Aladdin had to solve seven mysteries before getting the Magical Lamp which summons a powerful Genie. Here we are concerned about the first mystery.Aladdin was about to enter to a magi...
2018-11-08 19:18:34
217
原创 牛客网-数的变换
链接:https://ac.nowcoder.com/acm/contest/52/C题目描述给定一个序列数列, (ai互不相同)3种映射关系现在对于给定询问(x, k),求输入描述:多组读入数据T,对于每组数据,第1行一个整数n,q,n表示数列的大小,q为询问数第2行读入n个数a1,a2,…,an,表示数列中的数接下来q行,每行读入2个整...
2018-11-08 19:07:47
472
原创 快速排序
介绍:快速排序从它名字可以看出它的特点就是速度快、效率高,由于排序效率在同为O(N*logN)的几种排序方法中效率较高的一种。方法: 快排的工作过程其实比较简单,三步走: 选择基准值 pivot 将数组分成两个子数组:小于基准值的元素和大于基准值的元素。这个过程称之为 partition 对这两个子数组进行快速排序。 合并结果。 过程:...
2018-10-31 17:19:46
182
原创 ACM18级招新赛Round#1题解
A#include <bits/stdc++.h>using namespace std;int main(){ int t; cin>>t; while(t--) { long long a,b; cin>>a>>b; cout<<b<<endl; } }B#includ...
2018-10-20 19:50:36
942
原创 数论之神 牛客国庆集训派对Day5
链接:https://www.nowcoder.com/acm/contest/205/L来源:牛客网 题目描述终于活成了自己讨厌的样子。这是她们都还没长大的时候发生的故事。那个时候,栗子米也不需要为了所谓的爱情苦恼。她们可以在夏日的午后,花大把的时间去研究生活中一些琐碎而有趣的事情,比如数论。有一天西柚柚问了栗子米一个题,她想知道中有多少不同的数,这些不同的数字里面第k大...
2018-10-05 14:20:18
747
原创 Wannafly挑战赛19 -A队列Q
链接:https://www.nowcoder.com/acm/contest/131/A来源:牛客网 题目描述ZZT 创造了一个队列 Q。这个队列包含了 N 个元素,队列中的第 i 个元素用 Qi 表示。Q1 表示队头元素,QN 表示队尾元素。队列中的元素是 N 的一个全排列。ZZT 需要在这个队列上执行 P 次操作,操作分两种:FIRST X: 将元素 X 移到队头。LAST...
2018-07-23 10:28:09
302
原创 Manacher(马拉车)算法总结
Manacher算法 在介绍算法之前,首先介绍一下什么是回文串,所谓回文串,简单来说就是正着读和反着读都是一样的字符串,比如abba,noon等等,一个字符串的最长回文子串即为这个字符串的子串中,是回文串的最长的那个。 计算字符串的最长回文字串最简单的算法就是枚举该字符串的每一个子串,并且判断这个子串是否为回文串,这个算法的时间复杂度为O(n^3)的,显然无法令人满意,稍微优化的一个算...
2018-07-21 17:31:32
3867
2
原创 一个人的旅行 HDU - 2066
http://acm.hdu.edu.cn/showproblem.php?pid=2066虽然草儿是个路痴(就是在杭电待了一年多,居然还会在校园里迷路的人,汗~),但是草儿仍然很喜欢旅行,因为在旅途中 会遇见很多人(白马王子,^0^),很多事,还能丰富自己的阅历,还可以看美丽的风景……草儿想去很多地方,她想要去东京铁塔看夜景,去威尼斯看电影,去阳明山上看海芋,去纽约纯粹看雪景,去巴黎喝咖啡写...
2018-07-20 10:50:25
236
原创 飞行路线 HYSBZ - 2763
Alice和Bob现在要乘飞机旅行,他们选择了一家相对便宜的航空公司。该航空公司一共在n个城市设有业务,设这些城市分别标记为0到n-1,一共有m种航线,每种航线连接两个城市,并且航线有一定的价格。Alice和Bob现在要从一个城市沿着航线到达另一个城市,途中可以进行转机。航空公司对他们这次旅行也推出优惠,他们可以免费在最多k种航线上搭乘飞机。那么Alice和Bob这次出行最少花费多少? 数据...
2018-07-19 16:55:08
1085
1
原创 扩展欧几里德
定理一:如果d = gcd(a, b),则必能找到正的或负的整数k和l,使d = a*x+ b*y。定理二:若gcd(a, b) = 1,则方程ax ≡ c (mod b)在[0, b-1]上有唯一解。定理三:若gcd(a, b) = d,则方程ax ≡ c (mod b)在[0, b/d - 1]上有唯一解。证明:上述同余方程等价于ax + by = c,如果有解,两边同除以d,就有a/d...
2018-07-18 10:48:29
255
原创 Wannafly挑战赛19-A
链接:https://www.nowcoder.com/acm/contest/131/A来源:牛客网题目描述ZZT 创造了一个队列 Q。这个队列包含了 N 个元素,队列中的第 i 个元素用 Qi 表示。Q1 表示队头元素,QN 表示队尾元素。队列中的元素是 N 的一个全排列。 ZZT 需要在这个队列上执行 P 次操作,操作分两种: FIRST X: 将元素 X 移到队头。 LAST X: 将...
2018-07-11 10:42:01
241
原创 拓扑排序
#include<iostream>#include<string.h>#include<cstdio>using namespace std;const int maxn = 1000 + 5;int a[maxn][maxn];int c[maxn];int topo[maxn];int n, m, t,ans=0;bool dfs(i...
2018-05-09 21:38:15
188
原创 北京师范大学第十六届程序设计竞赛决赛-重现赛 G题
驼峰命名法是起变量名的一种规范,大致来说是用混合的大小写字母来构成变量名,在这个问题里你可以假设具体规则如下: 1.每个变量名由至少2个单词拼接构成,且每个单词长度至少为2; 2.每个单词的首字母必须大写,其他位置必须小写(除了变量名的第一个单词允许全部小写外)。 但是SK同学的英语很差,看到长长的变量名就很难脑补出是由哪些单词组成的,因此看驼峰...
2018-05-02 19:20:12
341
原创 poj-1019
A single positive integer i is given. Write a program to find the digit located in the position i in the sequence of number groups S1S2...Sk. Each group Sk consists of a sequence of positive i...
2018-04-18 17:18:20
367
原创 51Nod - 1267
给出N个整数,你来判断一下是否能够选出4个数,他们的和为0,可以则输出"Yes",否则输出"No"。Input第1行,1个数N,N为数组的长度(4 <= N <= 1000) 第2 - N + 1行:Aii(-10^9 <= Aii <= 10^9)Output如果可以选出4个数,使得他们的和为0,则输出"Yes",否则输出"No"。Sample Input5-11-
2018-04-11 20:13:45
222
原创 51Nod - 1183
编辑距离,又称Levenshtein距离(也叫做Edit Distance),是指两个字串之间,由一个转成另一个所需的最少编辑操作次数。许可的编辑操作包括将一个字符替换成另一个字符,插入一个字符,删除一个字符。例如将kitten一字转成sitting:sitten (k->s)sittin (e->i)sitting (->g)所以kitten和sitting的编辑距离是3。俄罗...
2018-04-11 10:45:24
263
原创 POJ - 2109
看到1<=p<10101 我就去想大数操作了,后来看了discuss原来double完全可以放。类型 长度 (bit) 有效数字 绝对值范围float 32 6~7 10^(-37) ~ 10^38double 6...
2018-04-09 16:41:38
328
原创 二进制枚举
某君有 nn 个互不相同的正整数,现在他要从这 nn个正整数之中无重复地选取任意个数,并仅通过加法凑出整数 XX。求某君有多少种不同的方案来凑出整数 XX。输入格式第一行,输入两个整数 n,X(1 \leq n \leq 20, 1 \leq X \leq 2000)n,X(1≤n≤20,1≤X≤2000)。接下来输入 nn 个整数,每个整数不超过 100100。输出格式输出一个整数,表示能凑出 ...
2018-04-09 10:48:40
352
ACM数学+数论大全
2019-03-12
ACM计算几何大全
2019-03-12
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人