ControlBear的博客_CSDN博客-数论,动态规划,数论---组合数学领域博主

原创 LeetCode 5.最长回文子串[区间DP]

题目：给定一个字符串 s，找到 s 中最长的回文子串。你可以假设s 的最大长度为 1000。示例 1：输入: "babad"输出: "bab"注意: "aba" 也是一个有效答案。示例 2：输入: "cbbd"输出: "bb"来源：力扣（LeetCode）链接：https://leetcode-cn.com/problems/longest-palindr...

2019-11-08 18:23:07 80

原创 HDU5663 Hillan and the girl[莫比乌斯反演]

U - Hillan and the girl HDU - 5663题意:给T组数据，每组数据有一个n和一个m，求，其中题解:我们反过来，先求出 gcd(i,j) 是平方数的组合个数，因为感觉这个更好算，然后用 n*m 减去这些个数就是答案了。设 f(x) 为 gcd(i,j)=x 的个数F(x)

2017-09-04 20:26:10 575

原创 POJ3604 Professor Ben[线性筛/暴力]

Professor Ben POJ3604 题意:给 T 组数据，每组数据含有一个 n，让你求出(d(x)是指 x 的因子个数)题解:一看到题目，我们会比较倾向于这样的打表。void initial(){ sum[1]=1; for (int i=2 ; i<N ; ++i) {

2017-09-01 16:14:22 346

原创 LightOJ1052 String Growth[矩阵快速幂]

F - Problem F LightOJ1052 题意:一串字符串，它只有 a 和 b 组成，假设第 i 个字符串是 abab 那么第 i+1个字符串则为 b(ab)b(ab) 即下一个字符串，是由上一个字符串，通过将 a 变为 b，将 b 变为 ab，得到的。给你第 N 个字符串的跟第 M 个字符串的长度，求出第 K 个字符

2017-09-01 10:12:57 261

原创线性筛 [约数个数/约数和]

刚才手动推了一下用线性筛筛约数个数和约数和，就顺便写篇博客记录一下。不过网上应该也有不少人推过了。

2017-08-28 19:52:58 4929 4

原创 SPOJ VLATTICE Visible Lattice Points[莫比乌斯反演]

A - Visible Lattice Points SPOJ - VLATTICE 题意:题解:#pragma comment(linker, "/STACK:102400000,102400000")#include#include#include#include#include#include#includ

2017-08-21 16:49:55 268

原创 BZOJ2818 Gcd[莫比乌斯反演]

E - Gcd HYSBZ2818BZOJ2818题解:首先根据题意，设 f(i) 为 gcd(x,y)=i 的对数。对应的设( d=k*j [k>=1] 因为总是忘记整除左大还是右大)F(j)我们可以很容易求出来，就是，因为 F(j) 代表在 n 里面所有 gcd(x,y)=i 其中 i 是 j 的倍数的所有情况。

2017-08-21 11:18:04 458

原创 HDU6156 Palindrome Function[数位DP]

Palindrome Function HDU6156题意:给T组数据，每组数据含有 L，R，l，r，求出函数其中f(i,j)表示题解:#pragma comment(linker, "/STACK:102400000,102400000")#include#include

2017-08-20 09:08:51 271

原创 HDU6129 Just do it[组合数学]

Just do it HDU - 6129 题意:给出T组数据，每组含有一个n和m，n代表有n个数，m代表进行m次操作。接下来给出这n个数的值a[i]。操作是b[i]=a[1]^a[2]^...^a[i]，求出m次操作后b[1]到b[n]的值，每次操作后，用b[i]覆盖a[i]。题解:我们以 n=4 a={4,6,5,7

2017-08-16 21:10:15 259

原创 BZOJ3884 上帝与集合的正确用法[指数循环节]

上帝与集合的正确用法 HYSBZ - 3884题解:根据指数循环节的式子，我们可以对指数不断的进行递归取模，因为phi(2)==1，所以最后必然能递归到指数为0的情况，所以必然有结果。而且因为是不断的进行平方，所以必然符合B>=phi(C)这个条件，可以使用这个式子来计算。#pragma c

2017-08-16 16:57:27 271

原创 HDU6127 Hard challenge[计算几何]

Hard challenge HDU - 6127 题意:给T组数据，每组数据有一个n，代表n个点，接下来n行，每行给出一个x和y还有val，代表每个点的位置和值。两个点连起来的值等于两点的val相乘，不存在两点同时存在于一条与原点相连的线上。问从原点处划一条直线，不与这n个点相交(即点都不能存在于这条直线上)，能经过最大的边和是多少。

2017-08-16 15:18:08 320

原创 HDU6092 Rikka with Subset[母函数]

Rikka with Subset HDU - 6092 题意:给出一个T，代表T组数据。接下来每组数据含有一个n和m，表示这个序列有n个数，这n个数的和是m。再接下来给出m+1个数，表示这n个数任意组合(不重复组合)得到的和的方案数。根据所给的n，m，Bi求出原序列A的值。如果有多组答案，输出字典序最小的一组。，

2017-08-09 09:57:46 257

原创 HDU6061 RXD and functions[NTT]

HDU6061 RXD and functions

2017-08-02 17:16:23 489

原创 HDU6030 Happy Necklace[矩阵快速幂]

Happy Necklace HDU - 6030 题意:小Q要弄一条项链给女朋友，而弄这条项链的要求，就是要每素数长度里，红色的珠子的个数都不能少于蓝色珠子的个数(就是说红色个数小于等于蓝色个数)，然后给你n个珠子(红色蓝色都可以)，问弄成符合这样条件的一条项链，能有多少种方案。题解:要保证每个素数长度都符合这个条件

2017-07-30 21:13:12 285

原创 HDU6050 Funny Function[矩阵快速幂]

F - Funny Function 题意:给出一个函数由题目给的公式可以求出其他项，给出一个n和m，根据这个求出。结果取模1e9+7题解:根据前几项1 1 3 5 11 21 43 85 171 341 683 ......

2017-07-28 11:10:05 346

原创 ZOJ3626 Treasure Hunt Ⅰ[树形DP]

N - Treasure Hunt I ZOJ - 3626 题意:CC在一个有n个城市的地方住着，而他住在K这座城市，他有m天的时间出去外面收集宝物并且回到自己的城市，如果超过m天回不来，他就会被杀死。每座城市都有相连的城市，n座城市只有n-1条路。题解:树形DP+01背包因为必须返回

2017-07-18 09:56:04 400

原创 HDU1520 Anniversary party[树形DP]

A - Anniversary party HDU - 1520 题意:给一个n，代表有n个人，接下里n行给出一个第i个人的欢乐值，再接下来n-1行给出两两之间的关系(有直接联系，就是表示一个是上司一个是下属)。要求不能跟其有直接联系的人参加这个聚会(上司跟下属不能同时参加这个聚会)，问这个聚会最大的欢乐值是多少。题解:

2017-07-18 09:34:42 332

原创 Codeforces327C Magic Five[组合数学]

C. Magic Fivetime limit per test1 secondmemory limit per test256 megabytesinputstandard inputoutputstandard outputThere is a long plate s containing n digits.

2017-07-14 15:47:30 315

原创 HDU2152 Fruit[母函数/背包]

B - Fruit HDU - 2152

2017-07-14 15:21:00 280

原创 ZOJ3856 Goldbach[FFT]

GoldbachZOJ3856

2017-07-08 15:46:01 565

原创 Codeforces229A Shifts[二分]

A. Shiftstime limit per test2 secondsmemory limit per test256 megabytesinputstandard inputoutputstandard outputYou are given a table consisting of n rows and m

2017-07-07 11:33:47 351

原创 HDU1402 A * B Problem Plus[FFT]

B - A * B Problem Plus HDU - 1402 题意:计算A*B(A,B均不超过5e4位)。题解:计算A*B 可以把AB的各项化成多项式来相乘把这两个多项式相乘起来就能得到A*B的结果了。用FFT可以快速的求出两多项式相乘的结果，然后就可以得到每一位的的结果。求得A*B的结果之后，做进

2017-07-03 16:15:59 279

原创 UVALive6886 Golf Bot[FFT]

M - Golf Bot UVALive - 6886 Do you like golf? I hate it. I hate golf so much that I decided tobuild the ultimate golf robot, a robot that will never miss a shot. I simplyplace it over the ball

2017-07-03 15:46:15 406

原创 HDU5456 Matches Puzzle Game[数位DP]

E - Matches Puzzle Game HDU - 5456

2017-07-01 20:24:36 314

原创 HDU4734 F(x)[数位DP]

D - F(x) HDU - 4734 题意:求出存在[0,B]中x使得F(x)小于等于F(A)的数有多少个。题解:首先求出F(A)的值，用这个值作为标记的初始值，第一维代表当前位置，第二维代表剩余值。而第二维的取值，最多能达到9*9*2^8，计算过程中，如果sum小于0，代表不可能因为要保证必须比F(A)小，直接返回0，而减的过程中，直接用i*(1&lt;&lt;pos) (pos为当前位置，i...

2017-07-01 10:17:20 275

原创 HDU3709 Balanced Number[数位DP]

C - Balanced Number HDU - 3709 题意:给T组数据，每组数据给出一个x和y，求出[x,y]当中满足以其中一位为平衡点，满足在平衡点的左右两边的力矩大小一致。题解:这道题要注意前导零的问题，因为不记录前导零，如果左界也不是0，会多出len-1种情况(len是这个数共有多少位)，那么如果左界是0

2017-07-01 09:47:43 328

原创 HDU3555 Bomb[数位DP]

B - Bomb HDU - 3555 题解:只要记录下三个状态，0号状态代表前面没有4，也没有49，1号状态代表前一位是4，那么只要当前位是9，就可以直接把状态转成2号状态，即当前的数已经存在49后面的数可以任意枚举。#pragma comment(linker, "/STACK:102400000,1024

2017-07-01 09:00:16 229

原创 HDU2089 不要62[数位DP]

A - 不要62 HDU - 2089 题解:统计[n,m]中不含4 或者 62的数有多少个,最基础的数位DP题目,首先判断有4就直接continue,然后判断如果前一个是6,下一个是2也continue,剩下的都是可行的.#include#include#include#include#include#include#include#

2017-06-30 23:42:31 230

原创 HDU5446 Unknown Treasure[中国剩余定理+Lucas]

J - Unknown Treasure HDU - 5446 题意:T组数据, 给你n,m,k ,接下来给你k个素数求C(n,m)%(这k个素数的积)题解:题目给我的M 是由(p1*p2*...*pk)组成的,给我的第一反应是中国剩余定理.那么照着这个方向,我们就必须得得到同余式.我们的同余式则是又 x%pi=ai 那么我们的x

2017-05-18 15:02:24 384

原创 ZOJ3203 Light Bulb[三分/推公式]

Light BulbTime Limit:1000MS Memory Limit:32768KBSubmit StatusDescriptionCompared to wildleopard's wealthiness, his brother mildleopard is rather poor. His house is na

2017-01-08 22:38:09 2413 8

原创 POJ1064 Cable master[二分]

Cable master POJ - 1064 Time limit1000 msMemory limit10000 kBOSLinuxSourceNortheastern Europe 2001题意:给出绳子的个数N，再给出需要切割出K条绳子，接下来N行，给出每条绳子的长度，求切割出来的长

2016-12-30 12:45:08 340

原创 POJ1064 Cable master[二分]

Time limit1000 msMemory limit10000 kBOSLinuxSourceNortheastern Europe 2001Cable master POJ - 1064 题意:给N条绳子,将其分成等长的K条绳子,下面给出N行,每行代表第i条绳子的长度

2016-12-30 12:40:11 203

原创 HDU4407 Sum[容斥定理]

D - SumTime Limit:1000MS Memory Limit:32768KB 64bit IO Format:%I64d & %I64uSubmit Status Practice HDU 4407DescriptionXXX is puzzled with the question below: 1, 2, 3

2016-11-09 15:31:33 305

原创 LightOJ1028 1090 1138 Trailing Zeroes (I)(Ⅱ)(Ⅲ)[算数基本定理/打表/打表+二分]

C - Trailing Zeroes (I)Time Limit:2000MS Memory Limit:32768KB 64bit IO Format:%lld & %lluSubmit StatusDescriptionWe know what a base of a number is and what the prop

2016-11-09 15:09:43 621

原创 HDU5768 Lucky7[中国剩余定理+容斥定理]

E - Lucky7Time Limit:1000MS Memory Limit:65536KB 64bit IO Format:%I64d & %I64uSubmit StatusDescriptionWhen ?? was born, seven crows flew in and stopped beside him. In its c

2016-08-07 08:29:16 914

原创 HDU1796 How many integers can you find[容斥定理]

D - How many integers can you findTime Limit:5000MS Memory Limit:32768KB 64bit IO Format:%I64d & %I64uSubmit StatusDescription Now you get a number N, and a M-integers set

2016-07-31 14:57:58 263

原创 HDU5733 tetrahedron[计算几何]

K - tetrahedronTime Limit:1000MS Memory Limit:65536KB 64bit IO Format:%I64d & %I64uSubmit Status Practice HDU 5733DescriptionGiven four points ABCD, if ABCD is a tetrahed

2016-07-27 14:46:40 595

原创 HDU1166 敌兵布阵[线段树/树状数组]

X - 敌兵布阵Time Limit:1000MS Memory Limit:32768KB 64bit IO Format:%I64d & %I64uSubmit StatusDescriptionC国的死对头A国这段时间正在进行军事演习，所以C国间谍头子Derek和他手下Tidy又开始忙乎了。A国在海岸线沿直线布置了N个工兵营地,Dere

2016-07-24 19:12:39 283

原创 SGU106 The equation[扩展欧几里德算法]

A - The equationTime Limit:250MS Memory Limit:4096KB 64bit IO Format:%I64d & %I64uSubmit StatusDescriptionThere is an equation ax + by + c = 0. Given a,b,c,x1,x2,

2016-07-24 18:55:04 931

原创 HDU5301 Buildings [图形]

B - BuildingsTime Limit:2000MS Memory Limit:131072KB 64bit IO Format:%I64d & %I64uSubmit Status Practice HDU 5301DescriptionYour current task is to make a ground plan for

2016-07-24 16:41:27 495