![](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 · 666 阅读 · 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 · 560 阅读 · 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 · 416 阅读 · 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 · 944 阅读 · 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 · 629 阅读 · 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 · 327 阅读 · 0 评论 -
POJ 3372 —— Candy Distribution
#include#include#include#include#include#include#includeusing 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 · 455 阅读 · 0 评论 -
UVA 10288 - Coupons
求 n/1+n/2+...n/n 用最简分数表示/*ID: BryatnMarkLANG: JAVATASK: 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 · 571 阅读 · 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 · 665 阅读 · 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 · 691 阅读 · 0 评论 -
POJ 2103 Jackpot
/*ID: BryatnMarkLANG: JAVATASK: 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 · 592 阅读 · 0 评论 -
POJ System
/*ID: BryatnMarkLANG: JAVATASK: 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 · 410 阅读 · 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 · 926 阅读 · 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 · 421 阅读 · 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 · 390 阅读 · 0 评论 -
HDU 4611
#include using namespace std;typedef long long LL;templateT gcd(T a, T b) { return b ? gcd(b, a % b) : a;}templateT 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 · 510 阅读 · 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 · 578 阅读 · 0 评论 -
hdu 4447
多项式 计算 模拟 多项式乘法 和 多项式除法就行了,多项式除法可以参考高精度除法#include#include#include#include#include#include#include#include#include#include#includeusing namespace std;typedef long long LL;const int原创 2013-06-29 06:10:39 · 829 阅读 · 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 · 709 阅读 · 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 · 488 阅读 · 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 · 684 阅读 · 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 · 461 阅读 · 0 评论 -
UVA 11246 K-Multiple Free set
/*ID: BryatnMarkLANG: JAVATASK: 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 · 450 阅读 · 0 评论 -
UVA 10236 The Fibonacci Primes
/*ID: BryatnMarkLANG: JAVATASK: 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 · 484 阅读 · 0 评论 -
UVA 10162 - Last Digit
/*ID: BryatnMarkLANG: JAVATASK: 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 · 982 阅读 · 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 · 471 阅读 · 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 · 645 阅读 · 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 · 444 阅读 · 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 · 496 阅读 · 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 · 555 阅读 · 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 · 1457 阅读 · 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 · 568 阅读 · 0 评论 -
poj2154 polya theory
解法与ploya入门题一样,但要用 euler 函数优化#include#include#include#include#includeusing namespace std;typedef long long ll;templateT operator % (T a,T b){ return a>=b?a%b:a;}const int maxn = 111111;b原创 2013-03-27 00:24:20 · 421 阅读 · 0 评论 -
zoj 3547
这题没啥好说的,就是1^4+2^4+3^4+...+n^4 这个公式找出来就行了,接下来就是容斥原理进行计算,hdu 4488 已经告诉我们怎么去递推系数了zoj 3547#include#include#include#include#include#include#includeusing namespace std;typedef long long ll;co原创 2013-04-04 17:06:46 · 529 阅读 · 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 · 515 阅读 · 0 评论 -
UVA 11806 容斥原理
#include#include#include#include#include#include#include#include#include#include#include#include#include#includeusing namespace std;const int maxn = 511;const int mod = 1000007;int C[原创 2013-05-04 14:07:04 · 516 阅读 · 0 评论 -
SGU 473
/*ID: BryatnMarkLANG: JAVATASK: 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 · 532 阅读 · 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 · 878 阅读 · 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 · 475 阅读 · 0 评论 -
poj 2417 lightoj 1325 求解离散对数 BABY STEP
poj 2417 HASH light oj 1325 map 红黑树硬搞#include#include#include#include#include#include#include#include#include#include#includeusing namespace std;typedef long long ll;const in原创 2013-04-19 15:25:41 · 464 阅读 · 0 评论