- 博客(34)
- 资源 (4)
- 收藏
- 关注
RSS订阅原创 codeforces 682 C
题意给出一棵树，dist(u,v)>au的点为sad点，dis(u,v)为从u到子结点v上的权值之和，au为结点u上的权值。计算删除多少个叶子结点，使得该树不存在sad点。dfs 从根结点开始计算，假设dist>au 那么 以该结点作为父节点的子树都不符合条件。假设distimport java.util.ArrayList;import java.util.Lin
2016-06-19 16:29:15
632
原创 hdu 1154
给出两个点所在直线 然后求出直线和多边形的交线长度#include #include #include #include #include #include using namespace std;const double eps=1e-8;const int maxn=1e3+5;struct Point{ double x,y; Point(doub
2015-09-19 21:47:51
319
原创 hdu 4502
题意：就是给出一个矩阵,矩阵上有非负的数值。然后你需要从(1,1)->(n,m)把过程中的值加起来。然后求出路径上最大的路线值的和。/*********************************************************************************************************************************
2015-09-10 17:59:59
366
原创 uva 11889
给出a,c 求出最小的b .利用质因数分解的思想拆开a和c的因子，然后凑出b#include #include #include #include #include #include using namespace std;typedef long long ll;const int maxn=1e4+5;int prime[maxn],len=0;bool vis[m
2015-09-02 17:18:42
346
原创 poj 1845
题目:求出A^B的因子和modp的值思路:利用sigma函数计算即可。sum=(p1^(k1*B+1)-1)/(p1-1)*... 然后就是将/(p1-1)转化成逆元去计算 就可以避免/运算。 但是这边有一个问题是 如果p|p1-1 那么逆元无解。wa情况。 也就只能另外计算了。p1=k*p+1 因为对于乘法同余,所以值为(k1*B+1)#incl
2015-08-27 00:12:42
325
原创 4002 大数 java
import java.util.*;import java.io.*;import java.math.*;public class Main{static int[] pri=new int[1000];static int len=0,anslen;static BigInteger n,up,temp;static boolean[] vis=new boo
2015-08-26 11:00:54
313
原创 素数测试
费马小定理: p是质数,则 0 二次探测定理: 0一次的miller-rabin如果p是质数,则偶数 p-1=m*2^q。求出m和q。然后先计算x^m的值.然后计算temp= (x^m)*2*k%p的值temp=1 x！=1 和x！=p-1 ，p就不是素数 。不满足二次探测temp!=1 把取模后的值作为新的x 循环计算知道k=q的时
2015-08-18 10:17:10
388
原创 同余式和中国剩余定理
x≡b1(mod a1) 1x≡b2(mod a2) 2x≡b3(mod a3) 4x≡b4(mod a4) 5x=k1*a1+b1≡b2(mod a2)k1*a1≡b2-b1(mod a2)裴蜀定理：对于方程 a*x+b*y=n；有整数解的充分必要条件是（n % gcd（a，b）== 0）， 即n能够被a和b的最大公约数整除，其实
2015-08-17 10:25:41
710
原创 sgu 106
给出a,b,c,x1,x2,y1,y2然后你要做的就是找出在矩形框内的整数点有多少个。初始解(x*(-c)/g,y*(-c)/g)然后对于x轴方向上的增量是 b/g我写成几何了 我去 。。。好好的数论题。先用扩展欧几里德求出初始解,获取到在离矩形框最近的点的x坐标。 然后计算直线和x=x1,x2，y=y1和y2上的点。然后取中间两个。判断两个点是否都在矩形框上。如果存在
2015-08-16 01:30:19
336
原创 欧几里得
欧几里得gcd(a,b)=gcd(b,a%b)证明:1. d=gcd(a,b) => d|a, d|b 令r=a-kb 所以 d|r 也就是 r 和 b 存在d的公约数。 假设 d=gcd( b,a%b ) , d|b ,d|r ==> d|( a-kb ) d|a 所以最大公约数满足 2. d=gcd(a,b) , a=
2015-08-14 11:20:31
415
原创 hdu 2588
题目给你 一个n和m 然后你要去求gcd(x,n)>=m的个数保证1对于gcd(x,n)>=m可以看做是 gcd(x,n)=m,gcd(x,n)=m+1.....gcd(x,n)=n; m必为n的倍数也就是 gcd(x/m,n/m)=1 也就是n/m内和n/m互质的数的个数。只要计算范围内的因子值的欧拉函数即可#include #include #include #i
2015-08-12 10:09:38
443
原创 hdu 3501
给出一个n求出与n公约数大于1的数的和。对于一个a与n互质, 那么n-a与a也互质。不存在a与n互质,n-a=a的情况。也就是说a 和n-a总是成对存在 那么phi[n]是偶数。那么和n互质数的和就是 n*phi[n]/2#include #include #include #include #include #include using namespace std;t
2015-08-11 21:03:32
350
原创 LightOJ 1197
给出一个区间 让你求区间内数的素数数量。#include #include #include #include #include #include using namespace std;const int maxn=50000;int prime[maxn],primelen;bool vis[maxn*2+5];void makeprime(){ prime
2015-08-10 15:37:13
636
2
原创 UVA 11426
给出代码段int G=0; for(int i=1; i for(int j=i+1; j { G+=gcd(i,j); }printf("%d\n",G);求G的值#include #include #include #include #i
2015-08-09 13:02:59
301
原创 LightOj Mysterious Bacteria
给出x=b^p要让p尽量大。给出的x范围为负。质因数分解，然后计算幂指数间的最大公约数。如果是负数结果不能为偶数#include #include #include #include #include using namespace std;const int maxn=1e6+5;int len=0;long long prime[maxn];bool status[m
2015-08-07 14:11:11
299
原创 hdu1124
求N!后面0的个数。产生0来源是2和5 2的话 每隔一个就有一个5的话 比较难得。所以关键就在于n！里面5的个数#include #include #include #include #include #include #include using namespace std;const int maxn=1e9+1;typedef long long ll;
2015-08-06 15:45:05
373
原创 hdu1111
求 X = a0 + a1B + a2B2 + ...+ anBn 中a0,a1,....an的值 其中X,Bi是复数。在这边要注意复数没有取模运算符。只能利用秦九韶公式 X=((((an)*b+an-1)*b+an-2)*b)*b...)*b+a0如果非复数 那我们就可以使用 进制取模的方法 来写。然并没有。只能是在每一层枚举规定范围内的值，然后找到符合要求的值#in
2015-08-06 13:41:01
548
原创 hdu1014
seed(x+1) = [seed(x) + STEP] % MODSTEP=3 MOD=5 出现的循环串是 0, 3, 1, 4, 2.假设gcd(step,mod)=g ,令step=k*g,mod=m*g; 那么x=f*k*g-(f*k*g)/(m*g)*(m*g)=g*( f*k % m)*m ) {0序列是 {0,step ,2step,3step,4step,..
2015-08-04 17:02:05
432
原创 hdu2601
求出n=i*j+i+j的对数是的(0#include #include #include #include #include #include using namespace std;typedef long long ll;const int maxn=1e6;ll prime[maxn],len=0;bool status[maxn];void pri(){
2015-08-04 09:08:49
299
原创 hdu1999 可以筛法做
好失败，竟然用打表才过的。#include #include #include #include #include #include using namespace std;const int maxn=1e6+5;int len=0;int prime[maxn],a[maxn];bool status[maxn];bool s[maxn]={0,1,0,1,1,0,1
2015-08-01 23:31:24
410
原创 hdu1593
深深地为自己的智商悲哀。说不定以后遇到这种情况都不能机智逃走。。。。。。直接往不同方向跑，直接错。应该先让自己和敌人和圆心处于三点共线。在速度允许的情况下，自己离开圆心但是始终保持三点共线。在也就是角速度一致。v2/R=v1/r 即 r=v1*R/v2如果r增大 就不能继续维持三点共线了。这个时候就只能直线跑，这个时候就看时间的问题了。你的路程就会使(R-r)
2015-08-01 11:11:37
369
原创 hdu1722
模拟下切的版块的数量，会发现其实切得过程中会出现切痕重叠。去重就好了#include #include #include #include #include #include #include #include #include #include #include #include using namespace std;int gcd(int a,int b){
2015-08-01 09:47:42
349
原创 hdu1220
给出一个n的立方体。然后你去拆成体积为1的立方体。对于拆出来的任意两个立方体。两个立方体之间的交点可能有 0,1,2,4只要所有的减去4的就是答案。四的 一列由n-1个,每个面有 N^2,总共上左下三个面#include #include #include #include #include #include #include #include #incl
2015-07-31 11:18:42
527
原创 hdu1141
求最大的n使得满足n!bit可以达到256超过long long然而n!=n*...*1.两边取对数,就是log(n)+log(n-1)+...+log(1)然后暴力求解。#include #include #include #include #include #include using namespace std;const double eps=1e-5;
2015-07-30 15:33:25
375
原创 求多边形的质心hdu1115
给出点。然后求多边形的质心。要求多边形的质心,就得先要知道三角形的质心要怎么求。三角形的质心(A+B+C)/3.可以用向量来证明。然后我们还得知道质点系的质心公式:在质量均匀分布在质点上。如果质量均匀分布在顶点上，那么这个公式是正确的。但是我们要求的是平面。所以质量均匀分布在平面上，可以将平面分成N-2个三角形(多边形三角剖分)。每个三角形的质心就可以求出来。这样平
2015-07-30 01:02:22
3916
原创 hdu1018
让你求阶乘的位数。对于位数为x的数的范围是 10^(x-1)解得x=log10(n)+1;然后用对数的性质对n进行分解#include #include #include #include #include #include using namespace std;int main(){ int t,n; scanf("%d",&t);
2015-07-29 08:32:15
259
原创 java高精度
1002/****************Author:linhanCreate Time:2015-07-23 20:57:10File Name:Main.java ************************* */import java.math.*;import java.util.*;public class Main{ public static vo
2015-07-23 23:03:26
205
原创 Harmonic Number (II)
给出一段求调和级数的值。数据量是int 直接循环爆炸。可能就是找规律了。以往都是靠答案找规律。都很少在计算的过程里面找规律。也算是学习了。针对27这个案例i:1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16....j:27 13 9 6 5 4 3 3 3 2 2 2 2 1 1 1
2015-07-22 14:17:49
356
原创 Goldbach`s Conjecture
给出哥德巴赫猜想:任意个偶数都可以表示成两个素数相加。数据范围是10^7开数组存素数表(傻逼的我1e6和664600大小分不清楚。。。。。然后对于判断是不是素数用bool就好了 不需要int型#include #include #include #include #include #include using namespace std;const int maxn=
2015-07-21 20:50:40
266
原创 Leading and Trailing
computer n^k leading three digits.n^k=a.lg(a)=k*lg(a).k*lg(a)=integer+decimal10^lg(a)=a=10^(k*lg(a))=10^(integer+decimal)=10^integer*10^decimal10^integer*10 算出的整十的倍数10^decimal才是计算真正结果的数字
2015-07-21 14:23:30
292
原创 Sigma Function
题目意思:给出sigma函数的意义：sigmafunction(x)=x的因子之和x=p1^(e1)*p2^e2*...*pk^eksigmafunction(x)=(p1^(e1+1)-1)/(p1-1)*(p2^(e2+1)-1)/(p2-1)*...*(pk^(ek+1)-1)/(pk-1).求1~n之间的偶数值的个数。对于指数为0的项,(px-1)/(px-1)=1
2015-07-20 17:37:59
823
原创 Aladdin and the Flying Carpet
题目意思:给你一个整数a和整数b 然后你需要求出存在多少个c*d=a (c!=d) 的组合数量 输出个数就好了。如果直接循环遍历的话会超时。然后就是利用唯一分解定理。筛素数表之前想说找出大于b的最小a因子 感觉除了搜索的话，在分解的时候实现麻烦。然后就是 算出总的因子数/2，然后减去小于b的。#include #include #include #inclu
2015-07-20 15:19:27
407
Scala函数式编程.pdf
2019-05-18
Scala编程 第3版.pdf
2019-05-16
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人