C语言数学函数相关操作

今天主要介绍C语言关于数学函数几个的相关操作。我用的编译器是tcc。实在没必要用一些大型的编译器，tcc小巧实用，实为居家旅行必备之神器！

1.求任意两个数之间的素数，我加了一个素数判断的函数，很实用。code:

/*求任意两数之间的素数*/
/*long 为 2^31-1=2,147,483,647*/
#include <stdio.h>
#include <math.h>  /*double sqrt(double)*/
int isPrime(long num); 
void AtoBPrime(long n1,long n2);
int main()
{
　　 if(isPrime(102))
 　　{
 　　　　 printf("Prime");
 　　}
 　　else
 　　{
  　　　　printf("Not Prime!");
　　 }
　　 AtoBPrime(100,10000);
 　　return 0;
}
/*判断是否是素数，1-是素数，0-不是素数，输入为任意整数*/
int isPrime(long num)
{
　　 int flag = 0;
　　 int mov = 2;
 　　long high;
 　　if(num <= 1)
 　　{
　　 　　 flag = 0;
　　 }
　　 if(2 == num)
 　　{
  　　flag = 1;
 　　}
 　　high = (long)sqrt(num);
 　　for(; mov<high; mov++)
　　 {
 　　 if(num % mov == 0)
 　　{
   　　flag = 0;
  　　 mov = num;
  　　}
  　　else if(high != mov+1)
  　　{
  　　 continue;
 　　 }
 　　 else
 　　 {
 　　  flag = 1;
 　　 }
50　　  }
 
 　　return flag;
}
/*求任意两个数之间的素数*/
void AtoBPrime(long n1,long n2)
{
 long mov;
 int cnt=0;
 printf("the array prime:\n");
 for(mov=n1;mov<n2;mov++)
 {
  if(isPrime(mov))
  {
   cnt++;
   if(cnt%5 == 0)
   {
    printf("\n");
   }
   printf("%d ",mov);
  }
 }
}

2.回文数。

/*Palindrome:回文数,指正读，反读都是一样的数*/
#include <stdio.h>
#include <stdlib.h> /*char* ltoa (long, char*, int)*/
#include <string.h> /*size_t strlen(char *)*/
#include <math.h>
int isPalindrome(long num); 
void AtoBPalindrome(long n1,long n2);
int main()
{
     if(isPalindrome(101))
     {
          printf("Palindrome\n");
     }
     else
     {
         printf("Not Palindrome!\n");
     }
     AtoBPalindrome(10,10000);
     return 0;
}
/*判断是否是素数，1-是回文数，0-不是回文数，输入为任意整数*/
int isPalindrome(long num)
{
 int i,len,half;
 int flag=0;
 char arr[20];  /*字符型数组，存放数的各个位上的数字*/
 ltoa(num,arr,10);
 /***************************
 char* ltoa (long, char*, int);  
 long：指定装换的长整形数字
 char*:接受字符串
 int:进制属性
 ****************************/
 len = strlen(arr);
 half = len/2;
 for(i=0;i<=half;i++)
 {
  if(arr[i] != arr[--len])
  {
   flag = 0;
   break;
  }
  if(i >= half)
  {
   flag = 1;
  }
  else
  {
   flag = 0;
  }
 }
 return flag;
}
/*打印三层回文数*/
void AtoBPalindrome(long n1,long n2)
{
 long mov;
 for(mov=n1 ; mov<n2 ;mov++)
 {
  if(isPalindrome(mov)&&isPalindrome(mov*mov)&&isPalindrome(mov*mov*mov))
  {
   printf("n=%d\n",mov);
  }
 }
}

3.神奇的6174

/*神奇6174，找到一个四位数，从大到小排列和从小到大排列所得数的差为它自己*/
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
int isSame(long num); 
void AtoBSame(long n1,long n2);
int main()
{
 AtoBSame(1000,10000);
 return 0;
}
int isSame(long num)
{
 int flag = 0;
 int arr[4];
 int i=0,j,temp;
 int n1,n2;//重新组合数字
 /*解析数字*/
 arr[0] = num/1000;
 arr[1] = num/100%10;
 arr[2] = num/10%10;
 arr[3] = num%10;
 
 /*排序,从大到小*/
 for(;i<3;i++)
 {
  for(j=i+1;j<4;j++)
  {
   if(arr[i]<arr[j])
   {
    temp = arr[j];
    arr[j] = arr[i];
    arr[i] = temp;
   }
  }
 }
 n1 = arr[0]*1000+arr[1]*100+arr[2]*10+arr[3];
 n2 = arr[3]*1000+arr[2]*100+arr[1]*10+arr[0];
 /*两个数是相差是否是一个数*/
 if(num == n1-n2)
 {
  flag = 1;
 }
 return flag;
}
/*求能够满足条件的特殊数*/
void AtoBSame(long n1,long n2)
{
 long mov;
 for(mov=n1;mov<n2;mov++)
 {
  if(isSame(mov))
  {
   printf("num=%d\n",mov);
  }
 }
}

4.数字e，与最大公约数，最小公倍数。

/*实用的数字操作*/
#include <stdio.h>
#include <math.h>
float efun();/*求的数字e*/
long divisor(long n1,long n2);/*最大公约数*/
long multiple(long n1,long n2);/*最小公倍数*/
int main()
{
 float e=0.0f;
 long d,m;
 
 e = efun();
 printf("%8f\n",e);
 /*最大公约数与最小公倍数测试语句*/
 d = divisor(9,100);
 m = multiple(10,101);
 printf("d=%d\nm=%d",d,m);
 return 0;
}
float efun()
{
 float e=1.0f,n=1.0f;
 int i = 1;
 while(1/n > 1e-10)
 {
  e += 1/n;
  i++;
  n = i*n;
 }
 return e;
}
/*求最大公约数函数*/
long divisor(long n1,long n2)
{ 
 long t,c,d; 
 /*说明n1比较大*/
 if(n1<n2)
 {
  t = n1;
  n1 = n2;
  n2 = t;
 }
 /*求最大公约数*/
 c = n1%n2;
 while(c != 0)
 {
  n1 = n2;
  n2 = c;
  c = n1%n2;
 }
 d = n2;/*d为最大公约数*/
 return d;
}
/*求小公倍数函数*/
long multiple(long n1,long n2)
{
 long c,m;
 c = n1*n2;
 m = c/divisor(n1,n2);
 return m;
}

5.牛顿法解方程。

/*牛顿法解方程*/
#include <stdio.h>
#include <math.h>  /*fabs(float)*/
float Newton(float a , float b , float c , float d);
int main()
{
 float x;
 x = Newton(1,2,3,4);
 printf("x=%.5f\n",x);
 return 0;
}
/*牛顿解法函数形式*/
float Newton(float a , float b , float c , float d)
{
 float x=1.0f,x0,f,f1;
 do
 {
  x0 = x;
  f = ((a*x0+b)*x0+c)*x0+d;  /*函数本身*/
  f1 = (3*a*x0+2*b)*x0+c;  /*导数*/
  x = x0 - f/f1;    /*迭代公式*/
 }while(fabs(x-x0) >= 1e-5);
 return x;
}
 

6.计算系统运行时间

/*计算系统运行时间*/
#include <stdio.h>
#include <conio.h>/*kbhit(void)*/
/*********************************************
kbhit(void)
功 能及返回值:检查当前是否有键盘输入，若有则返回一个非0值，否则返回0
用 法：int kbhit(void); 
*********************************************/
#include <stdlib.h>/*void _sleep (unsigned long);tcc*/
typedef struct time
{
 int hour;
 int minute;
 int second;
}time;
int main()
{
 time t = {0,0,0};
 while(!kbhit())
 {
  _sleep(1000);/*1s=1000ms*/
  if(t.second == 59)
  {
   t.minute = t.minute + 1;
   if(t.minute == 60)
   {
    t.hour = t.hour+1;
    t.minute = 0;
   }
   t.second = 0;
  }
  else
  {
   t.second = t.second+1;
  }
  printf("%d:%d:%d\n",t.hour,t.minute,t.second);
 }
}