C语言算法集

  • O3_383018
    了解作者
  • 14.3KB
    文件大小
  • rar
    文件格式
  • 0
    收藏次数
  • VIP专享
    资源类型
  • 0
    下载次数
  • 2022-05-29 08:41
    上传日期
C语言算法集 第一章 适合初学算法和计算机编程的学生
CH1.rar
  • CH1
  • ABINT0.DAT
    84B
  • ABAND0.C
    660B
  • AGSDL.C
    733B
  • ATLVS0.C
    303B
  • AGMIV0.C
    744B
  • AGGJE0.DAT
    160B
  • ACHOL.C
    1.4KB
  • ABAND.C
    1.6KB
  • ALDLE0.DAT
    185B
  • AGMQR0.DAT
    486B
  • AGRAD0.C
    504B
  • ACJDN.C
    2.3KB
  • AGMIV0.DAT
    245B
  • AGMIV.C
    346B
  • AGAUS0.C
    431B
  • ABINT0.C
    603B
  • ALDLE.C
    1.7KB
  • AGJDN.C
    1.7KB
  • ATLVS0.DAT
    120B
  • ACGAS.C
    2.5KB
  • AGSDL0.DAT
    80B
  • ACHOL0.C
    418B
  • ACJDN0.C
    873B
  • Atrde0.c
    313B
  • ALDLE0.C
    471B
  • AGGJE.C
    1.6KB
  • AGMQR0.C
    793B
  • ATRDE0.DAT
    100B
  • AGAUS.C
    1.7KB
  • AGRAD.C
    1.1KB
  • AGGJE0.C
    736B
  • ACGAS0.DAT
    152B
  • AGJDN0.DAT
    144B
  • AGSDL0.C
    399B
  • ACJDN0.DAT
    284B
  • ATRDE.C
    631B
  • ABAND0.DAT
    416B
  • AGMQR.C
    579B
  • ABINT.C
    1.2KB
  • ACGAS0.C
    692B
  • ACHOL0.DAT
    148B
  • AGJDN0.C
    496B
  • ATLVS.C
    1KB
  • AGRAD0.DAT
    76B
  • AGAUS0.DAT
    72B
内容介绍
#include "stdlib.h" #include "stdio.h" int acgas(ar,ai,n,br,bi) int n; double ar[],ai[],br[],bi[]; { int *js,l,k,i,j,is,u,v; double p,q,s,d; js=malloc(n*sizeof(int)); for (k=0;k<=n-2;k++) { d=0.0; for (i=k;i<=n-1;i++) for (j=k;j<=n-1;j++) { u=i*n+j; p=ar[u]*ar[u]+ai[u]*ai[u]; if (p>d) {d=p;js[k]=j;is=i;} } if (d+1.0==1.0) { free(js); printf("err**fail\n"); return(0); } if (is!=k) { for (j=k;j<=n-1;j++) { u=k*n+j; v=is*n+j; p=ar[u]; ar[u]=ar[v]; ar[v]=p; p=ai[u]; ai[u]=ai[v]; ai[v]=p; } p=br[k]; br[k]=br[is]; br[is]=p; p=bi[k]; bi[k]=bi[is]; bi[is]=p; } if (js[k]!=k) for (i=0;i<=n-1;i++) { u=i*n+k; v=i*n+js[k]; p=ar[u]; ar[u]=ar[v]; ar[v]=p; p=ai[u]; ai[u]=ai[v]; ai[v]=p; } v=k*n+k; for (j=k+1;j<=n-1;j++) { u=k*n+j; p=ar[u]*ar[v]; q=-ai[u]*ai[v]; s=(ar[v]-ai[v])*(ar[u]+ai[u]); ar[u]=(p-q)/d; ai[u]=(s-p-q)/d; } p=br[k]*ar[v]; q=-bi[k]*ai[v]; s=(ar[v]-ai[v])*(br[k]+bi[k]); br[k]=(p-q)/d; bi[k]=(s-p-q)/d; for (i=k+1;i<=n-1;i++) { u=i*n+k; for (j=k+1;j<=n-1;j++) { v=k*n+j; l=i*n+j; p=ar[u]*ar[v]; q=ai[u]*ai[v]; s=(ar[u]+ai[u])*(ar[v]+ai[v]); ar[l]=ar[l]-p+q; ai[l]=ai[l]-s+p+q; } p=ar[u]*br[k]; q=ai[u]*bi[k]; s=(ar[u]+ai[u])*(br[k]+bi[k]); br[i]=br[i]-p+q; bi[i]=bi[i]-s+p+q; } } u=(n-1)*n+n-1; d=ar[u]*ar[u]+ai[u]*ai[u]; if (d+1.0==1.0) { free(js); printf("err**fail\n"); return(0); } p=ar[u]*br[n-1]; q=-ai[u]*bi[n-1]; s=(ar[u]-ai[u])*(br[n-1]+bi[n-1]); br[n-1]=(p-q)/d; bi[n-1]=(s-p-q)/d; for (i=n-2;i>=0;i--) for (j=i+1;j<=n-1;j++) { u=i*n+j; p=ar[u]*br[j]; q=ai[u]*bi[j]; s=(ar[u]+ai[u])*(br[j]+bi[j]); br[i]=br[i]-p+q; bi[i]=bi[i]-s+p+q; } js[n-1]=n-1; for (k=n-1;k>=0;k--) if (js[k]!=k) { p=br[k]; br[k]=br[js[k]]; br[js[k]]=p; p=bi[k]; bi[k]=bi[js[k]]; bi[js[k]]=p; } free(js); return(1); } 
评论
    相关推荐