![](https://img-blog.csdnimg.cn/20201014180756928.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数学
文章平均质量分 73
BryantMark
这个作者很懒,什么都没留下…
展开
-
HDU 4609
实质上是两个多项式相乘注意去重 #include #include #include #include #include using namespace std; typedef long long LL; const double pi = acos(-1.0); const int maxn = 271111; const double eps = 1e-6; struct Com原创 2013-07-29 01:15:27 · 672 阅读 · 0 评论 -
UVA 12009 - Avaricious Maryanna
#include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long LL; typedef unsigned long long ULL; typedef vector原创 2013-05-28 23:32:39 · 571 阅读 · 0 评论 -
UVA 11426 - GCD - Extreme (II)
#include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long LL; typedef unsigned long long ULL; typedef vector原创 2013-05-29 13:35:09 · 425 阅读 · 0 评论 -
hdu 3369
求递推系数 i^k 变化到 (i+7)^k #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #原创 2013-05-15 18:54:59 · 957 阅读 · 0 评论 -
poj 3471 整除与剩余
给出一个多项式 f(x)=x^n+An-1*x^(n-1)+.....+A0 求其所有零点(包括)重根, 由韦达定理的扩展我们知道所有的整数零点的绝对值都是A0的约数,可以枚举约数,然后计算f(x)/(x-a) a为当前枚举的约数,判断是否整除,如果A0=0,则我们令f(x) ,中每一项均除以x,直到常数项不为0,然后对常数项的绝对值进行质因子分解然后求取所有因数 f(x)=(x-x1原创 2013-05-13 21:17:19 · 639 阅读 · 0 评论 -
worldfinals 2008 Always an integer
利用差分序列 p(n) 为K 次多项式 p(0)p(1) p(2) .........p(k) 均整除D 其一阶 差分序列 dp(0),dp(1)....dp(n) 均能整除D,类推可知 d^kp(0),d^kp(1).....可整除D,且差分序列最后 变为常数 #include #include #include #include #include #原创 2013-04-22 16:10:17 · 338 阅读 · 0 评论 -
POJ 3372 —— Candy Distribution
#include #include #include #include #include #include #include using namespace std; typedef vector VI; typedef vector VVI; typedef long long ll; vector::const_iterator it; #define L ch[rt][0] #define原创 2013-05-13 12:18:08 · 464 阅读 · 0 评论 -
UVA 10288 - Coupons
求 n/1+n/2+...n/n 用最简分数表示 /* ID: BryatnMark LANG: JAVA TASK: AC */ import java.util.*; import java.io.*; import java.math.*; import java.util.*; public class Main { public static void main(String ar原创 2013-05-13 12:48:51 · 579 阅读 · 0 评论 -
SGU 261 Discrete Roots (经典K次剩余)
首先求p的原根g, g满足g^(phi(n)/pi) mod n 均不等于1,令g^y=x,g^t=a 得到 g^yk=g^t(modp) (a=0特判) 此式 可转化为 k*y=t(mod(phi(p))); 即 k*y=t(mod(p-1)),接下来就是求g^t=a(mod(p)) ,用babystep 求离散对数的方法求出t ,然后解模线性方程 k*y=t(mod(p-1)), 得到所有原创 2013-05-10 14:05:45 · 682 阅读 · 0 评论 -
UVA 11542 - Square
#include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long LL; const int maxn = 511; bitsetvis; int p[maxn],size; void i原创 2013-05-26 22:22:26 · 701 阅读 · 0 评论 -
POJ 2103 Jackpot
/* ID: BryatnMark LANG: JAVA TASK: AC */ import java.util.*; import java.io.*; import java.math.*; import java.util.*; public class Main { public static void main(String args[]) throws Exception {原创 2013-05-22 23:47:00 · 606 阅读 · 0 评论 -
POJ System
/* ID: BryatnMark LANG: JAVA TASK: AC */ import java.util.*; import java.io.*; import java.math.*; import java.util.*; public class Main { public static void main(String args[]) throws Exception {原创 2013-05-31 08:29:34 · 421 阅读 · 0 评论 -
LightOJ 1342 - Aladdin and the Magical Sticks
#include #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long LL; typedef unsigned long long ULL; typed原创 2013-05-30 12:57:58 · 938 阅读 · 0 评论 -
HDU 4608
#include #include #include #include #include #include #include #include using namespace std; typedef vector VI; typedef pair PII; const int maxn = 111111; char s[maxn]; int arr[maxn]; int main原创 2013-07-29 01:12:24 · 432 阅读 · 0 评论 -
HDU 4602
EASY PROBLEM #include #include #include using namespace std; typedef long long LL; const int maxn = 111111; const int mod = (int) (1e9 + 7); char s[maxn]; int num[maxn]; LL POW(LL a, LL b) { L原创 2013-07-29 01:10:27 · 396 阅读 · 0 评论 -
HDU 4611
#include using namespace std; typedef long long LL; template T gcd(T a, T b) { return b ? gcd(b, a % b) : a; } template T lcm(T a, T b) { return a / gcd(a, b) * b; } LL cnt, n, a, b; LL cal(int x,原创 2013-07-25 21:01:55 · 520 阅读 · 0 评论 -
UVA 1426 - Discrete Square Roots
枚举N的因子 合并模线性方程 暴力找答案 /* * test.cpp * * Created on: 2013-7-3 * Author: Administrator */ #include #include #include #include #include #include #include using namespace std;原创 2013-07-21 15:09:25 · 590 阅读 · 0 评论 -
hdu 4447
多项式 计算 模拟 多项式乘法 和 多项式除法就行了,多项式除法可以参考高精度除法 #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long LL; const int原创 2013-06-29 06:10:39 · 840 阅读 · 0 评论 -
HDU 3486
#include #include const int N = 1111; char buf[N], op[N]; int lch[N], rch[N]; int h[N]; int build(char s[], int x, int y, int lch[], int rch[], char op[], int &nc) { int i, u, c1, c2; int p = 0;原创 2013-06-04 13:17:17 · 716 阅读 · 0 评论 -
HDU 4196 Remoteland
#include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long LL; typedef unsigned long long原创 2013-05-31 20:30:08 · 505 阅读 · 0 评论 -
light oj 1274 - Beating the Dataset
#include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long LL; typedef unsigned long long原创 2013-05-30 21:31:36 · 693 阅读 · 1 评论 -
URAL 1306. Sequence Median
#include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long LL; typedef unsigned long long原创 2013-05-30 16:50:39 · 472 阅读 · 0 评论 -
UVA 11246 K-Multiple Free set
/* ID: BryatnMark LANG: JAVA TASK: AC */ import java.util.*; import java.io.*; import java.math.*; import java.util.*; public class Main { public static void main(String args[]) throws E原创 2013-05-23 20:54:04 · 459 阅读 · 0 评论 -
UVA 10236 The Fibonacci Primes
/* ID: BryatnMark LANG: JAVA TASK: AC */ import java.util.*; import java.io.*; import java.math.*; import java.util.*; public class Main { public static void main(String args[]) throws Exception {原创 2013-05-23 19:55:17 · 494 阅读 · 0 评论 -
UVA 10162 - Last Digit
/* ID: BryatnMark LANG: JAVA TASK: AC */ import java.util.*; import java.io.*; import java.math.*; import java.util.*; public class Main { public static void main(String args[]) throws Exception {原创 2013-05-23 18:09:35 · 991 阅读 · 0 评论 -
HDU 2837 calculation
容易证明 命题 若A^X>M,则 X>=phi(M) 的反例只有 A=2,X=2,M=6 而此特例中 2^1%6=2 2^2%6=4 2^3%6= 2, 当X=1时已经开始循环,故此特例不必特判 #include #include #include #include #include #include #include #include #include #include原创 2013-05-16 17:55:35 · 484 阅读 · 0 评论 -
hdu 4408 最小生成树计数
经典题了 #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define clr(a,v) memset(a,v,sizeo原创 2013-04-30 19:43:53 · 654 阅读 · 0 评论 -
hdu 4305 生成树计数
此题建图稍显麻烦,我先固定一个点,然后对斜率排序,斜率则写成了分数形式,然后嘛就是上模板 #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #de原创 2013-04-30 19:42:02 · 451 阅读 · 0 评论 -
ural 1627
与论文题没什么区别 #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace s原创 2013-04-30 19:38:52 · 502 阅读 · 0 评论 -
生成树计数 spoj 104 highways
参照基尔霍夫矩阵的定义建图 #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using n原创 2013-04-30 19:34:23 · 561 阅读 · 0 评论 -
数值积分 (辛普森公式 辛普森自适应法则) UVA
double F(double x) { } double simpson(double l, double r) { double m = (l + r) * 0.5; return (F(l) + 4 * F(m) + F(r)) * (r - l) / 6.0; } double asr(double l, double r, double eps, double A) { doubl原创 2013-04-23 20:41:33 · 1467 阅读 · 0 评论 -
SGU 502
白痴做法 : 状态压缩DP 其实直接搜更快 #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long ll; typedef pair pii; con原创 2013-04-15 20:27:46 · 577 阅读 · 0 评论 -
poj2154 polya theory
解法与ploya入门题一样,但要用 euler 函数优化 #include #include #include #include #include using namespace std; typedef long long ll; template T operator % (T a,T b) { return a>=b?a%b:a; } const int maxn = 111111; b原创 2013-03-27 00:24:20 · 426 阅读 · 0 评论 -
zoj 3547
这题没啥好说的,就是1^4+2^4+3^4+...+n^4 这个公式找出来就行了,接下来就是容斥原理进行计算,hdu 4488 已经告诉我们怎么去递推系数了 zoj 3547 #include #include #include #include #include #include #include using namespace std; typedef long long ll; co原创 2013-04-04 17:06:46 · 537 阅读 · 0 评论 -
hdu 2449
赤裸裸的高斯消元 切勿忘记选取主元素,hdu 的数据可以让不选取主元素的程序wrong(看来出现了一些恶心人的小分数) import java.util.*; import java.math.*; import java.io.*; public class Main { public static void main(String[] args) { new Task().ma原创 2013-03-27 03:00:44 · 525 阅读 · 0 评论 -
UVA 11806 容斥原理
#include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; const int maxn = 511; const int mod = 1000007; int C[原创 2013-05-04 14:07:04 · 525 阅读 · 0 评论 -
SGU 473
/* ID: BryatnMark LANG: JAVA TASK: AC */ import java.util.*; import java.io.*; import java.math.*; import java.util.*; public class Solution { public static void main(String args[]) throws Exceptio原创 2013-05-17 13:01:40 · 541 阅读 · 0 评论 -
lightoj 1236 Pairs Forming LCM 与 liver archive 6153 LCM Pair Sum
第一题求的是有多少正整数对 (i,j) (1 解决方法 质因数分解,求出每个质因数的指数,c1,c2,c3..... cm ans=(c1*2+1)*(c2*2+1)....(cm*2+1)+1,ans/=2; 第二题求的是所有正整数对的和,n很大很大,最后结果取模。 解决方法 设质因数依次为 p1,p2,p3,p4..... pm 指数相应为 a1,a2,a原创 2013-03-06 20:59:33 · 888 阅读 · 0 评论 -
扩展BABT STEP hdu 2815 poj 3243
http://hi.baidu.com/aekdycoin/item/236937318413c680c2cf29d4 这里已经说的很清楚了 A^X = B MOD(C) 若 gcd(A,C)!=1 令 d = gcd(A,C) B/=d, C/=d 重复此过程 我们得到 d1,d2,d3,d4............原创 2013-04-19 21:34:08 · 489 阅读 · 0 评论 -
poj 2417 lightoj 1325 求解离散对数 BABY STEP
poj 2417 HASH light oj 1325 map 红黑树硬搞#include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef long long ll; const in原创 2013-04-19 15:25:41 · 474 阅读 · 0 评论