curry的博客

类欧, 万欧

..

..

筛法

前言:某天,作者想出一道 数论+斜率优化 的毒瘤题,但是他失败了,可是得出了一些漂亮的结论.

1.一句话定义：狄利克雷卷积是形如f(x)=g(x)∗h(x)f(x)=g(x)*h(x)f(x)=g(x)∗h(x)的式子（*是狄利克雷卷积符号，我也不知道手写叉乘行不行 ），它表示的含义为f(x)=∑d∣xg(d)∗h(x/d)f(x)=\sum_{d|x} g(d)*h(x/d)f(x)=∑d∣x​g(d)∗h(x/d).2.性质：狄利克雷卷积满足交换律，结合律。如：$f(x)*g(...

这是本人第一篇博客，写的不好请尽力吐槽。Time Limit: 5 Sec Memory Limit: 128 MBSubmit: 6612 Solved: 3190[Submit][Status][Discuss]Description给出正整数n和k，计算j(n, k)=k mod 1 + k mod 2 + k mod 3 + … + k mod n的值其中k mod i表示...

描述Hanks博士是BT（Bio-Tech，生物技术）领域的知名专家，他的儿子名叫Hankson。现在，刚刚放学回家的Hankson正在思考一个有趣的问题。今天在课堂上，老师讲解了如何求两个正整数c1和c2的最大公约数和最小公倍数。现在Hankson认为自己已经熟练地掌握了这些知识，他开始思考一个“求公约数”和“求公倍数”之类问题的“逆问题”，这个问题是这样的：已知正整数a0,a1,b0,b1...

给定整数n,g,求$q^{\sum_{k|n} C_n^k }\mod\ 999911659.999911659=234679*35617欧拉定理，中国剩余定理，Lucas定理杂糅.代码很友好：#include<cstdio>#include<cstring>#include<algorithm>using namespace std;typed...

Problem DescriptionIn our daily life we often use 233 to express our feelings. Actually, we may say 2333, 23333, or 233333 … in the same meaning. And here is the question: Suppose we have a matrix ca...

简单数论

lucas定理：当m∈[1,n]m\in [1,n]m∈[1,n]时，Cnm≡C⌊ n/p⌋⌊ m/p⌋∗Cn mod pm mod pC_n^m≡C_{\left\lfloor\ n/p \right\rfloor}^{\left\lfloor\ m/p \right\rfloor}*C_{n ~ mod ~ p}^{m ~ mod ~...

给定一个 1~n 的排列 p1,p2,…,pnp_1,p_2,…,p_np1​,p2​,…,pn​，可进行若干次操作，每次选择两个整数 x,y，交换 px,py。设把p1,p2,…,pnp_x,p_y。设把 p_1,p_2,…,p_npx​,py​。设把p1​,p2​,…,pn​ 变成单调递增的排列 1,2,…,n 至少需要 m 次交换。求有多少种操作方法可以只用 m 次交换达到上述目标。因为结果...

描述在一条数轴上有N家商店，它们的坐标分别为 A[1]~A[N]。现在需要在数轴上建立一家货仓，每天清晨，从货仓到每家商店都要运送一车商品。为了提高效率，求把货仓建在何处，可以使得货仓到每家商店的距离之和最小。输入格式第一行一个整数N，第二行N个整数A[1]~A[N]。输出格式一个整数，表示距离之和的最小值。样例输入46 2 9 1样例输出12数据范围与约定对于100%的数...

Description　　FGD正在破解一段密码，他需要回答很多类似的问题：对于给定的整数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。（...

描述石头游戏在一个 n 行 m 列 (1≤n,m≤8) 的网格上进行，每个格子对应一种操作序列，操作序列至多有10种，分别用0~9这10个数字指明。操作序列是一个长度不超过6且循环执行、每秒执行一个字符的字符串。每秒钟，所有格子同时执行各自操作序列里的下一个字符。序列中的每个字符是以下格式之一：数字09：表示拿09个石头到该格子。NWSE：表示把这个格子内所有的石头推到相邻的格子，N表示上...

Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 7745 Solved: 4295[Submit][Status][Discuss]Description　　有一个球形空间产生器能够在n维空间中产生一个坚硬的球体。现在，你被困在了这个n维球体中，你只知道球面上n+1个点的坐标，你需要以最快的速度确定这个n维球体的球心坐标，以便于摧毁这个球形...

Time Limit: 1000MS Memory Limit: 30000KDescription有N个相同的开关，每个开关都与某些开关有着联系，每当你打开或者关闭某个开关的时候，其他的与此开关相关联的开关也会相应地发生变化，即这些相联系的开关的状态如果原来为开就变为关，如果为关就变为开。你的目标是经过若干次开关操作后使得最后N个开关达到一个特定的状态。对于任意一个开关，最多只能进行一次开关...

原文矩阵的秩的定义：存在K阶子式不为0，对任意K+1阶子式均为0，则k即为矩阵的秩。向量组的秩的定义：向量组的极大线性无关组所包含向量的个数，称为向量组的秩。其次再弄清楚3个定理：1，矩阵A的行列式不为0的充要条件是A的行(列)向量线性无关2，无关组加分量仍无关3， r个n维列向量组线性无关的充要条件是这r个n维列向量组所构成的矩阵至少存在一个r阶子式不为0好了，简略证明过程开始,我...

数论中因数个数函数或者说是d（n）是个非常有趣的函数。它的意思是某个数正因数的个数。比如d（24）=8因为24有8个正因数：1,2,3,4,6,8,12和24。而数学中一个正整数n的阶乘表示为n!并且定义如下n!=123…n 另一个有趣的函数AF（n）（Again Factorial的缩写）定义如下AF（n）= 1!2!3!..*n! 给你一个n，你要做的就是算出d（AF（n））的值。输入格...

Description脸哥最近在玩一款神奇的游戏，这个游戏里有 n 件装备，每件装备有 m 个属性，用向量zi(aj ,…,am) 表示 (1 <= i <= n; 1 <= j <= m)，每个装备需要花费 ci，现在脸哥想买一些装备，但是脸哥很穷，所以总是盘算着怎样才能花尽量少的钱买尽量多的装备。对于脸哥来说，如果一件装备的属性能用购买的其他装备组合出（也就是说脸哥可...

Description有一张 n×m 的数表，其第 i 行第 j 列（1 <= i <= n, 1 <= j <= m）的数值为能同时整除 i 和 j 的所有自然数之和。给定 a , 计算数表中不大于 a 的数之和。Input输入包含多组数据。输入的第一行一个整数Q表示测试点内的数据组数接下来Q行，每行三个整数n，m，a(|a| < =10^9)描述一组数...

因为不同，所以存在。（由异或得到的想法）Problem DescriptionXOR is a kind of bit operator, we define that as follow: for two binary base number A and B, let C=A XOR B, then for each bit of C, we can get its value by che...

描述给定一个多项式(ax+by)^k, 请求出多项式展开后xn*ym项的系数。输入格式共一行，包含5 个整数，分别为 a ，b ，k ，n ，m，每两个整数之间用一个空格隔开。输出格式输出共1 行，包含一个整数，表示所求的系数，这个系数可能很大，输出对10007 取模后的结果。样例输入1 1 3 1 2样例输出3数据范围与约定对于30% 的数据，有 0 ≤k ≤10 ；对于...

DescriptionConsider two natural numbers A and B. Let S be the sum of all natural divisors of A^B. Determine S modulo 9901 (the rest of the division of S by 9901).InputThe only line contains the two...