UVA 10123(欧拉定理 + 大数类)
来源:互联网 发布:地税局数据采集手续 编辑:程序博客网 时间:2024/06/18 03:47
题目链接:
https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=1154
题解:
公式:ans = 1 + C(2,n) + C(4,n);
贴一个大数类的模板
大数类:
class integer{ friend ostream& operator<<(ostream&, const integer&); friend istream& operator>>(istream&, integer&); friend int compare(const integer, const integer); friend bool operator<(const integer, const integer); friend bool operator<=(const integer, const integer); friend integer operator<<(integer, int); friend integer operator+(integer, integer); friend integer operator-(integer, integer); friend integer operator*(integer, integer); friend integer operator/(integer, integer); friend integer operator^(integer, int);public: integer() { }; integer(int a) { signbit = (a >= 0 ? PLUS : MINUS); a = abs(a); while (a) { digits.push_back(a % 10); a /= 10; } }; ~integer() { }; void zero_justify(void);private: // 表示大数的结构。数位按低位到高位的顺序存放。存储的是数位值,而不是数位的字符值。 // 数字 1 用 ASCII 码为 1 的字符,而不是 ‘1’ 来表示。 vector < char > digits; // 存放数位值。 int signbit; // 符号位。};integer operator+(integer, integer);integer operator-(integer, integer);// 重载输出符号<<。ostream& operator<<(ostream& os, const integer& n){ // 当为自然数时,不输出“+”号。 os << (n.signbit > 0 ? "" : "-"); for (int i = n.digits.size() - 1; i >= 0; i--) // 存放的是数位值,输出时要转换为数位的字符值。 os << (char)('0' + n.digits[i]); return os;}// 重载输入符号>>。istream& operator>>(istream& in, integer& n){ // 读入表示大整数的字符串。 string s; in >> s; if (s.length() == 0) { // 若长度为0,表示读入的是一个空字符串,应该要给予错误提示,但对于本题来说,不 // 存在这种情况,故不考虑。若要完善,需要进一步处理此种情况。 } // 在解析之前,理应判断该字符串是否是一个合法的大数表示,即开始一位为符号位,为“+”或 // “-”,正整数“+”可以忽略,后面应该全部是数字字符,且不包含前导0。由于本题不存在此种情 // 况,故不考虑,若要完善需进一步处理此种情况。 // 解析符号位。 if (isdigit(s[0])) n.signbit = PLUS; else { n.signbit = (s[0] == '+' ? PLUS : MINUS); s.erase(s.begin()); } // 解析数字位,从低位到高位存放于动态数组中。 n.digits.clear(); for (int i = s.length() - 1; i >= 0; i--) n.digits.push_back((char)(s[i] - '0')); return in;}// 移除大数运算产生的前导0。void integer::zero_justify(void){ for (int i = digits.size() - 1; i >= 1; i--) if (digits[i] == 0) digits.erase(digits.begin() + i); else break; if (digits.size() == 1 && digits[0] == 0) signbit = PLUS;}// 比较两个数的大小,若相等,返回 0,若 a 大于 b,返回 1,若 a 小于 b,返回 -1。int compare(const integer a, const integer b){ // 若 a 为正,b 为负,则有 a 大于 b。 if (a.signbit == PLUS && b.signbit == MINUS) return PLUS; // 若 a 为负,b 为正,则 a 小于 b。 if (a.signbit == MINUS && b.signbit == PLUS) return MINUS; // 若两数符号相同,则 a 的数位长度大于 b 的数位长度,若 a 为正数,则有 a 大于 b,若 a // 为负数,有 a 小于 b,则将 a 的符号位乘以1作为结果返回即可。 if (a.digits.size() > b.digits.size()) return PLUS * a.signbit; // 若两数符号相同,若a的数位长度小于b的数位长度,若 a 为正数,则 a 小于 b,若 a 为负数, // 则 a 大于 b,则将 a 的符号位乘以 -1 作为结果返回即可。 if (a.digits.size() < b.digits.size()) return MINUS * a.signbit; // 两个数的数位长度相等,符号相等,则从低位到高位逐个比较数位值的大小。 int marker = a.digits.size() - 1; while (marker >= 0) { if (a.digits[marker] > b.digits[marker]) return PLUS * a.signbit; if (a.digits[marker] < b.digits[marker]) return MINUS * a.signbit; marker--; } return 0;}// 小于比较符。bool operator<(const integer a, const integer b){ return compare(a, b) < 0;}// 小于等于比较符。bool operator<=(const integer a, const integer b){ return compare(a, b) <= 0;}// 左移操作。b 表示左移的位数。即将 a 乘以 10^b。integer operator<<(integer a, int b){ integer c = a; // 若 a 等于 0,则移位后仍为 0。 if (a.digits.size() == 1 && a.digits[0] == 0) return c; // 向左移动 b 位,相当于在数位数组前插入 b 个 0。 for (int i = 0; i < b; i++) c.digits.insert(c.digits.begin(), (char)(0)); return c;}// 加法运算。从低位往高位进行,每一位相加产生的进位向高位传递,注意负数的处理,可以相应的// 转变为减法来简化操作。integer operator+(integer a, integer b){ integer c; // 两数符号位相等,则最终符号位和原来的加数符号位相同。 if (a.signbit == b.signbit) c.signbit = a.signbit; else { // 若两个加数为一正一负,则将其转换为减法操作。 if (a.signbit == MINUS) { a.signbit = PLUS; c = b - a; } else { b.signbit = PLUS; c = a - b; } return c; } // 从动态数组的开始即低位开始相加,逐渐往前进位,如果有的话。 int carry = 0; // 进位。 int maxsize = max(a.digits.size(), b.digits.size()); int marker = 0; // 为两数添加前导 0,以使得数位相同,便于计算。 while (a.digits.size() < maxsize) a.digits.push_back(char(0)); while (b.digits.size() < maxsize) b.digits.push_back(char(0)); // 逐位相加。 while (marker < maxsize) { int t = a.digits[marker] + b.digits[marker] + carry; carry = t / 10; c.digits.push_back((char)(t % 10)); marker++; } // 若最后一次进位为 1,则需在最高位加 1。 if (carry) c.digits.push_back((char)(1)); c.zero_justify(); return c;}// 减法。为了保证借位能够成功终止,最好确保被减数不小于减数。同样可以在适当情况转变为加法来简化操作。integer operator-(integer a, integer b){ integer c; // 当某一个数为负数时,转换为加法。 if (a.signbit == MINUS || b.signbit == MINUS) { b.signbit = MINUS * b.signbit; c = a + b; return c; } // 两个都为正数,若 a < b,则尽量保证被减数大于减数。 if (a < b) { c = b - a; c.signbit = MINUS; return c; } else c.signbit = PLUS; int borrow = 0; // 借位。 int maxsize = a.digits.size(); // 被减数大于减数。 int marker = 0; // 计数器。 // 为减数添加前导 0,便于计算。因为先前已经使得 a > b 且 a 和 b 均为正数。 while (b.digits.size() < maxsize) b.digits.push_back(char(0)); // 从低位开始逐位相减,不够的向高位借位。 while (marker < maxsize) { int v = a.digits[marker] - borrow - b.digits[marker]; if (v < 0) { v += 10; borrow = 1; } else borrow = 0; c.digits.push_back((char)(v % 10)); marker++; } c.zero_justify(); return c;}// 乘法,采用一行一行算的方法比单纯的反复加法更快且编写代码不会变复杂许多。每次操作都把第一个// 乘数左移一个数位,然后把左移后的数乘以x的积加到最终结果中,其中x是第二个乘数对应的数字。尽// 管用反复加法来实现这一步看起来不够精巧,但由于对于第二个乘数的每个数字,内层循环最多执行 9// 次,实际上并没有浪费太多的时间,左移一个数位等价于用它乘以进位制的基数,在十进制里就是乘以 10。integer operator*(integer a, integer b){ integer c = integer(0), row = a; row.signbit = PLUS; for (int i = 0; i < b.digits.size(); i++) { for (int j = 0; j < b.digits[i]; j++) c = (c + row); row = row << 1; } c.signbit = a.signbit * b.signbit; c.zero_justify(); return c;}// 除法运算。将余数和商不断左移,余数加上被除数的下一个数字,然后尝试减去除数。integer operator/(integer a, integer b){ integer c, row; row.signbit = PLUS; // 商的符号为 a 的符号乘以 b 的符号。 c.signbit = a.signbit * b.signbit; int v; // 从数a的最高位开始除。 for (int i = a.digits.size() - 1; i >= 0; i--) { row = row << 1; row.digits[0] = a.digits[i]; // 当除数小于被除数,减去除数,对应位的商增 1。 v = 0; while (b <= row) { v++; row = row - b; } // 将商左移 1 位,刚得到的商即为最低位。 c = c << 1; c.digits[0] = (char)(v); } c.zero_justify(); return c;}// 幂计算。幂就是反复作乘法,利用:a^n = a^(n / 2) * a^(n / 2) * a^(n mod 2),可以减少乘// 的次数。integer operator^(integer a, int b){ if (b == 0) return integer(1); if (b == 1) return a; if (b == 2) return a * a; return ((a ^ (b / 2)) ^ 2) * (a ^ (b % 2));}
AC代码:
#include <iostream>#include <vector>#include <algorithm>using namespace std;#define PLUS 1#define MINUS -1class integer{ friend ostream& operator<<(ostream&, const integer&); friend istream& operator>>(istream&, integer&); friend int compare(const integer, const integer); friend bool operator<(const integer, const integer); friend bool operator<=(const integer, const integer); friend integer operator<<(integer, int); friend integer operator+(integer, integer); friend integer operator-(integer, integer); friend integer operator*(integer, integer); friend integer operator/(integer, integer); friend integer operator^(integer, int);public: integer() { }; integer(int a) { signbit = (a >= 0 ? PLUS : MINUS); a = abs(a); while (a) { digits.push_back(a % 10); a /= 10; } }; ~integer() { }; void zero_justify(void);private: vector < char > digits; int signbit;};integer operator+(integer, integer);integer operator-(integer, integer);ostream& operator<<(ostream& os, const integer& n){ os << (n.signbit > 0 ? "" : "-"); for (int i = n.digits.size() - 1; i >= 0; i--) os << (char)('0' + n.digits[i]); return os;}istream& operator>>(istream& in, integer& n){ string s; in >> s; if (s.length() == 0) { } if (isdigit(s[0])) n.signbit = PLUS; else { n.signbit = (s[0] == '+' ? PLUS : MINUS); s.erase(s.begin()); } n.digits.clear(); for (int i = s.length() - 1; i >= 0; i--) n.digits.push_back((char)(s[i] - '0')); return in;}void integer::zero_justify(void){ for (int i = digits.size() - 1; i >= 1; i--) if (digits[i] == 0) digits.erase(digits.begin() + i); else break; if (digits.size() == 1 && digits[0] == 0) signbit = PLUS;}int compare(const integer a, const integer b){ if (a.signbit == PLUS && b.signbit == MINUS) return PLUS; if (a.signbit == MINUS && b.signbit == PLUS) return MINUS; if (a.digits.size() > b.digits.size()) return PLUS * a.signbit; if (a.digits.size() < b.digits.size()) return MINUS * a.signbit; int marker = a.digits.size() - 1; while (marker >= 0) { if (a.digits[marker] > b.digits[marker]) return PLUS * a.signbit; if (a.digits[marker] < b.digits[marker]) return MINUS * a.signbit; marker--; } return 0;}bool operator<(const integer a, const integer b){ return compare(a, b) < 0;}bool operator<=(const integer a, const integer b){ return compare(a, b) <= 0;}integer operator<<(integer a, int b){ integer c = a; if (a.digits.size() == 1 && a.digits[0] == 0) return c; for (int i = 0; i < b; i++) c.digits.insert(c.digits.begin(), (char)(0)); return c;}integer operator+(integer a, integer b){ integer c; if (a.signbit == b.signbit) c.signbit = a.signbit; else { if (a.signbit == MINUS) { a.signbit = PLUS; c = b - a; } else { b.signbit = PLUS; c = a - b; } return c; } int carry = 0; int maxsize = max(a.digits.size(), b.digits.size()); int marker = 0; while (a.digits.size() < maxsize) a.digits.push_back(char(0)); while (b.digits.size() < maxsize) b.digits.push_back(char(0)); while (marker < maxsize) { int t = a.digits[marker] + b.digits[marker] + carry; carry = t / 10; c.digits.push_back((char)(t % 10)); marker++; } if (carry) c.digits.push_back((char)(1)); c.zero_justify(); return c;}integer operator-(integer a, integer b){ integer c; if (a.signbit == MINUS || b.signbit == MINUS) { b.signbit = MINUS * b.signbit; c = a + b; return c; } if (a < b) { c = b - a; c.signbit = MINUS; return c; } else c.signbit = PLUS; int borrow = 0; int maxsize = a.digits.size(); int marker = 0; while (b.digits.size() < maxsize) b.digits.push_back(char(0)); while (marker < maxsize) { int v = a.digits[marker] - borrow - b.digits[marker]; if (v < 0) { v += 10; borrow = 1; } else borrow = 0; c.digits.push_back((char)(v % 10)); marker++; } c.zero_justify(); return c;}integer operator*(integer a, integer b){ integer c = integer(0), row = a; row.signbit = PLUS; for (int i = 0; i < b.digits.size(); i++) { for (int j = 0; j < b.digits[i]; j++) c = (c + row); row = row << 1; } c.signbit = a.signbit * b.signbit; c.zero_justify(); return c;}integer operator/(integer a, integer b){ integer c, row; row.signbit = PLUS; c.signbit = a.signbit * b.signbit; int v; for (int i = a.digits.size() - 1; i >= 0; i--) { row = row << 1; row.digits[0] = a.digits[i]; // 当除数小于被除数,减去除数,对应位的商增 1。 v = 0; while (b <= row) { v++; row = row - b; } // 将商左移 1 位,刚得到的商即为最低位。 c = c << 1; c.digits[0] = (char)(v); } c.zero_justify(); return c;}integer operator^(integer a, int b){ if (b == 0) return integer(1); if (b == 1) return a; if (b == 2) return a * a; return ((a ^ (b / 2)) ^ 2) * (a ^ (b % 2));}int main(){ int s; integer n; scanf("%d", &s); while(s--) { cin >> n; cout << ((n ^ 4) - integer(6) * (n ^ 3) + integer(23) * (n ^ 2) - integer(18) * n + integer(24)) / integer(24) << endl; } return 0;}
0 0
- UVA 10123(欧拉定理 + 大数类)
- 计算几何 ( 平面欧拉定理+大数运算 )——UVA 10213
- uva 1342 欧拉定理(计算几何模板)
- uva 1342 - That Nice Euler Circuit(欧拉定理)
- UVA 10200 (欧拉定理打表)
- fzoj1759(欧拉定理)
- 【hdu2462】【数论】【欧拉函数+欧拉定理+大数快速幂】The Luckiest number
- UVa 10692 Huge Mods(快速幂)(欧拉定理)
- 【UVA】1342 - That Nice Euler Circuit(几何+欧拉定理)
- UVA 10213 How Many Pieces of Land? 欧拉定理
- UVa 10692 - Huge Mod (数论 + 欧拉定理扩展)
- UVa 10213 How Many Pieces of Land (欧拉定理)
- 费马小定理&&欧拉定理
- 威尔逊定理,费马小定理,欧拉定理
- SUM 大数取余 欧拉定理 二项展开 数论
- 大数取模:一般取模+技巧取模+快速幂取模+欧拉函数(费马小定理)
- poj3358数论(欧拉定理)
- POJ 3358 (欧拉定理)
- Select()函数的介绍
- Linux的发展史
- hdoj5301
- MVP设计模式初步学习
- Java主要特性
- UVA 10123(欧拉定理 + 大数类)
- sdut oj1138 数据结构上机测试2-1:单链表操作A
- UIsearchBar的placeholder字体颜色和大小
- IPhone 自定义 UITableViewCell 行缩进
- 转载 感谢原作者 矩阵分解在推荐系统中的应用:NMF和经典SVD实战
- NOI 2016
- tjut 4622
- Buffer is too small
- 了解PCB