
===数学相关===
文章平均质量分 67
DQSSS
蒟蒻
展开
-
【codevs10月月赛(第三次)】day 1
codevs第三次月赛T1正解,T2部分分,T3迭代深搜…… 结果:T1:100,T2:40,T3:40 rank25…其实T3的想法挺接近70分算法的,然而我还是打的暴力…还是太弱T1:Cww的作业打表可得:fn=n/2(下取整)于是答案是2∑k^2(1<=k<=n)%10007公式可知,答案是(n(n+1)(2n+1)/3) %10007带除法的模,要求3的逆元,也就是求3x=1(%1000原创 2015-10-06 11:18:42 · 856 阅读 · 0 评论 -
【bzoj2440】[中山市选2011]完全平方数 莫比乌斯反演+容斥原理
Description小 X 自幼就很喜欢数。但奇怪的是,他十分讨厌完全平方数。他觉得这些 数看起来很令人难受。由此,他也讨厌所有是完全平方数的正整数倍的数。然而 这丝毫不影响他对其他数的热爱。 这天是小X的生日,小 W 想送一个数给他作为生日礼物。当然他不能送一 个小X讨厌的数。他列出了所有小X不讨厌的数,然后选取了第 K个数送给了 小X。小X很开心地收下了。 然而现在小 W 却记原创 2016-01-20 20:05:27 · 648 阅读 · 0 评论 -
【poj3233】Matrix Power Series 矩阵+快速幂
DescriptionGiven a n × n matrix A and a positive integer k, find the sum S = A + A2 + A3 + … + Ak.InputThe input contains exactly one test case. The first line of input contains three positive integers原创 2016-01-11 20:42:14 · 585 阅读 · 0 评论 -
【bzoj2186】[Sdoi2008]沙拉公主的困惑 数论 线性筛逆元
Description 大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为1到N的阶乘,但是,政府只发行编号与M!互质的钞票。房地产第一大户沙拉公主决定预测一下大富翁国现在所有真钞票的数量。现在,请你帮助沙拉公主解决这个问题,由于可能张数非常大,你只需计算出对R取模后的答案即可。R是一个质数。Input第一行为两个整数T,R。R<=10^9+10,T<=10000,表原创 2016-01-14 21:19:53 · 912 阅读 · 0 评论 -
【bzoj1257】[CQOI2007]余数之和sum 数论乱搞
Description给出正整数n和k,计算j(n, k)=k mod 1 + k mod 2 + k mod 3 + … + k mod n的值,其中k mod i表示k除以i的余数。例如j(5, 3)=3 mod 1 + 3 mod 2 + 3 mod 3 + 3 mod 4 + 3 mod 5=0+1+0+3+3=7Input输入仅一行,包含两个整数n, k。Output输出仅一行,即j(n,原创 2016-01-15 12:47:56 · 1131 阅读 · 0 评论 -
【bzoj2190】[SDOI2008]仪仗队 线性筛欧拉函数
Description 作为体育委员,C君负责这次运动会仪仗队的训练。仪仗队是由学生组成的N * N的方阵,为了保证队伍在行进中整齐划一,C君会跟在仪仗队的左后方,根据其视线所及的学生人数来判断队伍是否整齐(如下图)。 现在,C君希望你告诉他队伍整齐时能看到的学生人数。Input 共一个数N。Output 共一个数,即C君应看到的学生人数。Sample Input4Sampl原创 2016-01-15 14:52:04 · 590 阅读 · 0 评论 -
【bzoj1013】[JSOI2008]球形空间产生器sphere 高斯消元
Description有一个球形空间产生器能够在n维空间中产生一个坚硬的球体。现在,你被困在了这个n维球体中,你只知道球面上n+1个点的坐标,你需要以最快的速度确定这个n维球体的球心坐标,以便于摧毁这个球形空间产生器。Input第一行是一个整数,n。接下来的n+1行,每行有n个实数,表示球面上一点的n维坐标。每一个实数精确到小数点后6位,且其绝对值都不超过20000。Output有且只有一行,依次给原创 2016-03-16 20:03:17 · 634 阅读 · 0 评论 -
【bzoj1923】[Sdoi2010]外星千足虫 高斯消元解xor方程组
DescriptionInput第一行是两个正整数 N, M。 接下来 M行,按顺序给出 Charles 这M次使用“点足机”的统计结果。每行 包含一个“01”串和一个数字,用一个空格隔开。“01”串按位依次表示每只虫 子是否被放入机器:如果第 i 个字符是“0”则代表编号为 i 的虫子未被放入,“1” 则代表已被放入。后面跟的数字是统计的昆虫足数 mod 2 的结果。 由于 NASA的实验机器精确原创 2016-03-17 08:53:19 · 841 阅读 · 0 评论 -
【bzoj2115】[Wc2011] Xor 高斯消元+dfs
DescriptionInput第一行包含两个整数N和 M, 表示该无向图中点的数目与边的数目。 接下来M 行描述 M 条边,每行三个整数Si,Ti ,Di,表示 Si 与Ti之间存在 一条权值为 Di的无向边。 图中可能有重边或自环。Output仅包含一个整数,表示最大的XOR和(十进制结果),注意输出后加换行回车。Sample Input5 7 1 2 2 1 3 2 2 4 1 2 5 1 4原创 2016-03-17 10:26:15 · 619 阅读 · 0 评论 -
【bzoj3994】[SDOI2015]约数个数和 莫比乌斯反演
Description设d(x)为x的约数个数,给定N、M,求 Input输入文件包含多组测试数据。第一行,一个整数T,表示测试数据的组数。 接下来的T行,每行两个整数N、M。OutputT行,每行一个整数,表示你所求的答案。Sample Input27 45 6Sample Output110121HINT1<=N, M<=500001<=T<=50000SourceRound 1 感谢y原创 2016-02-23 21:13:46 · 1191 阅读 · 0 评论 -
CH Round #65 - MFOI杯水题欢乐赛 day2 树形DP+SG函数+线段树合并+DP
题目太长不粘贴了,链接在这: CH Round #65 - MFOI杯水题欢乐赛 day2题目链接 -> 密码:exdb 题解T1算叶子之间的距离算成深度相加忘了减去LCA了,我傻XT2不会,弃题,没学过SG函数T3dfs10分,dp不会,我傻X总结:10分滚粗,人弱就不该作死T1题意:给你两棵树,从a树中选两个不同叶子,b树中选两个不同的叶子,连接起来组成一个环,问所有不同方案这种环的大小总和原创 2016-03-23 10:50:11 · 904 阅读 · 0 评论 -
【bzoj1898】[Zjoi2004]Swamp 沼泽鳄鱼 矩阵乘法+DP
Description潘塔纳尔沼泽地号称世界上最大的一块湿地,它地位于巴西中部马托格罗索州的南部地区。每当雨季来临,这里碧波荡漾、生机盎然,引来不少游客。为了让游玩更有情趣,人们在池塘的中央建设了几座石墩和石桥,每座石桥连接着两座石墩,且每两座石墩之间至多只有一座石桥。这个景点造好之后一直没敢对外开放,原因是池塘里有不少危险的食人鱼。豆豆先生酷爱冒险,他一听说这个消息,立马赶到了池塘,想做第一个在桥原创 2016-03-14 08:57:12 · 723 阅读 · 0 评论 -
【bzoj3329】Xorequ 矩阵乘法+数位DP
DescriptionInput第一行一个正整数,表示数据组数据 ,接下来T行 每行一个正整数NOutput2*T行 第2*i-1行表示第i个数据中问题一的解,第2*i行表示第i个数据中问题二的解,Sample Input11Sample Output12HINTx=1与x=2都是原方程的根,注意第一个问题的解不要mod 10^9+71<=N<=10^18 1<=T<=1000SourceB原创 2016-03-15 14:18:52 · 1198 阅读 · 0 评论 -
【bzoj1009】[HNOI2008]GT考试 矩阵+kmp+DP
Description阿申准备报名参加GT考试,准考证号为N位数X1X2….Xn(0<=Xi<=9),他不希望准考证号上出现不吉利的数字。他的不吉利数学A1A2…Am(0<=Ai<=9)有M位,不出现是指X1X2…Xn中没有恰好一段等于A1A2…Am. A1和X1可以为0Input第一行输入N,M,K.接下来一行输入M位的数。 100%数据N<=10^9,M<=20,K<=1000 40%数据N<=原创 2016-03-15 16:33:33 · 2484 阅读 · 2 评论 -
【tsinsen A1490】osu!(乔明达) 矩阵+线段树
试题来源 2013中国国家集训队第二次作业问题描述 osu!是一个基于《押忍!战斗!应援团》《精英节拍特工》《太鼓达人》等各种音乐游戏做成的一款独特的PC版音乐游戏。游戏中,玩家需要根据音乐的节奏,通过鼠标点击或敲击按键合成一首歌曲。 一张osu!的地图是由若干个“音”排列而成的。在本题中,对于每个音我们只需要考虑成功点击和错过(miss)这两种情况。对于一张osu!地图,玩家的完成情况可原创 2016-04-25 20:15:09 · 1602 阅读 · 1 评论 -
【弱校胡策】2016.4.25 (bzoj4108、tsinsenA1339、tsinsenA1490)上下界网络流+树状数组+线段树+矩阵乘法+概率与期望
出题人说T1是大水题,于是我放弃了几乎一眼的上下界网络流,开始想MST 几个小时之后然后出题人告诉我正解是上下界网络流 卧槽说好的大水题呢 T1题意好像不对2333 T3竟然是期望!暴力我都不会 T2一眼看去只会特殊数据,好吧我就写了这些。后来发现暴力可以写一下?算了算了弃疗2333fqk毒瘤出题人!!雪题意一个序列,每个点只能往它后面的点走,有一个费用,要求走K次不重复走完所有点,求最小花费原创 2016-04-25 20:53:31 · 903 阅读 · 0 评论 -
【弱校胡策】2016.4.14 (bzoj2164)最短路+状压DP+矩阵乘法+高斯消元+树链剖分+线段树+背包DP
cyyz&qhyz&lwyz&gryz弱校胡策 命题人:cyyz ws_fqkT3暴力写挫了 50+10+0滚粗辣!奇妙的约会(appointment.cpp/c/pas)【问题描述】DQS和sxb在网上结识后成为了非常好的朋友,并且都有着惊人 的OI水平。在NOI2333的比赛中,两人均拿到了金牌,并保送进入 HU/PKU。于是两人决定在这喜大普奔的时刻进行面基。 NOI2333参赛选手众多原创 2016-04-14 22:12:08 · 5058 阅读 · 7 评论 -
【弱校胡策】2016.4.19 LCA+LCT+莫比乌斯反演+SAM+启发式合并
弱校胡策题解 命题人:Loi_DQS 2016.4.19前言来自出题人的吐槽:T1的题目来源是去年十月份做NOIP模拟题和lcyz(聊城一中)胡策(其实也不算胡策,从他们那里要的题)的T3,T2是去年五月份学长带着我们在tyvj举办的有奖赛(http://www.tyvj.cn/Contest/187 and http://www.tyvj.cn/Contest/192)的某题。T3是上周六(20原创 2016-04-19 20:24:15 · 2514 阅读 · 4 评论 -
【bzoj2694】Lcm 莫比乌斯反演
Description对于任意的>1的n gcd(a, b)不是n^2的倍数 也就是说gcd(a, b)没有一个因子的次数>=2Input一个正整数T表示数据组数 接下来T行 每行两个正整数 表示N、MOutputT行 每行一个整数 表示第i组数据的结果Sample Input42 43 36 58 3Sample Output2428233178HINTHINTT <= 10000N, M原创 2016-01-25 21:35:13 · 1468 阅读 · 0 评论 -
【codevs4355】王的对决(简单数论) 莫比乌斯反演
题目描述 Description/*在神秘而遥远的T国度生活着一群神奇的物种——DD!DD们有着独特的生活习惯,他们总是喜欢居住在tree上,而且只有上层住了DD下层才能再住。并且每层住的DD最多不超过2n2个。。。这一天,DD跳下来了!(也就是说刚才的条件一点用都没有~_~)*/他们遇到了一只自称为S的dd。但是本身性质为S的DD们不服,于是他们决定进行摔跤,摔跤获胜者成为王!(由于DD和dd们的原创 2016-01-25 15:29:26 · 1152 阅读 · 0 评论 -
【NOI2012】随机数生成器
bzoj2875第一眼看是矩阵乘法+高精,顺便借此机会学了一下压八位高精然后一直WA…拍到最后拍不下去了,索性重打了一遍,不过没用高精,用的类似快速幂的快速乘,同样解决了溢出问题~代码: (快速乘+矩阵乘法)原创 2015-09-19 21:06:46 · 853 阅读 · 0 评论 -
【Smart OJ P2485】Sum of product 简单数学
描述 Descript.对于A1,A2,…..,AN,求 的值。输入 Input第1 行,1 个整数N。 第2 行,N 个整数A1,A2,…..,AN。输出 Output1 个整数,表示所求的值。样例 Sample输入数据31 2 3输出数据36备注 Hint• 对于30% 的数据,3 <= N <= 500; • 对于60% 的数据,3 <= N <= 5000; •原创 2015-10-28 19:23:09 · 1025 阅读 · 0 评论 -
【bzoj1008】越狱 组合数学
Description监狱有连续编号为1…N的N个房间,每个房间关押一个犯人,有M种宗教,每个犯人可能信仰其中一种。如果相邻房间的犯人的宗教相同,就可能发生越狱,求有多少种状态可能发生越狱Input输入两个整数M,N.1<=M<=10^8,1<=N<=10^12Output可能越狱的状态数,模100003取余Sample Input2 3Sample Output6HINT6种状态为(000)(原创 2015-10-29 13:38:26 · 1271 阅读 · 2 评论 -
【lcyz】noip模拟题
题目来自友好外交的聊城一中~由于我是自己闲的做的,LOI做的人不多……战绩:T1:100,T2:20(数组开小+lld),T3:0(输入有问题【据说?】+lld)后来把T2T3傻X错误改了然后AK了,不得不说眼瞎真是害死人!!注意!!!这里的代码是我120分的代码!!请注意T3的输入【据说?】!! 还有就是不要像我一样傻X全开longlong,会T!!! 出题人语文我给十分!!原创 2015-10-15 22:00:16 · 734 阅读 · 0 评论 -
【codevs1060】搞笑世界杯 DP+数学
题目描述 Description随着世界杯小组赛的结束,法国,阿根廷等世界强队都纷纷被淘汰,让人心痛不已. 于是有人组织了一场搞笑世界杯,将这些被淘汰的强队重新组织起来和世界杯一同比赛.你和你的朋友欣然去购买球票.不过搞笑世界杯的球票出售方式也很特别,它们只准备了两种球票.A 类票——免费球票 B 类票——-双倍价钱球票.购买时由工作人员通过掷硬币决定,投到正面的买A类票, 反面的买B类票.并且由原创 2015-10-14 16:38:12 · 936 阅读 · 0 评论 -
【NOIP2006普及】数列 快速幂+找规律
题目描述 Description给定一个正整数k(3≤k≤15),把所有k的方幂及所有有限个互不相等的k的方幂之和构成一个递增的序列,例如,当k=3时,这个序列是:1,3,4,9,10,12,13,…(该序列实际上就是:30,31,30+31,32,30+32,31+32,30+31+32,…)请你求出这个序列的第N项的值(用10进制数表示)。例如,对于k=3,N=100,正确答案应该是981。输入原创 2015-10-29 21:00:27 · 1270 阅读 · 2 评论 -
【NOIP2014】解方程 玄学哈希+乱搞
我觉得我考的话最多敲个高精了……n=0 => n%p=0,显然是对的。n%p=0 => n=0,显然是错的。但p多选几个,就极有可能是对的……所以说我们可以对输入给定的数模好几个质数,然后枚举1~m看看有几个能使方程成立,这样能拿70 。满分做法:n%p=(n+kp)%p,所以每次只枚举[0,p-1],然后把倍数筛掉即可。p不能选太大,否则会T。p不能选太小,否则会WA。p不能选太少,否则筛不干净…原创 2015-11-02 15:22:38 · 1275 阅读 · 0 评论 -
【bzoj3239】Discrete Logging BSGS
DescriptionGiven a prime P, 2 <= P < 231, an integer B, 2 <= B < P, and an integer N, 2 <= N < P, compute the discrete logarithm of N, base B, modulo P. That is, find an integer L such thatBL == N (mod原创 2016-01-15 21:36:53 · 674 阅读 · 0 评论 -
【bzoj2242】[SDOI2011]计算器 快速幂+exgcd+BSGS
Description你被要求设计一个计算器完成以下三项任务: 1、给定y,z,p,计算Y^Z Mod P 的值; 2、给定y,z,p,计算满足xy≡ Z ( mod P )的最小非负整数; 3、给定y,z,p,计算满足Y^x ≡ Z ( mod P)的最小非负整数。Input输入包含多组数据。第一行包含两个正整数T,K分别表示数据组数和询问类型(对于一个测试点内的所有数据,询问类型相同)。原创 2016-01-15 21:02:28 · 639 阅读 · 0 评论 -
【bzoj2705】[SDOI2012]Longge的问题 欧拉函数
DescriptionLongge的数学成绩非常好,并且他非常乐于挑战高难度的数学问题。现在问题来了:给定一个整数N,你需要求出∑gcd(i, N)(1<=i <=N)。Input一个整数,为N。Output一个整数,为所求的答案。Sample Input6Sample Output15HINT【数据范围】对于60%的数据,0Sourceround1 day1先贴一个hzwer的,感觉说的好有道原创 2016-01-19 16:01:58 · 590 阅读 · 0 评论 -
【bzoj1101】[POI2007]Zap 莫比乌斯反演
DescriptionFGD正在破解一段密码,他需要回答很多类似的问题:对于给定的整数a,b和d,有多少正整数对x,y,满足x<=a,y<=b,并且gcd(x,y)=d。作为FGD的同学,FGD希望得到你的帮助。Input第一行包含一个正整数n,表示一共有n组询问。(1<=n<= 50000)接下来n行,每行表示一个询问,每行三个正整数,分别为a,b,d。(1<=d<=a,b<=50000)Outp原创 2016-01-19 19:38:52 · 803 阅读 · 0 评论 -
【codevs3223】素数密度 埃氏筛法
题目描述 Description给定区间[L, R](L <= R <= 2147483647,R-L <= 1000000),请计算区间中素数的个数。输入描述 Input Description两个数L和R输出描述 Output Description一行,区间中素数的个数样例输入 Sample Input2 11样例输出 Sample Output5数据范围及提示 Data Size & H原创 2015-10-29 08:09:12 · 1178 阅读 · 1 评论 -
【bzoj2154/2693】Crash的数字表格/jzptab 莫比乌斯反演
Description今天的数学课上,Crash小朋友学习了最小公倍数(Least Common Multiple)。对于两个正整数a和b,LCM(a, b)表示能同时被a和b整除的最小正整数。例如,LCM(6, 8) = 24。回到家后,Crash还在想着课上学的东西,为了研究最小公倍数,他画了一张N*M的表格。每个格子里写了一个数字,其中第i行第j列的那个格子里写着数为LCM(i, j)。一个4原创 2016-01-20 11:35:49 · 1003 阅读 · 0 评论 -
【bzoj2301】[HAOI2011]Problem b 莫比乌斯反演
Description对于给出的n个询问,每次求有多少个数对(x,y),满足a≤x≤b,c≤y≤d,且gcd(x,y) = k,gcd(x,y)函数为x和y的最大公约数。Input第一行一个整数n,接下来n行每行五个整数,分别表示a、b、c、d、kOutput共n行,每行一个整数表示满足要求的数对(x,y)的个数Sample Input22 5 1 5 11 5 1 5 2Sample Outpu原创 2016-01-20 16:26:22 · 721 阅读 · 0 评论 -
【bzoj2818】Gcd 欧拉函数
Description给定整数N,求1<=x,y<=N且Gcd(x,y)为素数的 数对(x,y)有多少对.Input一个整数NOutput如题Sample Input4Sample Output4HINThint对于样例(2,2),(2,4),(3,3),(4,2)1<=N<=10^7Source湖北省队互测题目真是简单粗暴…∑p是素数∑i<=n∑j<=ngcd(i,j)=p\sum_{p是素原创 2016-01-21 10:45:54 · 855 阅读 · 0 评论 -
【bzoj2820】YY的GCD 莫比乌斯反演
Description神犇YY虐完数论后给傻×kAc出了一题 给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对 kAc这种傻×必然不会了,于是向你来请教…… 多组输入Input第一行一个整数T 表述数据组数 接下来T行,每行两个正整数,表示N, MOutputT行,每行一个整数表示第i组数据的结果Sample Input210 10100 10原创 2016-01-21 11:18:35 · 1799 阅读 · 0 评论 -
【bzoj2005】[Noi2010]能量采集 莫比乌斯反演
Description栋栋有一块长方形的地,他在地上种了一种能量植物,这种植物可以采集太阳光的能量。在这些植物采集能量后,栋栋再使用一个能量汇集机器把这些植物采集到的能量汇集到一起。 栋栋的植物种得非常整齐,一共有n列,每列有m棵,植物的横竖间距都一样,因此对于每一棵植物,栋栋可以用一个坐标(x, y)来表示,其中x的范围是1至n,表示是在第x列,y的范围是1至m,表示是在第x列的第y棵。 由于能量原创 2016-01-21 15:01:35 · 1327 阅读 · 0 评论 -
【模板】【笔记】数论各种姿势等等
数论的东西比较烦…最近脑子有点乱,学一点就赶紧记下来,以免忘了……先贴个线性筛:void shai(){ for(int i = 2,tot = 0;i <= n;i ++) { if(!vis[i]) pri[++ tot] = i;//当i是素数时 for(int j = 1,m;j <= tot && (m = i * pri[j]) <= n原创 2016-01-18 20:57:19 · 1325 阅读 · 6 评论