栏目分类:
子分类:
返回
名师互学网用户登录
快速导航关闭
当前搜索
当前分类
子分类
实用工具
热门搜索
名师互学网 > IT > 面试经验 > 面试问答

zoj 2573 DNA Problem

面试问答 更新时间: 发布时间: IT归档 最新发布 模块sitemap 名妆网 法律咨询 聚返吧 英语巴士网 伯小乐 网商动力

zoj 2573 DNA Problem

#include<cstdio>#include<cstdlib>#include<cmath>#include<algorithm>#include<cstring>using namespace std;struct mat{unsigned long long a[51][51],n;mat(){memset(a,0,sizeof(a));}}A,B,C;int n,m,a,b;int MAT1[2][55],MAT2[2][55];mat mult(mat A,mat B){mat C;int n=A.n;C.n=n;for (int i=0;i<n;i++)for (int j=0;j<n;j++){C.a[i][j]=0;for (int k=0;k<n;k++)C.a[i][j]=(C.a[i][j]+A.a[i][k]*B.a[k][j])%2;}return C;}mat power(mat A,int n){mat ans,bas=A;ans.n=A.n;for (int i=0;i<A.n;i++)for (int j=0;j<A.n;j++){if (i==j) ans.a[i][j]=1;else ans.a[i][j]=0;}while (n>0){if (n&1) ans=mult(ans,bas);bas=mult(bas,bas);n>>=1;}return ans;}int main(){while (scanf("%d%d%d%d",&n,&m,&a,&b)>0){int tmp=m/(a+b);A.n=n;B.n=n;for (int i=0;i<n;i++)for (int j=0;j<n;j++) scanf("%d",&A.a[j][i]);for (int i=0;i<n;i++)for (int j=0;j<n;j++) scanf("%d",&B.a[j][i]);mat Q1=power(A,a),Q2=power(B,b);mat Q3=mult(Q1,Q2);mat Q4=power(Q3,tmp);mat Q5=power(A,min(m-tmp*(a+b),a)),Q6=power(B,max(m-tmp*(a+b)-a,0));mat Q7=mult(Q4,Q5);mat Q8=mult(Q7,Q6);for (int i=0;i<n;i++){char st[10];scanf("%s",st);if (st[0]=='a') MAT1[0][i]=1;else MAT1[0][i]=0;}memset(MAT2,0,sizeof(MAT2));for (int i=0;i<n;i++)for (int k=0;k<n;k++)MAT2[0][i]=(MAT2[0][i]+MAT1[0][k]*Q8.a[k][i])%2;if (MAT2[0][0]==1) printf("alive");else printf("dead");for (int i=1;i<n;i++)if (MAT2[0][i]==1) printf(" alive");else printf(" dead");puts("");}}
转载请注明:文章转载自 www.mshxw.com
本文地址:https://www.mshxw.com/it/378047.html
我们一直用心在做
关于我们 文章归档 网站地图 联系我们

版权所有 (c)2021-2022 MSHXW.COM

ICP备案号:晋ICP备2021003244-6号