Mathematics
文章平均质量分 51
Irish_Moonshine
沐月
展开
-
PAT (Advanced Level) 1128 N Queens Puzzle (20分)【化简】
PAT (Advanced Level) 1128 N Queens Puzzle (20分)首先声明,这个题直接双循环暴力就可以过,但是我觉得这样就没什么意思了。复杂度已经超过两亿 200∗1000∗1000200*1000*1000200∗1000∗1000 ,O(N2∗K)O({N}^{2}*K)O(N2∗K) 可能是浙大的OJ比较强劲。本算法复杂度 O(NK)O(NK)O(NK)首先解决的问题就是保证每行与每列的元素唯一,给出的数据就已经是行唯一了,对于列开一个标记数组解决一下。对于列上原创 2020-08-06 21:52:11 · 168 阅读 · 0 评论 -
Hihocoder 编程练习赛69 C 阶乘问题【质因数分解】
时间限制:20000ms 单点时限:1000ms 内存限制:256MB 描述 给定 n, k,求一个最大的整数 m,使得 km 是 n! 的约数输入 第一行两个正整数 n, k2 ≤ n,k ≤ 109输出 输出最大的 m样例输入 5 2 样例输出 3http://hihocoder.com/contest/offers69/problem/3打素数表表打...原创 2018-07-22 14:27:44 · 356 阅读 · 0 评论 -
HDU 4565 So Easy!【共轭构造+矩阵幂】
http://acm.hdu.edu.cn/showproblem.php?pid=4565注意取值范围: 0<a,m<215,(a−1)2<b<a2,0<b,n<2310<a,m<215,(a−1)2<b<a2,0<b,n<2310< a, m < 2^{15}, (a-1)^2< b < a^2, 0 <原创 2018-08-03 16:46:17 · 158 阅读 · 0 评论 -
牛客网暑期ACM多校训练营(第三场)H Diff-prime Pairs【规律】
https://www.nowcoder.com/acm/contest/141/H#include <iostream>#include <string.h>#include <stdio.h>using namespace std;const int maxn = 1e7 + 10;int p[maxn], is[maxn];int n;v...原创 2018-07-26 19:21:32 · 291 阅读 · 0 评论 -
牛客网暑期ACM多校训练营(第一场)B Symmetric Matrix【建模】
https://www.nowcoder.com/acm/contest/139/B根据题目中给出的条件,可以把这个问题抽象成所有点度为2的图的种类数量。 如果度都为2,那么必然成环,有一个或多个环的可能。点为iii的时候,方程可以写成 dp[i][j]=(i−1)∗dp[i−1]+sigma(k<n−2)(n−1)!f(k)/k!/2dp[i][j]=(i−1)∗dp[i−1]...原创 2018-07-25 16:20:55 · 216 阅读 · 0 评论 -
HDU 2091 空心三角形【water】
Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 49657 Accepted Submission(s): 13458Problem Description把一个字符三角形掏空,就能节省材料成本,减轻重量,但关键是为了追求另一种视觉原创 2017-11-20 00:06:57 · 226 阅读 · 0 评论 -
ACdream 1069 无耻的出题人【规律】
听到X神要参加比赛,只会Fibnacci数的出题人被吓得哭晕在厕所。为了防止X神AK(ALL KILL)比赛题目,无耻的出题人只好在题面上做些手脚(加密)。其中一道题的题目描述如下:hjxh dwh v vxxpde,mmo ijzr yfcz hg pbzrxdvgqij rid stl mc zspm vfvuu vb uwu spmwzh.比赛的时候X神对题面研究了四个多小时,终于发现了一丝规原创 2017-12-01 21:45:08 · 401 阅读 · 0 评论 -
上海高校金马五校赛 C 序列变换【全排列】
链接:https://www.nowcoder.com/acm/contest/91/C 来源:牛客网时间限制:C/C++ 2秒,其他语言4秒 空间限制:C/C++ 131072K,其他语言262144K 64bit IO Format: %lld 题目描述 给定两个长度为n的序列,ai, bi(1<=i<=n), 通过3种魔法使得序列a变换为序列b,也就是ai=bi(...原创 2018-05-04 00:26:49 · 282 阅读 · 0 评论 -
第13届景驰-埃森哲杯广东工业大学ACM程序设计大赛 F 等式【因数分解】
链接:https://www.nowcoder.com/acm/contest/90/F 来源:牛客网时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 32768K,其他语言65536K 64bit IO Format: %lld 题目描述 给定n,求1/x + 1/y = 1/n (x<=y)的解数。(x、y、n均为正整数)输入描述: 在第一行输入一个正...原创 2018-04-04 12:15:27 · 178 阅读 · 0 评论 -
CF 38 D Vasya the Architect【递推+重心公式】
Once Vasya played bricks. All the bricks in the set had regular cubical shape. Vasya vas a talented architect, however the tower he built kept falling apart.Let us consider the building process. Vas...原创 2018-04-03 22:09:50 · 258 阅读 · 0 评论 -
第十四届浙江财经大学程序设计竞赛 B Bazinga【循环结+规律+快速幂】
链接:https://www.nowcoder.com/acm/contest/89/B 来源:牛客网 这里写图片描述 时间限制:C/C++ 3秒,其他语言6秒 空间限制:C/C++ 262144K,其他语言524288K 64bit IO Format: %lld 题目描述 There is a formula: gcd(a, b) is the greatest comm...原创 2018-03-28 12:52:58 · 361 阅读 · 0 评论 -
康托展开【Template】
#include<iostream>#include<cstring>#include<cstdio>#include<algorithm>#include<vector>#include<cmath>#include<map>#include<set>#include&原创 2018-03-28 12:44:13 · 190 阅读 · 0 评论 -
第十四届浙江财经大学程序设计竞赛 I Interesting Set【组合数】
链接:https://www.nowcoder.com/acm/contest/89/I 来源:牛客网时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 262144K,其他语言524288K 64bit IO Format: %lld 题目描述 Mr.Frog is researching integers. He converts two integers X a...原创 2018-03-28 12:43:38 · 230 阅读 · 0 评论 -
上下火车【二分+矩阵快速幂】
问题 J: 上下火车 时间限制: 1 Sec 内存限制: 125 MB 提交: 28 解决: 16 [提交][状态][讨论版][命题人:外部导入] 题目描述 火车从始发站(称为第1站)开出,在始发站上车的人数为a,然后到达第2站,在第2站有人上、下车,但上、下车的人数相同,因此在第2站开出时(即在到达第3站之前)车上的人数保持为a人。从第3站起(包括第3站)上、下车的人数有一定规律:...原创 2018-03-17 23:24:54 · 268 阅读 · 0 评论 -
NOIP2005普及组第4题 循环【高精度+规律】
问题 C: NOIP2005普及组第4题 循环 时间限制: 1 Sec 内存限制: 128 MB 提交: 4 解决: 2 [提交][状态][讨论版][命题人:外部导入] 题目描述 乐乐是一个聪明而又勤奋好学的孩子。他总喜欢探求事物的规律。一天,他突然对数的正整数次幂产生了兴趣。 众所周知,2的正整数次幂最后一位数总是不断的在重复2,4,8,6,2,4,8,6……我们说2的正整...原创 2018-03-15 21:47:58 · 1446 阅读 · 0 评论 -
牛客网暑期ACM多校训练营(第三场)J Distance to Work【二分+圆交】
https://www.nowcoder.com/acm/contest/141/J给出一个多边形,多边形内区域为城市。 给出几个工作地点,一个工作地点给出一个比例p/q,表示相比于所有的可能,我这个工作范围恰好站总生活面积的p/q。我们可以二分圆的半径,使公共部分面积=p/q∗p/q∗p/q*城市面积。圆与多边形交可以套用板子。#include<bits/stdc++.h...原创 2018-07-27 13:34:09 · 194 阅读 · 0 评论 -
牛客网暑期ACM多校训练营(第二场)I car【规律】
https://www.nowcoder.com/acm/contest/140/I用四个数组存下这张图。 策略为交错放车,奇偶要特殊考虑。#include<bits/stdc++.h>using namespace std;#define ll long longconst int N=1e5+10;const int inf=2147483647;int L[N...原创 2018-07-27 18:12:46 · 178 阅读 · 0 评论 -
洛谷 P2257 YY的GCD【莫比乌斯反演】
https://www.luogu.org/problemnew/show/P2257一年前留下的莫比乌斯反演的坑,竟发现一年后还是不懂(这不是废话嘛)!但是我觉得搞一搞还是可以的。实际上算是第一次做莫比乌斯反演问题。初步感觉解决这类问题的办法就是推公式,当然不是说是这么简单,日后发现有问题还将更改。一般来说将一个枚举求和问题转换成一个整数分块求前缀和的问题。FFF与fff的关...原创 2018-08-14 06:18:04 · 290 阅读 · 0 评论 -
2018 Multi-University Training Contest 8 1001 Character Encoding【容斥】
http://acm.hdu.edu.cn/contests/contest_showproblem.php?pid=1001&amp;cid=809题意:每个数字的取值范围为0到n-1,共有m个数字,求总和为k的方案数。因为k的上限为1e5,那么我们不妨把k看做k个1,然后通过m-1个隔板来分割出m个数字。但是这么做会有一个问题,那就是数字可能为0。那么不妨将数字的取值范围改成1-n,m...原创 2018-08-16 05:21:44 · 398 阅读 · 0 评论 -
2018 Multi-University Training Contest 7 1010 Sequence【整数分块+矩阵幂】
http://acm.hdu.edu.cn/contests/contest_showproblem.php?pid=1010&cid=808题意:在扩展斐波纳挈的基础上加了一个变量P/nP/nP/n。求第nnn项的取值。分析:考虑将每一种P/nP/nP/n进行矩阵快速幂,也就是进行了整数分块处理。对于每一个整数块可以使用矩阵快速幂,然后维护A,BA,BA,B用作下一次的矩阵快速幂...原创 2018-08-13 18:34:53 · 274 阅读 · 0 评论 -
HDU 4556 Stern-Brocot Tree【法里数列,欧拉函数】
http://acm.hdu.edu.cn/showproblem.php?pid=4556只观察树的左半边会发现全是真分数,每一层就是分母不超过1-n的所有真分数的表示。这种树是用来构造法里数列的;法里数列表示的是分母不超过n的所有真分数的数量。也可以理解成是欧拉函数的前缀和。我们可以考虑用类似素数筛的方式来快速求得欧拉函数的值。#include<bits/stdc...原创 2018-08-03 14:32:39 · 472 阅读 · 0 评论 -
牛客网暑期ACM多校训练营(第八场)G Counting regions【几何公式推导】
https://www.nowcoder.com/acm/contest/146/G公式推导参考自:凸n边形的对角线最多能将其内部分成几个区域。(n−1)(n−2)(n2−3n+12)/24(n−1)(n−2)(n2−3n+12)/24(n-1)(n-2)(n^2-3n+12)/24#include<iostream>#include<cstring>#in...原创 2018-08-12 05:33:41 · 716 阅读 · 0 评论 -
2018 Multi-University Training Contest 6 L Pinball【受力分析】
http://acm.hdu.edu.cn/contests/contest_showproblem.php?cid=807&amp;pid=1012将物体的重力加速度分解。 一部分用来做弹起运动,一部分用来做匀加速直线运动。 可以计算得出两个方向上的运动周期,两者之商为答案。#include &lt;iostream&gt;#include &lt;stdio.h&gt;#inc...原创 2018-08-12 05:05:23 · 440 阅读 · 0 评论 -
2018 Multi-University Training Contest 6 A oval-and-rectangle【精度】
http://acm.hdu.edu.cn/contests/contest_showproblem.php?pid=1001&cid=807用定积分推导出答案是 2.0*b+a*pi 但是答案要求忽略掉第七位,所以不可以直接用double(四舍五入)。 要么转换成长整型,要么将double第七位的隐患解决掉(感觉可以二分,很啰嗦但很有效)。#include <iostr...原创 2018-08-12 04:52:50 · 308 阅读 · 0 评论 -
2015 ACM/ICPC Asia Regional Shenyang Online & HDU 5451 Best Solver【循环节+数论||共轭矩阵幂+广义斐波那契降幂】
http://acm.hdu.edu.cn/showproblem.php?pid=5451通过一道名为So easy!的共轭构造矩阵幂取模的问题来联想到这个题,不难发现应该用同样的套路来进行共轭构造。但是问题来了,指数过分强大,取模数超小,于是一开始考虑欧拉降幂再矩阵快速幂。数据一大,发现答案全错!后请教得知无理数无法欧拉降幂。那么就只有一种可能了——有规律可循!找循环节,其实欧...原创 2018-08-12 04:23:34 · 406 阅读 · 0 评论 -
牛客网暑期ACM多校训练营(第五场)J plan【分类】
https://www.nowcoder.com/acm/contest/143/J答案无非那么几种,要么都用2,要么都用3,要么3主2补,要么2主3补。#include<iostream>#include<cstring>#include<cstdio>#include<algorithm>#include<iomanip&g...原创 2018-08-12 03:20:17 · 298 阅读 · 0 评论 -
牛客网暑期ACM多校训练营(第五场)G max【数学】
https://www.nowcoder.com/acm/contest/143/G给出一个ccc,给出一个取值范围[1,n][1,n][1,n]。在范围中取两个数a,ba,ba,b使得 gcd(a,b)=cgcd(a,b)=cgcd(a,b)=c,最大化 a∗ba∗ba*b。答案相当于是c∗c∗x∗yc∗c∗x∗yc*c*x*y 其中 x,yx,yx,y 互质,那么最好的办法就是两者相差...原创 2018-08-12 03:09:21 · 307 阅读 · 0 评论 -
BZOJ 4805 欧拉函数求和【杜教筛】
https://www.lydsy.com/JudgeOnline/problem.php?id=4805求欧拉函数的前缀和,项数小于2e9。以目前的视野来看待杜教筛的话,感觉就像是将一个线性的式子,进一步优化,然后通过记忆化搜索来实现的一个过程。S(n)=∑i=1nϕ(i)S(n)=∑i=1nϕ(i)S(n)=\sum^{n}_{i=1}\phi(i) 然后卷积一下,通过套路式...原创 2018-08-14 16:05:30 · 321 阅读 · 0 评论 -
2018 Multi-University Training Contest 7 6390 GuGuFishtion【莫比乌斯反演】
http://acm.hdu.edu.cn/showproblem.php?pid=6390昨天一下午,今天一上午,我已经快被这个题恶心死了。 TLE,TLE永远都是TLE。不得不说,我觉得卡我卡的没什么道理,非常的不爽!卡内循环枚举变量类型我是真的不懂,发现知识盲区+1.可是就算是再不爽也要含泪挖坑。题意求: (∑a=1m∑b=1nϕ(ab)ϕ(a)ϕ(b))mod(p)(∑...原创 2018-08-14 13:33:04 · 321 阅读 · 0 评论 -
2018 Multi-University Training Contest 5 1005 Everything Has Changed【圆交】
http://acm.hdu.edu.cn/contests/contest_showproblem.php?pid=1005&cid=806题意:给出一个母圆,给出m个子圆来切割子圆,求最后母圆被切割后的周长。题目保证所有的子圆互不相交。如果子圆内切与母圆,那么他对于答案的贡献是他的圆周长。 如果子圆相离、内含、外切与母圆,贡献为0;注意有优弧与劣弧。#includ...原创 2018-08-06 18:05:14 · 180 阅读 · 0 评论 -
洛谷 P3327 [SDOI2015]约数个数和【莫比乌斯反演】
https://www.luogu.org/problemnew/show/P3327前置技能: d(ij)=∑x|i∑y|j[gcd(x,y)=1]d(ij)=∑x|i∑y|j[gcd(x,y)=1]d(ij)=\sum^{}_{x|i}\sum^{}_{y|j}[gcd(x,y)=1] 然后按照套路一顿乱推……好吧我是没推出来= = 最终结果: Ans=∑d=1min(n,m)μ(...原创 2018-08-14 08:30:51 · 282 阅读 · 0 评论 -
洛谷 P3455 [POI2007]ZAP-Queries【莫比乌斯反演】
https://www.luogu.org/problemnew/show/P3455这个题的题意与 【洛谷 P2257 YY的GCD】挺相似的,只不过现在要求去gcdgcdgcd为固定值的答案的数量。由于很像,公式还是很好推的:Ans=∑t=1min(a,b)μ(t)⌊atd⌋⌊btd⌋Ans=∑t=1min(a,b)μ(t)⌊atd⌋⌊btd⌋Ans=\sum^{min(a,b...原创 2018-08-14 06:41:04 · 212 阅读 · 0 评论 -
2017浙工大院赛预赛 J 裁缝大师【几何】
时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 131072K,其他语言262144K 64bit IO Format: %lld 题目描述 坤酱想把一块圆形的布裁成正多边形,于是请你告诉坤酱正多边形的几个顶点应在哪里? 为了方便表示,圆给出在坐标系中,正多边形的第一个顶点固定在该圆在平行于x轴正方向最远的位置上,请按顺时针顺序输出所有的顶点。 输入描述: 输入第一行给原创 2017-12-16 23:50:30 · 508 阅读 · 2 评论 -
2017 浙工业院赛预赛 C【推导+二分+逆元】
时间限制:C/C++ 3秒,其他语言6秒 空间限制:C/C++ 131072K,其他语言262144K 64bit IO Format: %lld 题目描述 给定一个序列数列 , (ai互不相同) 3种映射关系现在对于给定询问(x, k),求输入描述: 多组读入数据T, 对于每组数据, 第1行一个整数n,q,n表示数列的大小,q为询问数 第2行读入n个数a1,a2,…,an,表示原创 2017-12-18 23:54:48 · 239 阅读 · 0 评论 -
POJ 2356 Find a multiple【鸽巢原理】
Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 8457 Accepted: 3673 Special Judge DescriptionThe input contains N natural (i.e. positive integer) numbers ( N <= 10000 ). Each原创 2017-11-14 18:10:48 · 371 阅读 · 0 评论 -
HDU 1085 Holding Bin-Laden Captive!【母函数||多重背包】
题意:给你cnt[0]个一元硬币,cnt[1]个两元硬币,cnt[2]个五元硬币,问不能凑出来的第一个面额是多少。分析:不得不说遇到这个题的第一反应是多重背包265ms。一直跑到面额最大的数字,显然肯定能解决。另一种则是抽象成生成函数(母函数)加以解决64ms。母函数模拟方式: (1+x+x2+x3+.........xcnt1)∗(1+x2+x4+x6+.........xcnt2)∗(1+x5原创 2017-11-28 16:31:27 · 269 阅读 · 0 评论 -
G - Lining Up【几何】
“How am I ever going to solve this problem?” said the pilot. Indeed, the pilot was not facing an easy task. She had to drop packages at specific points scattered in a dangerous area. Furthermore, the p原创 2017-11-13 21:24:41 · 271 阅读 · 0 评论 -
A - A Funny Game 【博弈】
Alice and Bob decide to play a funny game. At the beginning of the game they pick n(1 <= n <= 10 6) coins in a circle, as Figure 1 shows. A move consists in removing one or two adjacent coins, leaving原创 2017-11-13 21:23:22 · 488 阅读 · 0 评论 -
蝴蝶定理【几何】
蝴蝶定理的英文是Butterfly Theorem,蝴蝶定理是古典欧氏平面几何的最精彩的结果之一。蝴蝶定理(Butterfly Theorem):设M为圆内弦PQ的中点,过M作弦AB和CD。设AD和BC各相交PQ于点X和Y,则M是XY的中点。 去掉中点的条件,结论变为一个一般关于有向线段的比例式,称为“坎迪定理”, 不为中点时满足:1/MY-1/MX=1/MQ-1/MP ,这对2,3均成立。证明:原创 2017-10-22 14:40:48 · 2541 阅读 · 0 评论 -
CF 332 B.Maximum Absurdity【DP】
Reforms continue entering Berland. For example, during yesterday sitting the Berland Parliament approved as much as n laws (each law has been assigned a unique number from 1 to n). Today all these laws原创 2017-10-21 17:38:42 · 524 阅读 · 0 评论