
数论
文章平均质量分 51
敬宅神社
(x^2 + (9/4)y^2 + z^2 - 1)^3 - x^2z^3 - (9/80)y^2z^3 = 0
展开
-
Uva 1635
题意:大概就是输出C(n,k)不能整除 m看来LRJ 的解析,开始不知道输出方式0.0一直WA ,还以为哪里爆了int,结果.............可以优化0.0代码:#include#include#include#include#includeusing namespace std;typedef unsigned long long UL原创 2015-05-11 20:20:23 · 340 阅读 · 0 评论 -
SGU 106
扩展欧几里得的应用但是精度确实是被卡成Dog#include#include#includeusing namespace std;typedef __int64 LL;LL gcd(LL a,LL b){return b==0?a:gcd(b,a%b);}void gcdex(LL a, LL b, LL &d, LL &x, LL &y) { if(!b) d = a,原创 2015-05-23 20:45:29 · 365 阅读 · 0 评论 -
LightOJ 1282
题意 :求 n 的 k 次方的 前 3位数 和 后三位数;后三位直接 快速幂去摸;前三位:( n^ k / 10^ ? )可以 直接 Lg(n) * k - int (Lg(n) * k )得出的是 N ^ K的除 10 ^ y后的小数 对 10 取对数的结果 ans;这时候反过来 pow (10, ans ) 得到这个数再 * 100 解决了#in原创 2015-05-19 20:39:33 · 342 阅读 · 0 评论 -
LightOJ 1259
题意:求一个偶数有多少个二元组质数相加正好等于其本身 哥德巴赫猜想0.0卡了空间,筛选质数用 bool 数组 0.0#include#include#include#include#include#include#include#include#include#define lson l , m, rt << 1#define rson m+1, r, rt << 1|原创 2015-05-19 21:14:13 · 397 阅读 · 0 评论 -
LightOJ 1341
题意 :一个数可以表示 N 对 数相乘(x1,x2)输入另一个数b 求有多少对x1 x2 俊大于等于 b求标准分解式的指数 因子数为 II(ai + 1)(表示不出数学符号.............)代码:#include#include#include#includeusing namespace std;typedef long long LL;co原创 2015-05-19 19:17:43 · 335 阅读 · 0 评论 -
Uva 10791
题意:输入N,求至少两个数 使其lcm值为 N#include#includeusing namespace std;typedef long long LL;const LL INF = 0x7fffffff;const int M = sqrt(INF + 0.5);LL P[M];int count;void Get(LL n) {count =原创 2015-05-07 19:41:27 · 330 阅读 · 0 评论 -
Uva 12716
题意:输入N , 求有多少对(a,b) gcd(a,b) == a ^ b然而 gcd(a,b) == a ^ b = c 可以得出 b = a - c;然后就是枚举c a 了#include#include#include#includeusing namespace std;const int maxn = 30000000 +131;typedef long原创 2015-05-07 20:47:09 · 278 阅读 · 0 评论 -
Light Oj 1236
题意:求1-n有多少对数的LCM正好为 n先求 N 的标准分解试p1 ^ a1 * p2 ^ a2 .......... pk ^ ak;有LCM 的可知 (a, b) 中必必有一个 含有 pi ^ ai;另一个则是 pi ^ 0 . . ...... pi ^ ai;所以一共有 NUM = II (ai +1 ) * 2 - 1又因为 对原创 2015-05-20 20:31:28 · 342 阅读 · 0 评论 -
Uva 10375
题意:求解C(m,n) / C(r,s) 利用每个数的唯一分解试,然后用一个辅助数组写出素数的指数最后在统一计算(开始手打快速幂,发现我SB,会出现负数)#include#include#include#includeusing namespace std;const int maxn = 5000;int num[maxn<<2],primes[maxn],e[m原创 2015-05-06 22:27:53 · 319 阅读 · 0 评论 -
Uva 12169
http://blog.youkuaiyun.com/u013451221/article/details/38497029感谢原博主提供的思想。带扩展二字的算法都可以出奇迹啊!#includeusing namespace std;typedef unsigned long long LL;LL num[300];void Gcd(LL a,LL b,LL &gcd,LL &x,LL &转载 2015-05-06 21:24:29 · 467 阅读 · 0 评论 -
UVA 11582
二元组有循环节,找出循环,在快速幂求出第几个,输出即可, (LL 大法好)#include#includeusing namespace std;typedef unsigned long long LL;const int maxn = 1000 + 5;int M[maxn];int F[maxn][maxn<<2];int pow_mod(LL a,LL b,int mod)原创 2015-05-06 20:12:09 · 323 阅读 · 0 评论 -
Uva 11181
题意:有n个人去超市买东西Pi 为 i 买东西的概率逛完后 得知有 r 人买了东西输出每人买东西的概率有R 人买了东西的事件为 E , 第 I 人买了东西的 事件为 Ai题目求的是 P(Ai | E ) = P ( Ai E ) / P(E)P(E)即为 正好有r 人买了东西的概率和而 一组r 人买了东西概率 为 Bi 第 i 人买了东西的 概率为 B原创 2015-05-12 21:46:00 · 272 阅读 · 0 评论 -
UVA 10491
题意:门里面放牛或者车有a 头牛, b两车;有人帮你开了c 个有牛的门求概率0.0分情况:你选的是车还是不是车的情况#include#includeusing namespace std;int main() { double a,b,c; while(cin >> a >> b >> c ) { printf("%.5lf\n",(a*b+b*(原创 2015-05-12 20:42:22 · 324 阅读 · 0 评论 -
Uva 1636
题意:第一枪开出去没有子弹要使后一枪木有子弹,怎么操作0.0所以两枪木有子弹是 00子串还有随机转一下, 转完后是0的概率为 NUM(0) / NUM(len )所以就是比NUM(00) / 【NUM(00) + NUM(01)】 和 NUM (0) / NUM(len)的大小0.0#include#include#includeusing namespace std原创 2015-05-12 20:17:38 · 334 阅读 · 0 评论 -
Uva 10820
就是求有多少组(x,y) 二元组 x ,y 互质先x 《 y 再 y 《 x 0.0phi(1-n)0.0:#includeusing namespace std;typedef long long LL;const int maxn = 50000 + 131;LL PHI[maxn];void Phi(int n,LL *phi) { for(int i = 2原创 2015-05-11 22:18:01 · 333 阅读 · 0 评论 -
LightOJ 1336
http://blog.youkuaiyun.com/lxpaopao/article/details/45417517只能说神代码,666666,佩服!!转载 2015-05-24 20:50:24 · 464 阅读 · 0 评论