2017年12月_Coco_T_

原创反演中的常见公式

完美的总结都是为了偷懒GCD的计数① ② ③ GCD之和① ② LCM之和① ② n是质数： n不是质数： n的非平方质因子i： n的其他质因子i（即n=i*p且i|p）：其它常见公式① 1~n中所有与n互质的数之和 ② 约数个数 ③ i^2前缀和

2017-12-30 10:27:56 3059

原创 bzoj4659 Lcm（反演）

Description给出A，B，考虑所有满足l<=a<=A，l<=b<=B，且不存在n>1使得n^2同时整除a和b的有序数对(a，b)，求其lcm(a，b)之和。答案模2^30。Input第一行一个整数T表示数据组数。接下来T行每行两个整数A，B表示一组数据。 T ≤ 2000,A,B ≤ 4 × 10^6Output对每组数据输出一行一个整数表示答案模2^30的值Sample Input5

2017-12-30 15:03:44 445

原创 bzoj2694 Lcm（反演）

Description 对于任意的>1的n gcd(a, b)不是n^2的倍数也就是说gcd(a, b)没有一个因子的次数>=2Input一个正整数T表示数据组数接下来T行每行两个正整数表示N、MOutputT行每行一个整数表示第i组数据的结果Sample Input4 2 4 3 3 6 5 8 3Sample Output24 28 233 178HINTT <=

2017-12-30 11:30:50 681

原创 bzoj2956 模积和（分块（平方的前缀和公式））

Description求∑∑((n mod i)*(m mod j))其中1<=i<=n,1<=j<=m,i≠j。Input第一行两个数n,m。Output一个整数表示答案mod 19940417的值Sample Input3 4Sample Output1样例说明　　答案为(3 mod 1) * (4 mod 2)+(3 mod 1) * (4 mod 3)+(3 mod 1) * (4 mo

2017-12-30 08:50:21 766

原创 luoguP3327 [SDOI2015]约数个数和（反演（关于d的结论））

题目链接分析：有一个非常显而易见的d的算法一开始我从这个式子下手，但是大家一眼就可以看出这是没有前途的实际上，这道题有一个很鬼的结论，在这里我们先给出：简单的证明以上是曲神给出的证明（不过说实话我是看不大懂的）所以以下用另一种方式证明一下：我们在这里先给出一个简单点的式子：显然，这个式子的前缀和就是上面的结论现在我们就来证明这个式子是对的：我

2017-12-30 07:02:35 300

原创 bzoj2226 [Spoj 5971] LCMSum（反演（互质数之和的公式））

DescriptionGiven n, calculate the sum LCM(1,n) + LCM(2,n) + .. + LCM(n,n), where LCM(i,n) denotes the Least Common Multiple of the integers i and n.InputThe first line contains T the number of test ca

2017-12-29 18:34:47 382

原创 bzoj2693 jzptab（反演LCM之和）

DescriptionInput一个正整数T表示数据组数接下来T行每行两个正整数表示N、MOutputT行每行一个整数表示第i组数据的结果Sample Input1 4 5Sample Output122HINTT <= 10000 N, M<=10000000分析：关于式子的反演：LCM之和（我实在不想再来一遍，不过听说舒老师有一种好方法，只能等ta传授了）所以我们就有了这

2017-12-29 15:26:54 396

原创 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 Input2 10 10 100 100 S

2017-12-29 14:32:39 400

原创 bzoj2154 Crash的数字表格（反演）

Description今天的数学课上，Crash小朋友学习了最小公倍数(Least Common Multiple)。对于两个正整数a和b，LCM(a, b)表示能同时被a和b整除的最小正整数。例如，LCM(6, 8) = 24。回到家后，Crash还在想着课上学的东西，为了研究最小公倍数，他画了一张N*M的表格。每个格子里写了一个数字，其中第i行第j列的那个格子里写着数为LCM(i, j)。一个

2017-12-28 19:22:34 305

原创 bzoj1257 [CQOI2007]余数之和sum（关于反演sqrt(n)复杂度的简单解释）

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,

2017-12-28 17:17:43 352

原创 bzoj2818 Gcd（反演）

Description给定整数N，求1<=x,y<=N且Gcd(x,y)为素数的数对(x,y)有多少对.Input一个整数NOutput如题Sample Input 4Sample Output 4HINT对于样例(2,2),(2,4),(3,3),(4,2)1<=N<=10^7 [Submit][Status][Discuss]分析：题

2017-12-28 14:59:50 349

原创 bzoj2705 [SDOI2012]Longge的问题（反演）

DescriptionLongge的数学成绩非常好，并且他非常乐于挑战高难度的数学问题。现在问题来了：给定一个整数N，你需要求出∑gcd(i, N)(1<=i <=N)。Input一个整数，为N。Output一个整数，为所求的答案。Sample Input6Sample Output15HINT【数据范围】对于6...

2017-12-27 20:11:28 233

原创 hdu1695 GCD（反演）

Problem DescriptionGiven 5 integers: a, b, c, d, k, you’re to find x in a…b, y in c…d that GCD(x, y) = k. GCD(x, y) means the greatest common divisor of x and y. Since the number of choices may be ver

2017-12-27 19:28:35 606

原创 hdu2588 GCD（反演）

Problem Description The greatest common divisor GCD(a,b) of two positive integers a and b,sometimes written (a,b),is the largest divisor common to a and b,For example,(1,2)=1,(12,18)=6. (a,b) can be

2017-12-27 17:11:40 245

原创 bzoj2879 [NOI2012]美食节（费用流）

题目链接分析：看完题面，我就想到了这道题实际上，这道题的建图就是这样：写完之后 —> 60 （一定要算准空间）边的数量比较多，所以如果一口气把所有的边都连上的话，时间是承受不了的由于我们跑一次spfa只能找出一次增广路，所以我们可以暂时不连不需要的边一开始，我们把所有厨师做倒数第１道菜与所有菜连好，然后找一条增广路，这条增广路上一定经过了一个点，表示第j个厨师做倒数第１道菜，

2017-12-26 19:10:39 273

原创 bzoj1901 Zju2112 Dynamic Rankings（可修改主席树）

Description给定一个含有n个数的序列a[1],a[2],a[3]……a[n]，程序必须回答这样的询问：对于给定的i,j,k，在a[i],a[i+1],a[i+2]……a[j]中第k小的数是多少(1≤k≤j-i+1)，并且，你可以改变一些a[i]的值，改变后，程序还能针对改变后的a继续回答上面的问题。Input第一行有两个正整数n(1≤n≤10000)，m(1≤m≤10000)。分

2017-12-26 17:40:26 611

原创可做题1（Code+第二次网络赛）

题目背景“codeplus比赛的时候在做什么？有没有空？能来解决停机问题吗？”qmqmqm这样问sublinekelzrip。当然，sublinekelzrip并不会停机问题，所以qmqmqm改为提出了另一个题目，现在请你帮助sublinekelzrip解决这个题目。题目描述这个问题是这样的：对于任何一个nn阶方阵，若任意从其中选择nn个不同行不同列的位置，其上的权值之和均相等，则我们称这个矩阵是巧

2017-12-25 20:14:31 616

原创可做题2（Code+12月网络赛）（扩欧+矩阵快速幂）

题目背景“codeplus比赛的时候在做什么？有没有空？能来解决丢番图方程问题吗？”sublinekelzrip这样问qmqmqm。当然，qmqmqm并不会丢番图方程问题，所以sublinekelzrip改为提出了另一个题目，现在请你帮助qmqmqm解决这个题目。题目描述这个问题是这样的：若一个数列aa满足条件an=an−1+an−2,n>3,而a1,a2a_n=a_{n-1}+a_{n-2},n>

2017-12-25 19:14:57 790

原创年前日常一览

之前有人说我的blog有点水，所以就把日常合到一起了这段日常的起源，就是平安夜的爆零赛12.24 爆零的一天，挺背，但是晚上收到了霖和黄的贺卡，心情好了不少12.25 早上卡点进教室，于是这件事耗完了一天的欧气上午上课一直被老师找茬，题也做不对，课也听得迷迷糊糊，中午赶工复习一下所学下午给霖写回信，写题效率不是很高感觉会考了之后，就找不到学竞赛的感觉了在曲神，ZYXZYXZYX和舒老

2017-12-25 18:25:15 266

原创可修改主席树

前言：这个就是传说中的树套树。。。主席树首先我们先简单的复习一下主席树（我的主席树入门）：主席数实际上就是一种权值线段树 ta能解决的一个经典问题就是区间第K大（后来知道也可以用整体二分解决）给定一个长度为n的序列，询问区间第K大我们建了n棵权值线段树，每一棵线段树维护序列1~i的信息查询的时候，计算区间元素个数t，与k比较，之后继续向左右结点查找建树过程：可修改给定

2017-12-25 17:45:30 1404

原创爆零赛T2（矩阵加速递推式）

前言：曲神：今天中午做了个噩梦，梦见给DP讲题，讲了半个小时。。。分析：终于在曲神，ZYXZYXZYX和舒老师的共同努力下，zz博主终于学会了70%做法动态规划我们首先就要设计状态：f[i][j]：表示到第i位为止，[i-j+1,i]这j位每一位上的颜色都不同（一共有j种颜色），这种情况下的方案数我们需要从f[i-1]转移到f[i] 有两种情况：（色差盲福利）得到的状态转移：

2017-12-25 15:31:41 376

原创复数（i？）

前言：之所以写这篇blog，是因为今天发生了一件足以使我骂街的事情平安夜，虽然离天黑还有5个小时，但是空气中已经弥漫着圣诞节姜饼的味道，人们不慌不忙的在晴好的阳光下漫步（往教室走去），每个人脸上洋溢的都是对节日的期待和对未来的幻想金黄的阳光透过窗户，洒到机房雪白的墙上，也映到微微发亮的电脑屏幕上原本不是很暖和的机房，本只能靠欢声笑语提升一点气氛，然而此时此刻却出奇的沉寂 loli的

2017-12-24 18:42:17 2380

原创爆零赛（12.24）

平安夜的hu测 mdzz最近就是有点颓，明知道这次成绩要上报，还是妥妥的爆了零。。。T1 分析：%60：做这道题之前，要了解一些复数运算手玩了r的1~8次方之后，发现i前面一定会有一个sqrt(7)，而常数项中一定没有sqrt(7) 实际上这就是复数的计算法则：而我们需要求的答案也是符合这个形式的（只有i前面有sqrt(7)）那我们就可以把i前面的sqrt(7)抹掉（但是如果

2017-12-24 18:25:09 329

原创 hdu1506 Largest Rectangle in a Histogram（单调栈）

Problem Description A histogram is a polygon composed of a sequence of rectangles aligned at a common base line. The rectangles have equal widths but may have different heights. For example, the figure

2017-12-24 11:14:35 320

原创 bzoj1113 [Poi2008]海报PLA（单调栈）

DescriptionN个矩形,排成一排. 现在希望用尽量少的矩形海报Cover住它们.Input第一行给出数字N,代表有N个矩形.N在[1,250000] 下面N行,每行给出矩形的长与宽.其值在[1,1000000000]2 1/2 PosteringOutput最少数量的海报数.Sample Input5 1 2 1 3 2 2 2 5 1 4 Sample Output4 分

2017-12-24 10:21:53 311

原创 poj2796 Feel Good（单调栈）

DescriptionBill is developing a new mathematical theory for human emotions. His recent investigations are dedicated to studying how good or bad days influent people’s memories about some period of life

2017-12-24 08:42:06 255

原创 poj3250 Bad Hair Day（单调栈）

题目链接分析：题目描述很迷。。。 bad hair day 英[bæd hɛə dei] 美[bæd hɛr de] n. “头发糟糕的日子”，暗示一整天都不顺利，心情不好。美国相当流行的一个新口语 [例句]All this fuss is because Carol is having a bad hair day. 就因为卡萝尔今天不太顺心，所以弄得大家鸡飞狗跳

2017-12-23 21:01:51 249

原创单调栈（单调队列的孪生兄弟？）

°Faust、 Mario_sz从题目就可以看出，单调栈和单调队列很相似首先我们给出简单介绍：简介定义：单调栈就是栈内元素单调递增或者单调递减的栈实际上，如果我们用一些生活情境来说明会更好（看到的一个很形象的例子）ZYXZYXZYX约好和zhhe0101一起下馆子，身为小受，ZYXZYXZYX当然要提前去定桌子啦但是那个饭店实在是太火了，大清早到了那里就已经排起了长长的队，

2017-12-23 19:20:19 343

原创 bzoj3879 SvT（后缀数组+单调栈）

Description(我并不想告诉你题目名字是什么鬼)有一个长度为n的仅包含小写字母的字符串S,下标范围为[1,n].现在有若干组询问,对于每一个询问,我们给出若干个后缀(以其在S中出现的起始位置来表示),求这些后缀两两之间的LCP(LongestCommonPrefix)的长度之和.一对后缀之间的LCP长度仅统计一遍.Input第一行两个正整数n,m,分别表示S的长度以及询问的次数.接下来一行有

2017-12-23 18:56:34 376

原创 poj1696 Space Ant（TuB）

DescriptionThe most exciting space discovery occurred at the end of the 20th century. In 1999, scientists traced down an ant-like creature in the planet Y1999 and called it M11. It has only one eye on...

2017-12-23 09:10:28 251

原创组合数学之四 —— 鸽巢原理

嘿，大家统统围过来，给你们看一个很棒的东西，那是一种很酷很辣很炫很Top的东西：基本知识如果要把n+1个物品放进n个盒子中，那么至少有一个盒子包含两个或更多东西这个就是鸽巢原理，本文完原理很简单，我们提出两个推论：如果将n个物体放入n个盒子并且没有一个盒子是空的，那么每个盒子恰好有一个物体如果将n个物体放入n个盒子中并且没有一个盒子被放入多于一个的物体，那么每个盒子里有一个物体看一下例

2017-12-22 19:45:01 2359

原创 hdu1398 Square Coins（生成函数）

Problem Description People in Silverland use square coins. Not only they have square shapes but also their values are square numbers. Coins with values of all square numbers up to 289 (=17^2), i.e., 1-

2017-12-22 17:41:42 198

原创 hdu1059 Dividing（生成函数||背包）

Problem Description Marsha and Bill own a collection of marbles. They want to split the collection among themselves so that both receive an equal share of the marbles. This would be easy if all the mar

2017-12-22 15:41:43 315

原创 poj2084 Game of Connections（Catalan）

DescriptionThis is a small but ancient game. You are supposed to write down the numbers 1, 2, 3, … , 2n - 1, 2n consecutively in clockwise order on the ground to form a circle, and then, to draw some s

2017-12-22 15:30:14 296

原创 poj3734 Blocks（指数生成函数）

DescriptionPanda has received an assignment of painting a line of blocks. Since Panda is such an intelligent boy, he starts to think of a math problem of painting. Suppose there are N blocks in a line

2017-12-21 20:41:09 939

原创 hdu1171 Big Event in HDU（生成函数）

Problem Description Nowadays, we all know that Computer College is the biggest department in HDU. But, maybe you don’t know that Computer College had ever been split into Computer College and Software

2017-12-21 18:29:44 329

原创 bzoj3028 食物（生成函数+Lucas）

Description明明这次又要出去旅游了，和上次不同的是，他这次要去宇宙探险！我们暂且不讨论他有多么NC，他又幻想了他应该带一些什么东西。理所当然的，你当然要帮他计算携带N件物品的方案数。他这次又准备带一些受欢迎的食物，如：蜜桃多啦，鸡块啦，承德汉堡等等当然，他又有一些稀奇古怪的限制：每种食物的限制如下：承德汉堡：偶数个可乐：0个或1个鸡腿：0个，1个或2个蜜桃多：奇数个鸡块：4的倍...

2017-12-20 20:38:56 1055

原创组合数学之三 —— 生成函数

前言：我觉得生成函数是一个很玄妙的东西知识铺垫我们在这里简单的介绍一下几种为人熟知的数列（给出公式及部分证明）：算术数列，其中的每一项比前一项大一个常数q（等差数列）几何数列，其中的每一项是前一项的常数q倍（等比数列）一 . 算术数列An=A0+nq Sn=( n-1 ) * A0+( q * n *（n+1）)/2二 . 几何数列An=A...

2017-12-20 19:02:14 16345 11

原创组合数学之二 —— 容斥原理及应用

这会是一个大型连载blog 内容主要来自《组合数学》+博主的思考容斥原理之前在“组合数学一”中提到过容斥原理我们在这里直接给出推论：设Ai表示在集合S中拥有特征Pi的元素子集，则集合S中至少具有性质P1,P2,P3,…,Pm之一的对象个数由下式表示：在这一节中，有一个经典问题：∞ 例一确定多重集合T={ 3*a，4*b，5*c }的10组合的数目这道题我们要充分利用容斥原理：

2017-12-20 15:48:15 3515 1

原创 bzoj2738 矩阵乘法（整体二分+二维树状数组）

Description　　给你一个N*N的矩阵，不用算矩阵乘法，但是每次询问一个子矩形的第K小数。Input 　　第一行两个数N,Q，表示矩阵大小和询问组数；　　接下来N行N列一共N*N个数，表示这个矩阵；　　再接下来Q行每行5个数描述一个询问：x1,y1,x2,y2,k表示找到以(x1,y1)为左上角、以(x2,y2)为右下角的子矩形中的第K小数。Output　　对于每组询问输出第K

2017-12-19 21:13:29 361

空空如也

空空如也