三角形数

Description

一定数目的点或圆在等距离的排列下可以形成一个等边三角形，这样的数被称为三角形数。比如10个点可以组成一个等边三角形，因此10是一个三角形数：

　　　　ｘ

　　　ｘ　ｘ

　　ｘ　ｘ　ｘ

　ｘ　ｘ　ｘ　ｘ

开始18个三角形数是1、3、6、10、15、21、28、36、45、55、66、78、91、105、120、136、153、171。

请求出第n个三角形数。

Input

输入数据包含多组测试用例。

每组测试用例为一个整数N，长度不超100。

Output

对于每组测试用例，输出对应的三角形数。

Sample Input

1

2

3

4

5

11111111

Sample Output

1

3

6

10

15

61728399382716

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdlib>
using namespace std;
const int MAXSIZE = 101;
void Add(char *str1, char *str2, char *str3);
void Minus(char *str1, char *str2, char *str3);
void Mul(char *str1, char *str2, char *str3);
void Div(char *str1, char *str2, char *str3);
int main()
{
//    freopen("in.txt","r",stdin);
char str1[MAXSIZE],str2[MAXSIZE],str3[MAXSIZE*MAXSIZE],str11[MAXSIZE],str22[MAXSIZE];
while(scanf("%s",str1)!=EOF)
{
int i=0;
memset(str2,'0',sizeof(str2));
memset(str3,'0',sizeof(str3));
memset(str11,'0',sizeof(str2));
memset(str22,'0',sizeof(str2));
while(str1[i])i++;
if((str1[i-1]-48)%2==0)
{
Div(str1,"2",str11);
Add(str1,"1",str2);
Mul(str11,str2,str3);
puts(str3);
}
else
{
Add(str1,"1",str2);
Div(str2,"2",str22);
Mul(str1,str22,str3);
puts(str3);
}
}
return 0;
}
void Add(char *str1, char *str2, char *str3)
{// str3 = str1 + str2;
int  i, j, i1, i2, tmp, carry;
int  len1 = strlen(str1), len2 = strlen(str2);
char  ch;
i1 = len1-1; i2 = len2-1;
j = carry = 0;
for( ; i1 >= 0 && i2 >= 0; ++j, --i1, --i2 ){
tmp = str1[i1]-'0'+str2[i2]-'0'+carry;
carry = tmp/10;
str3[j] = tmp%10+'0';
}
while( i1 >= 0 ){
tmp = str1[i1--]-'0'+carry;
carry = tmp/10;
carry = tmp/10;
str3[j++] = tmp%10+'0';
}
if( carry ) str3[j++] = carry+'0';
str3[j] = '\0';
for( i=0, --j; i < j; ++i, --j ){
ch = str3[i]; str3[i] = str3[j]; str3[j] = ch;
}
}
void Mul(char *str1, char *str2, char *str3){
int  i, j, i1, i2, tmp, carry, jj;
int  len1 = strlen(str1), len2 = strlen(str2);
char  ch;
jj = carry = 0;
for( i1=len1-1; i1 >= 0; --i1 ){
j = jj;
for( i2=len2-1; i2 >= 0; --i2, ++j ){
tmp =
(str3[j]-'0')+(str1[i1]-'0')*(str2[i2]-'0')+carry;
if( tmp > 9 ){
carry = tmp/10; str3[j] = tmp%10+'0';
}
else {
str3[j] = tmp+'0'; carry = 0;
}
}
if( carry ) {
str3[j] = carry+'0'; carry = 0; ++j;
}
++jj;
}
--j;
while( str3[j] == '0' && j > 0 ) --j;
str3[++j] = '\0';
for( i=0, --j; i < j; ++i, --j ){
ch = str3[i]; str3[i] = str3[j]; str3[j] = ch;
}
}
void Div(char *str1, char *str2, char *str3){
int  i1, i2, i, j, jj, tag, carry, cf, c[MAXSIZE];
int  len1 = strlen(str1), len2 = strlen(str2), lend;
char  d[MAXSIZE];
memset(c, 0, sizeof(c));
memcpy(d, str1, len2);
lend = len2; j = 0;
for( i1=len2-1; i1 < len1; ++i1 ){
if( lend < len2 ){
d[lend] = str1[i1+1]; c[j] = 0;
++j; ++lend;
}
else if( lend == len2 ){
jj = 1;
for( i=0; i < lend; ++i ){
if( d[i] > str2[i] ) break;
else if( d[i] < str2[i] ){
jj = 0; break;
}
}
if( jj == 0 ){
d[lend] = str1[i1+1]; c[j] = 0;
++j; ++lend;
continue;
}
}
if( jj==1 || lend > len2 ){
cf = jj=0;
while( d[jj] <= '0' && jj < lend ) ++jj;
if( lend-jj > len2 ) cf = 1;
else if( lend-jj < len2 ) cf = 0;
else{
i2 = 0; cf = 1;
for( i=jj; i < lend; ++i ){
if( d[i] < str2[i2] ){
cf = 0; break;
}
else if( d[i] > str2[i2] ){
break;
}
++i2;
}
}//else
while( cf ){
i2 = len2-1; cf = 0;
for( i=lend-1; i >= lend-len2; --i ){
d[i] = d[i]-str2[i2]+'0';
if( d[i] < '0' ){
d[i] = d[i]+10;    carry = 1;
--d[i-1];
}
else carry = 0;
--i2;
}
++c[j]; jj=0;
while( d[jj] <= '0' && jj < lend ) ++jj;
if( lend-jj > len2 ) cf = 1;
else if( lend-jj < len2 ) cf = 0;
else{
i2 = 0; cf = 1;
for( i=jj; i < lend; ++i ){
if( d[i] < str2[i2] ){
cf = 0; break;
}
else if( d[i] > str2[i2] ){
break;
}
++i2;
}
}//else
}//while
jj = 0;
while( d[jj] <= '0' && jj < lend ) ++jj;
for( i=0;i < lend-jj; ++i ) d[i] = d[i+jj];
d[i] = str1[i1+1];    lend = i+1;
++j;
}//else
}//for
i = tag = 0;
while( c[i] == 0 ) ++i;
for( ; i < j; ++i, ++tag ) str3[tag] = c[i]+'0';
str3[tag] = '\0';
}