ACM小组内部预定函数 数学问题：

最新推荐文章于 2024-06-08 09:53:18 发布
最新推荐文章于 2024-06-08 09:53:18 发布
阅读量672 收藏 1
点赞数
分类专栏： ACM 文章标签： math.h c null fft function up

 

ACM小组内部预定函数

Ver 2.0 by IcyFenix       

一、数学问题

1.精度计算——大数阶乘

语法：int result=factorial(int n);

参数：

n：

n 的阶乘

返回值：

阶乘结果的位数

注意：

 

 

本程序直接输出n!的结果，需要返回结果请保留long a[]

 

需要 math.h

源程序：

 

 

int factorial(int n)
{
long a[10000];
int i,j,l,c,m=0,w;

a[0]=1;
for(i=1;i<=n;i++)
    {
    c=0;
    for(j=0;j<=m;j++)
        {
        a[j]=a[j]*i+c;
        c=a[j]/10000;
        a[j]=a[j]%10000;
    }
    if(c>0) {m++;a[m]=c;}
}

w=m*4+log10(a[m])+1;
printf("/n%ld",a[m]);
for(i=m-1;i>=0;i--) printf("%4.4ld",a[i]);
return w;
}


2.精度计算——乘法（大数乘小数）

语法：mult(char c[],char t[],int m);

参数：

c[]：

被乘数，用字符串表示，位数不限

t[]：

结果，用字符串表示

m：

乘数，限定10以内

返回值：

null

注意：

 

 

需要 string.h

源程序：

 

 

void mult(char c[],char t[],int m)
{
    int i,l,k,flag,add=0;
    char s[100];
    l=strlen(c);
    for (i=0;i<l;i++)
        s[l-i-1]=c[i]-'0';

    for (i=0;i<l;i++)
           {
           k=s[i]*m+add;
           if (k>=10) {s[i]=k%10;add=k/10;flag=1;} else {s[i]=k;flag=0;add=0;}
           }
    if (flag) {l=i+1;s[i]=add;} else l=i;

    for (i=0;i<l;i++)
        t[l-1-i]=s[i]+'0';
    t[l]='/0';
}

 

3.精度计算——乘法（大数乘大数）

语法：mult(char a[],char b[],char s[]);

参数：

a[]：

被乘数，用字符串表示，位数不限

b[]：

乘数，用字符串表示，位数不限

t[]：

结果，用字符串表示

返回值：

null

注意：

 

 

空间复杂度为 o(n^2)

 

需要 string.h

源程序：

 

 

void mult(char a[],char b[],char s[])
{
    int i,j,k=0,alen,blen,sum=0,res[65][65]={0},flag=0;
    char result[65];
    alen=strlen(a);blen=strlen(b);

    for (i=0;i<alen;i++)
    for (j=0;j<blen;j++) res[i][j]=(a[i]-'0')*(b[j]-'0');

    for (i=alen-1;i>=0;i--)
        {
            for (j=blen-1;j>=0;j--) sum=sum+res[i+blen-j-1][j];
            result[k]=sum%10;
            k=k+1;
            sum=sum/10;
        }

    for (i=blen-2;i>=0;i--)
        {
            for (j=0;j<=i;j++) sum=sum+res[i-j][j];
            result[k]=sum%10;
            k=k+1;
            sum=sum/10;
        }
    if (sum!=0) {result[k]=sum;k=k+1;}

    for (i=0;i<k;i++) result[i]+='0';
    for (i=k-1;i>=0;i--) s[i]=result[k-1-i];
    s[k]='/0';

    while(1)
        {
        if (strlen(s)!=strlen(a)&&s[0]=='0')
            strcpy(s,s+1);
        else
            break;
        }
}

 

4.精度计算——加法

语法：add(char a[],char b[],char s[]);

参数：

a[]：

被乘数，用字符串表示，位数不限

b[]：

乘数，用字符串表示，位数不限

t[]：

结果，用字符串表示

返回值：

null

注意：

 

 

空间复杂度为 o(n^2)

 

需要 string.h

源程序：

 

 

void add(char a[],char b[],char back[])
{
    int i,j,k,up,x,y,z,l;
    char *c;
    if (strlen(a)>strlen(b)) l=strlen(a)+2; else l=strlen(b)+2;
    c=(char *) malloc(l*sizeof(char));
    i=strlen(a)-1;
    j=strlen(b)-1;
    k=0;up=0;
    while(i>=0||j>=0)
        {
            if(i<0) x='0'; else x=a[i];
            if(j<0) y='0'; else y=b[j];
            z=x-'0'+y-'0';
            if(up) z+=1;
            if(z>9) {up=1;z%=10;} else up=0;
            c[k++]=z+'0';
            i--;j--;
        }
    if(up) c[k++]='1';
    i=0;
    c[k]='/0';
    for(k-=1;k>=0;k--)
        back[i++]=c[k];
    back[i]='/0';
}

 

5.精度计算——减法

语法：sub(char s1[],char s2[],char t[]);

参数：

s1[]：

被减数，用字符串表示，位数不限

s2[]：

减数，用字符串表示，位数不限

t[]：

结果，用字符串表示

返回值：

null

注意：

 

 

默认s1>=s2，程序未处理负数情况

 

需要 string.h

源程序：

 

 

void sub(char s1[],char s2[],char t[])
{
    int i,l2,l1,k;
    l2=strlen(s2);l1=strlen(s1);
    t[l1]='/0';l1--;
    for (i=l2-1;i>=0;i--,l1--)
        {
        if (s1[l1]-s2[i]>=0)
            t[l1]=s1[l1]-s2[i]+'0';
        else
            {
            t[l1]=10+s1[l1]-s2[i]+'0';
            s1[l1-1]=s1[l1-1]-1;
            }
        }
    k=l1;
    while(s1[k]<0) {s1[k]+=10;s1[k-1]-=1;k--;}
    while(l1>=0) {t[l1]=s1[l1];l1--;}
loop:
    if (t[0]=='0')
        {
        l1=strlen(s1);
        for (i=0;i<l1-1;i++) t[i]=t[i+1];
        t[l1-1]='/0';
        goto loop;
        }
    if (strlen(t)==0) {t[0]='0';t[1]='/0';}
}

 

6.任意进制转换

语法：conversion(char s1[],char s2[],long d1,long d2);

参数：

s[]：

原进制数字，用字符串表示

s2[]：

转换结果，用字符串表示

d1：

原进制数

d2：

需要转换到的进制数

返回值：

null

注意：

 

 

高于9的位数用大写'A'～'Z'表示，2～16位进制通过验证

源程序：

 

 

void conversion(char s[],char s2[],long d1,long d2)
{
    long i,j,t,num;
    char c;
    num=0;
    for (i=0;s[i]!='/0';i++)
        {
        if (s[i]<='9'&&s[i]>='0') t=s[i]-'0'; else t=s[i]-'A'+10;
        num=num*d1+t;
        }
    i=0;
    while(1)
        {
        t=num%d2;
        if (t<=9) s2[i]=t+'0'; else s2[i]=t+'A'-10;
        num/=d2;
        if (num==0) break;
        i++;
        }
    for (j=0;j<i/2;j++)
        {c=s2[j];s2[j]=s[i-j];s2[i-j]=c;}
    s2[i+1]='/0';
}

 

7.最大公约数、最小公倍数

语法：resulet=hcf(int a,int b)、result=lcd(int a,int b)

参数：

a：

int a，求最大公约数或最小公倍数

b：

int b，求最大公约数或最小公倍数

返回值：

返回最大公约数（hcf）或最小公倍数（lcd）

注意：

 

 

lcd 需要连同 hcf 使用

源程序：

 

 

int hcf(int a,int b)
{
    int r=0;
    while(b!=0)
        {
        r=a%b;
        a=b;
        b=r;
        }
    return(a);
}

lcd(int u,int v,int h)
{
    return(u*v/h);
}

 

8.组合序列

语法：m_of_n(int m, int n1, int m1, int* a, int head)

参数：

m：

组合数C的上参数

n1：

组合数C的下参数

m1：

组合数C的上参数，递归之用

*a：

1～n的整数序列数组

head：

头指针

返回值：

null

注意：

 

 

*a需要自行产生

 

初始调用时，m=m1、head=0

 

调用例子：求C(m,n)序列：m_of_n(m,n,m,a,0);

源程序：

 

 

void m_of_n(int m, int n1, int m1, int* a, int head)
{
    int i,t;
    if(m1<0 || m1>n1) return;

    if(m1==n1)
        {
        for(i=0;i<m;i++) cout<<a[i]<<' '; // 输出序列
        cout<<'/n';
        return;
        }
    m_of_n(m,n1-1,m1,a,head); // 递归调用
    t=a[head];a[head]=a[n1-1+head];a[n1-1+head]=t;
    m_of_n(m,n1-1,m1-1,a,head+1); // 再次递归调用
    t=a[head];a[head]=a[n1-1+head];a[n1-1+head]=t;
}

 

9.快速傅立叶变换（FFT）

语法：kkfft(double pr[],double pi[],int n,int k,double fr[],double fi[],int l,int il);

参数：

pr[n]：

输入的实部

pi[n]：

数入的虚部

n，k：

满足n=2^k

fr[n]：

输出的实部

fi[n]：

输出的虚部

l：

逻辑开关，0 FFT，1 ifFT

il：

逻辑开关，0 输出按实部/虚部；1 输出按模/幅角

返回值：

null

注意：

 

 

需要 math.h

源程序：

 

 

void kkfft(pr,pi,n,k,fr,fi,l,il)
int n,k,l,il;
double pr[],pi[],fr[],fi[];
{
    int it,m,is,i,j,nv,l0;
    double p,q,s,vr,vi,poddr,poddi;
    for (it=0; it<=n-1; it++)
        {
         m=it; is=0;
        for (i=0; i<=k-1; i++)
            {j=m/2; is=2*is+(m-2*j); m=j;}
        fr[it]=pr[is]; fi[it]=pi[is];
        }
    pr[0]=1.0; pi[0]=0.0;
    p=6.283185306/(1.0*n);
    pr[1]=cos(p); pi[1]=-sin(p);
    if (l!=0) pi[1]=-pi[1];
    for (i=2; i<=n-1; i++)
        {
       p=pr[i-1]*pr[1];
       q=pi[i-1]*pi[1];
        s=(pr[i-1]+pi[i-1])*(pr[1]+pi[1]);
        pr[i]=p-q; pi[i]=s-p-q;
        }
    for (it=0; it<=n-2; it=it+2)
        {
       vr=fr[it]; vi=fi[it];
        fr[it]=vr+fr[it+1]; fi[it]=vi+fi[it+1];
        fr[it+1]=vr-fr[it+1]; fi[it+1]=vi-fi[it+1];
        }
    m=n/2; nv=2;
    for (l0=k-2; l0>=0; l0--)
        {
        m=m/2; nv=2*nv;
        for (it=0; it<=(m-1)*nv; it=it+nv)
            for (j=0; j<=(nv/2)-1; j++)
                {
               p=pr[m*j]*fr[it+j+nv/2];
                q=pi[m*j]*fi[it+j+nv/2];
                s=pr[m*j]+pi[m*j];
                s=s*(fr[it+j+nv/2]+fi[it+j+nv/2]);
                poddr=p-q; poddi=s-p-q;
                fr[it+j+nv/2]=fr[it+j]-poddr;
                fi[it+j+nv/2]=fi[it+j]-poddi;
                fr[it+j]=fr[it+j]+poddr;
                fi[it+j]=fi[it+j]+poddi;
                }
        }
    if (l!=0)
        for (i=0; i<=n-1; i++)
            {
           fr[i]=fr[i]/(1.0*n);
            fi[i]=fi[i]/(1.0*n);
            }
    if (il!=0)
            for (i=0; i<=n-1; i++)
            {
           pr[i]=sqrt(fr[i]*fr[i]+fi[i]*fi[i]);
            if (fabs(fr[i])<0.000001*fabs(fi[i]))
                {
               if ((fi[i]*fr[i])>0) pi[i]=90.0;
                else pi[i]=-90.0;
                }
            else
                pi[i]=atan(fi[i]/fr[i])*360.0/6.283185306;
            }
    return;
}

 

10.Ronberg算法计算积分

语法：result=integral(double a,double b);

参数：

a：

积分上限

b：

积分下限

function f：

积分函数

返回值：

f在（a,b）之间的积分值

注意：

 

 

function f(x)需要自行修改，程序中用的是sina(x)/x

 

需要 math.h

 

默认精度要求是1e-5

源程序：

 

 

double f(double x)
{
    return sin(x)/x; //在这里插入被积函数
}

double integral(double a,double b)
{
    double h=b-a;
    double t1=(1+f(b))*h/2.0;
    int k=1;
    double r1,r2,s1,s2,c1,c2,t2;
loop:
    double s=0.0;
    double x=a+h/2.0;
    while(x<b)
        {
        s+=f(x);
        x+=h;
        }
    t2=(t1+h*s)/2.0;
    s2=t2+(t2-t1)/3.0;
    if(k==1)
      {
        k++;h/=2.0;t1=t2;s1=s2;
        goto loop;
        }
    c2=s2+(s2-s1)/15.0;
    if(k==2){
        c1=c2;k++;h/=2.0;
        t1=t2;s1=s2;
        goto loop;
        }
    r2=c2+(c2-c1)/63.0;
    if(k==3){
        r1=r2; c1=c2;k++;
        h/=2.0;
        t1=t2;s1=s2;
        goto loop;
        }
    while(fabs(1-r1/r2)>1e-5){
        r1=r2;c1=c2;k++;
        h/=2.0;
        t1=t2;s1=s2;
        goto loop;
        }
    return r2;
}

 

11.行列式计算

语法：result=js(int s[][],int n)

参数：

s[][]：

行列式存储数组

n：

行列式维数，递归用

返回值：

行列式值

注意：

 

 

函数中常数N为行列式维度，需自行定义

源程序：

 

 

int js(s,n)
int s[][N],n;
{
    int z,j,k,r,total=0;
    int b[N][N];/*b[N][N]用于存放，在矩阵s[N][N]中元素s[0]的余子式*/
    if(n>2)
        {
        for(z=0;z<n;z++)
            {
            for(j=0;j<n-1;j++)
                 for(k=0;k<n-1;k++)
                        if(k>=z) b[j][k]=s[j+1][k+1];  else b[j][k]=s[j+1][k];
            if(z%2==0) r=s[0][z]*js(b,n-1); /*递归调用*/
            else r=(-1)*s[0][z]*js(b,n-1);
            total=total+r;
            }
        }
    else if(n==2)
       total=s[0][0]*s[1][1]-s[0][1]*s[1][0];
    return total;
}

 

12.求排列组合数

语法：result=P(long n,long m); / result=long C(long n,long m);

参数：

m：

排列组合的上系数

n：

排列组合的下系数

返回值：

排列组合数

注意：

 

 

符合数学规则：m<＝n

源程序：

 

 

long P(long n,long m)
{
    long p=1;
    while(m!=0)
        {p*=n;n--;m--;}
    return p;
}

long C(long n,long m)
{
    long i,c=1;
    i=m;
    while(i!=0)
        {c*=n;n--;i--;}
    while(m!=0)
        {c/=m;m--;}
    return c;
}

 

maomao625111654
关注
专栏目录
ACM小组内部预定函数（详细）
05-31
ACM小组内部预定函数，涵盖了ACm大量的函数方法，而且很详细。适用于入门的ACM学习者
ACM小组内部预定函数
戎码人生
04-07 3090
ACM小组内部预定函数 Ver 2.0 by IcyFenix        数学问题：       1.精度计算——大数阶乘 2.精度计算——乘法（大数乘小数） 3.精度计算——乘法（大数乘大数） 4.精度计算——加法 5.精度计算——减法 6.任意进制转换
精度计算-大数阶乘
滨州学院csdn高校俱乐部 冶近玮
05-21 1090
精度计算-大数阶乘 本算法的目的在于计算一个比较大的数的阶乘，由于得到的结果比较大，是现有的数据类型无法存储的，所以我决定将结果存储在一个long a[]数组中。 我们的思路是把每4位数看做数组的一个元素来存储，例如：个、十、百、千存在a[0]，万、十万、百万、千万存在a[1]以此类推。 下面是我的C语言实现过程： int factorial(int n) { long a[1000
C语言 使用FFT得到信号幅度谱
pan的博客
06-06 8842
前面的一篇文章我们介绍了使用DFT得到信号的幅度谱的方法，现在我们来看一下FFT实现信号幅度谱。这里我们使用的FFT源程序是徐士良老师的C语言算法程序——快速傅里叶变换。另外，本文也会对DFT以及FFT实现信号幅度谱计算时，所需要的计算量以及花费的时间进行比较。 1.FFT原理 FFT（Fast Fourier Transformation）是离散傅氏变换（DFT）的快速算法，也就是快速傅里叶变换...
C语言实现FFT并与DFT比较（徐士良老师编写的c语言算法程序）
Eric_The_Red的博客
06-03 3975
一 . FFT方法说明 计算n个采样点 P=p0, p1, ..........., pn-1 的傅里叶变换，可以归结为计算多项式 : f(x)= p0+p1 x+ p2 x 2 +...... +pn-1 x n-1 在各n次单位根1，w，w 2，…,w n-1上的值，即 f0=p0+p1+..........+pn-1 f1=p0+p1w+...
如何入门参与数学建模？
热门推荐
JiandaoStudio
12-25 7万+
数学建模感想 纪念逝去的大学数学建模：两次校赛，两次国赛，两次美赛，一次电工杯。从大一下学期组队到现在，大三下学期，时间飞逝，我的大学建模生涯也告一段落。感谢建模路上帮助过我的学长和学姐们，滴水之恩当涌泉相报，写下这篇感想，希望可以给学弟学妹们一丝启发，也就完成我的想法了。拙劣的文笔，也不知道写些啥，按顺序随便写写吧。 我是怎么选择建模的： 大一上，第一次听到数学建模其实是大一上学期，not...
《青少年编程与数学》课程方案：1、课程目标
最新发布
明月看潮生的博客
06-08 1319
《青少年编程与数学》课程旨在为10岁以上青少年提供编程与数学学习路径，强调真实学习体验，避免"玩具式"教育。课程分为初期和中后期阶段，初期培养软件工程师基础能力，中后期进阶到人工智能领域高级软件工程师。同时，课程计划让学生提前完成数学知识学习，提高学业成绩。通过实践、持续学习、跨学科融合，激发青少年的好奇心和求知欲，培养逻辑思维和自学能力，最终实现自主成长管理和社会价值体现。
教案ACM小组内部预定函数.pdf
10-06
【ACM小组内部预定函数】涉及的是一系列与算法竞赛相关的数学、字符串处理、计算几何、数论和图论问题，以及高精度运算专题。这些知识点对于编程竞赛和高性能计算至关重要。 1. **精度计算**： - **大数阶乘**：在...
ACM小组内部预定函数.pdf
10-07
这篇文档是关于ACM竞赛中常见的编程问题和算法的总结，主要涵盖了数学问题、字符串处理、计算几何、数论以及图论等多个领域。以下是对这些知识点的详细解释： 1. **精度计算**： - **大数阶乘**：用于计算大整数的...
"ACM内部预定函数汇总：数学、字符串处理、计算几何、数论、图论全面覆盖
总的来说，ACM内部预订函数是ACM小组的利器，提供了丰富的数学函数、字符串处理函数、计算几何函数、数论函数和图论函数，可以帮助ACM小组内部成员解决各种复杂的问题，提高解题效率，是ACM小组学习、训练和比赛中的...
频谱分析：基于python画出时域频域波形
taw19960426的博客
09-26 1万+
一，FFT解释 FFT（Fast Fourier Transformation）是离散傅氏变换（DFT）的快速算法。即为快速傅氏变换。它是根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅立叶变换的算法进行改进获得的。 二、基于python画出时域频域波形 1.python代码 依据快速傅里叶算法得到信号的频域，将其各个频率分量的幅值绘制成图。为便于计算，我们将采样频率8000次近似设为8192次 ...
K8s 很难么？带你从头到尾捋一遍，不信你学不会
民工哥的博客
02-23 9154
点击关注公众号，回复“1024”获取2TB学习资源！为什么要学习 Kubernetes?虽然 Docker 已经很强大了，但是在实际使用上还是有诸多不便，比如集群管理、资源调度、文件管理等...
10天智能锁项目实战第1天（了解单片机STM32F401RET6和C语言基础）
北豼不太皮的博客
02-21 3670
10天智能锁项目实战第1天（了解单片机STM32F401RET6和C语言基础）一、学习目标二、了解单片机STM32F401RET6三、C语言基础 一、学习目标 二、了解单片机STM32F401RET6 4、STM32F401RE特征 三、C语言基础 1.数据类型 常用2的次方： 2^7 = 128 2^8 = 256 2^15 = 32768 2^16= 65536 51单片机常见的数据类型： char: 占内存1字节 取值范围：-128~127 -2^7 ~ 2
Android开发（1） | Fragment 的应用——新闻应用
Jormungand_V的博客
02-21 7155
文章目录 news_item.xml： <TextView xmlns:android="http://schemas.android.com/apk/res/android" android:id="@+id/news_title" android:layout_width="match_parent" android:layout_height="wrap_content" android:lines="1" android:ellipsize="
系统学习 TypeScript（四）——变量声明的初步学习
编程三昧
02-25 1785
认识了 TypeScript 中的基础类型，接下来当然是变量声明的相关学习了，我在这里记录一下我学习过程中的一些总结。
清除浮动的五种方法
weixin_52212950的博客
02-22 2773
为什么要清除浮动？因为往往浮动后的子元素因为脱离了标准流，不能自动撑起父元素的高度，所以会对后续布局产生影响，对此清除浮动不如理解为清除浮动产生的影响更为合理。 例如：我们设置了两个盒子如图所示，粉色为父盒子，只有宽度没有高度，蓝色盒子有宽有高，粉色盒子的高由蓝盒子的高度撑开。 但是给蓝色的子盒子设置了左浮动后，脱离了标准流，无法再撑开粉盒子，可以看到粉盒子高度变成了0； 清除浮动有五种方法： 直接设置父元素的高度 额外标签法 单伪元素法 双伪元素法 ove.
HTML实现学习网站首页
做技术，坚持才是最重要的，一时的热情只是暂时的进步
02-20 4248
项目访问 reset.css /* http://meyerweb.com/eric/tools/css/reset/ v2.0 | 20110126 License: none (public domain) */ html, body, div, span, applet, object, iframe, h1, h2, h3, h4, h5, h6, p, blockquote, pre, a, abbr, acronym, address, big, cite, code, del
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值