<Sicily>Polynomial
一、题目描述
Given a polynomial and the value of the variable x, you task is to calculate the value of the polynomial.
二、输入
The first line of the input is a positive integer T. T is the number of test cases followed.
In each test case, the first line contains two integer n and x (0<=n<=50, 0<=x<=10000). The second line contains n+1 integers, representing the coefficients of a polynomial from power degree N down to power degree 0, each integer is no less than 0 and no more than 10000.
三、输出
The output of each test case should consist one line, containing the result of the polynomial.
例如:
输入:
1
2 5
1 2 3
输出:
38
四、解题思路
题意:这道题要求的是多项式值。如果只是简单的求多项式的值很简单。可题目有个条件是 (0<=n<=50, 0<=x<=10000),这样求出来的数已经超过了long long类型。
开始不知道怎么处理这个大数问题,在这里卡住了。后来同学说使用数组处理。以前没接触过使用输出保存大数的问题,当时也上不了网,于是自己思考一下,使用一个数组n表示n位十进制。每位表示十进制中的一位数。大数+(int类型)常数;大数*(int类型)常数。
1、多项式计算方法
记得高中学过求多项式值得方法能较少计算量:
如多项式:1*5^2+2*5^1+3*5^0
2、数组表示整数
例如439保存在数组bigInteger[3];bigInteger[2]=4;bigInteger[1]=3;bigInteger[0]=9;
1)数组表示的整数加上一个常数num
例如加上35。
先从最低位开始,例如bigInteger[0]+35=43;bigInteger[0]=43%10;进位数carryInt=43/10=4;
向高位逐步计算
bigInteger[1]+4=7,bigInteger[1]=7%10,carryInt=7/10。
当出现carryInt为0是停止计算,否则不断想高位计算。
计算的结果是:
bigInteger[2]=4;bigInteger[1]=7;bigInteger[0]=4;
2)数组表示的整数乘以一个常数num
例如bigInteger[2]=4;bigInteger[1]=7;bigInteger[0]=4;乘以4。
也是从最低位开始bigInteger[0]*4=16,bigInteger[0]=16%10=6;
进位数carryInt=16/10=1;
继续到下一位(高位):bigInteger[1]*4=28,28+carryInt=29,bigInteger[1]=29%10=9,
carryInt=29/10=2;
继续下一位:bigInteger[2]*4=16,16+carryInt=18,bigInteger[2]=18%10=8,
carryInt=18/10=1;
继续下一位:
bigInteger[3]开始至为0,所以bigInteger[3]=bigInteger[3]+carryInt=1。
最后计算结果为:bigInteger[3]=1,bigInteger[2]=8,bigInteger[1]=9,bigInteger[0]=6
五、代码
#include<iostream>
using namespace std;
const int BIT_NUM = 200; //支持大数的最大位数(200位)
int main()
{
int times;
cin >> times;
while(times--)
{
int bigInteger[BIT_NUM] = {0}; //用数组保存一个大数
int n, x, argue; //n-多项式的最高次数,x-位置数x的值,argue-每一项的系数
long long sum = 0; //多项式的结果
cin >> n >> x;
for(int k = 0; k < (n + 1); k++) //计算多项式,这里使用了一个公式例如1*5^2+2*5^1+3*5^0=38 可以写成(1*5+2)*5+3=38
{
cin >> argue;
int carryInt = 0; //进位(从低位算起,例如各位数是7乘以一个常数8,那么结果各位数是6,向高位数进5(carryInt))
for(int i = BIT_NUM; i > 0; i--) //一个大数乘以一个常数,从各位开始逐个与常数相乘,满十向高位进
{
bigInteger[i-1] = bigInteger[i-1] * x + carryInt;
carryInt = bigInteger[i-1] / 10;
bigInteger[i-1] %= 10;
}
//大数加上一个常数
int addCarry = 0;
bigInteger[BIT_NUM - 1] += argue; //大数加上一个常数argue,先个位数加
for(int i = BIT_NUM; i > 0; i--)
{
bigInteger[i-1] = bigInteger[i-1] + addCarry;
addCarry = bigInteger[i-1] / 10; //addCarry保存进制数
bigInteger[i-1] %= 10;
if(addCarry == 0) break; //如果出现不需要向高位数进的退出
}
}
int flag = false;
for(int i = 0; i < BIT_NUM; i++) //从高位第一个非0开始输出大数
{
if(bigInteger[i] != 0) {flag = true;}
if(flag) {cout << bigInteger[i];}
}
cout << endl;
}
return 0;
}
<Sicily>Polynomial的更多相关文章
- Polynomial Library in OpenCascade
Polynomial Library in OpenCascade eryar@163.com 摘要Abstract:分析幂基曲线即多项式曲线在OpenCascade中的计算方法,以及利用OpenSc ...
- sicily 中缀表达式转后缀表达式
题目描述 将中缀表达式(infix expression)转换为后缀表达式(postfix expression).假设中缀表达式中的操作数均以单个英文字母表示,且其中只包含左括号'(',右括号‘)’ ...
- sicily 1934. 移动小球
Description 你有一些小球,从左到右依次编号为1,2,3,...,n. 你可以执行两种指令(1或者2).其中, 1 X Y表示把小球X移动到小球Y的左边, 2 X Y表示把小球X移动到小球Y ...
- 周赛-Integration of Polynomial 分类: 比赛 2015-08-02 08:40 10人阅读 评论(0) 收藏
Integration of Polynomial Time Limit: 2000/1000MS (Java/Others) Memory Limit: 128000/64000KB (Java/O ...
- FZU 2215 Simple Polynomial Problem(简单多项式问题)
Description 题目描述 You are given an polynomial of x consisting of only addition marks, multiplication ...
- Project Euler 101 :Optimum polynomial 最优多项式
Optimum polynomial If we are presented with the first k terms of a sequence it is impossible to say ...
- 大数求模 sicily 1020
Search
- 【数论】UVa 10586 - Polynomial Remains
Problem F: Polynomial Remains Given the polynomial a(x) = an xn + ... + a1 x + a0, compute the remai ...
- Sicily 1510欢迎提出优化方案
这道题我觉得是除1000(A-B)外最简单的题了……不过还是提出一个小问题:在本机用gcc编译的时候我没包括string.h头文件,通过编译,为什么在sicily上却编译失败? 1510. Mispe ...
随机推荐
- Unity3D——加入剑痕效果(PocketRPG Trail插件)
首先非常感谢大家的支持,因为近期项目吃紧,所以更新的速度可能会有点慢!希望大家谅解,当然大家的支持是我最大的动力.我也会尽我所能写出更好的文章,当然因为本人是个新手并且工作的内容也不是unity3D. ...
- android studio 、 as 如何导入eclipse项目
安卓项目有两种,一种是eclipse开发的,一种的android studio开发的.有些在github开源的安卓项目,下载下来之后不知道该如何处理了. 这个是Eclipse安卓项目的目录结构. 这个 ...
- [JZOJ 100026] [NOIP2017提高A组模拟7.7] 图 解题报告 (倍增)
题目链接: http://172.16.0.132/senior/#main/show/100026 题目: 有一个$n$个点$n$条边的有向图,每条边为$<i,f(i),w(i)>$,意 ...
- BZOJ 1898 构造矩阵+矩阵快速幂
思路: T的最小公倍数是12 那么12以内暴力 整除12 的部分用矩阵快速幂 //By SiriusRen #include <cstdio> #include <cstring&g ...
- JavaScript实现记住密码功能
用js实现记住密码功能,但是前端记住密码不安全,最好还是不要用.我感觉这个记住密码应该是通过与后台建立一个会话来实现. 这个效果的测试地址在:http://ofoyou.com/blog/rePass ...
- winform控件命名规范对照表
WinForm Control 命名规范 数据类型 数据类型简写 标准命名举例 Label lbl lblMessage LinkLabel llbl llblToday Button btn btn ...
- MEF基本概念
基本概念: Managed Extensibility Framework 或 MEF 是一个用于创建可扩展的轻型应用程序的库,在.NET 4.0发布 Container:容器,使用Compositi ...
- Uncaught TypeError: undefined is not a function
index.html <script src="resources/sap-ui-core.js" id="sap-ui-bootstrap" data- ...
- hiho 1590 - 紧张的会议室。区间问题
题目链接 小Hi的公司最近员工增长迅速,同时大大小小的会议也越来越多:导致公司内的M间会议室非常紧张. 现在小Hi知道公司目前有N个会议,其中第i个会议的时间区间是(Si, Ei). 注意这里时间区间 ...
- SLAM概念学习之随机SLAM算法
这一节,在熟悉了Featue maps相关概念之后,我们将开始学习基于EKF的特征图SLAM算法. 1. 机器人,图和增强的状态向量 随机SLAM算法一般存储机器人位姿和图中的地标在单个状态向量中,然 ...