矩阵+acm竞赛模版.docVIP

下载本文档

0
0
约2.44千字
约 3页
2021-07-10 发布于安徽
举报
版权申诉

矩阵+acm竞赛模版.doc

1、本文档共3页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

define MAXN 100#define fabs(x) ((x)0?(x):-(x))#define zero(x) (fabs(x)1e-10)struct mat{??? int n,m;??? double data[MAXN][MAXN];};int mul(mat c,const mat a,const mat b){??? int i,j,k;??? if (a.m!=b.n)??? ??? return 0;??? c.n=a.n,c.m=b.m;??? for (i=0;ic.n;i++)??? ??? for (j=0;jc.m;j++)??? ??? ??? for (c.data[j]=k=0;ka.m;k++)??? ??? ??? ??? c.data[j]+=a.data[k]*b.data[k][j];??? return 1;}int inv(mat a){??? int i,j,k,is[MAXN],js[MAXN];??? double t;??? if (a.n!=a.m)??? ??? return 0;??? for (k=0;ka.n;k++){??? ??? for (t=0,i=k;ia.n;i++)??? ??? ??? for (j=k;ja.n;j++)??? ??? ??? ??? if (fabs(a.data[j])t)??? ??? ??? ??? ??? t=fabs(a.data[is[k]=i][js[k]=j]);??? ??? if (zero(t))??? ??? ??? return 0;??? ??? if (is[k]!=k)??? ??? ??? for (j=0;ja.n;j++)??? ??? ??? ??? t=a.data[k][j],a.data[k][j]=a.data[is[k]][j],a.data[is[k]][j]=t;??? ??? if (js[k]!=k)??? ??? ??? for (i=0;ia.n;i++)??? ??? ??? ??? t=a.data[k],a.data[k]=a.data[js[k]],a.data[js[k]]=t;??? ??? a.data[k][k]=1/a.data[k][k];??? ??? for (j=0;ja.n;j++)??? ??? ??? if (j!=k)??? ??? ??? ??? a.data[k][j]*=a.data[k][k];??? ??? for (i=0;ia.n;i++)??? ??? ??? if (i!=k)??? ??? ??? ??? for (j=0;ja.n;j++)??? ??? ??? ??? ??? if (j!=k)??? ??? ??? ??? ??? ??? a.data[j]-=a.data[k]*a.data[k][j];??? ??? for (i=0;ia.n;i++)??? ??? ??? if (i!=k)??? ??? ??? ??? a.data[k]*=-a.data[k][k];??? }??? for (k=a.n-1;k=0;k--){??? ??? for (j=0;ja.n;j++)??? ??? ??? if (js[k]!=k)??? ??? ??? ??? t=a.data[k][j],a.data[k][j]=a.data[js[k]][j],a.data[js[k]][j]=t;??? ??? for (i=0;ia.n;i++)??? ??? ??? if (is[k]!=k)??? ??? ??? ??? t=a.data[k],a.data[k]=a.data[is[k]],a.data[is[k]]=t;??? }??? return 1;}double det(const mat a){??? int i,j,k,sign=0;??? double b[MAXN][MAXN],ret=1,t;??? if (a.n!=a.m)??? ??? return 0;??? for (i=0;ia.n;i++)??? ??? for (j=0;ja.m;j++)??? ??? ??? b[j]=a.data[j];??? for (i=0;ia.n;i++){??? ??? if (zero(b)){??? ??? ??? for (j=i+1;ja.n;j++)??? ??? ??? ??? if (!zero(b[j]))?