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

zoj 2853 Evolution

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

zoj 2853 Evolution

#include<cstdio>#include<cstdlib>#include<cctype>#include<cstring>#include<cmath>#include<ctime>#include<vector>#include<map>#include<set>#include<list>#include<deque>#include<stack>#include<bitset>#include<functional>#include<numeric>#include<utility>#include<sstream>#include<iomanip>#include<algorithm>#include<limits>#include<iostream>#include<fstream>#include<string.h>#include<string>#include<queue> using namespace std;#define max(a,b) ((a>b)?a:b) #define min(a,b) ((a<b)?a:b) #define Max(a,b,c) max(a,max(b,c)) #define Min(a,b,c) min(a,min(b,c)) #define Pi acos(-1.0) #define Clear(a) memset(a,0,sizeof(a)) #define setfalse(a) memset(a,false,sizeof(a)) #define settrue(a) memset(a,true,sizeof(a)) #define clrstr(a) memset(a,'',sizeof(a)) #define Case(a) printf("Case %d: ",a) #define caseh(a) printf("Case #%d: ",a)#define getcase(a) scanf("%d",&a) #define caseloop(a,b) for(a=1;a<=b;a++) #define EPS 1e-9 #define LL long long #define MX 50000 #define PB(x) push_back(x) #define PP pop_back() #define PF(x) push_front(x) #define PPF(x) pop_front() #define IN(x) insert(x)#define PS(x) push(wax) #define P(x) pop() template<class T>inline void checkmin(T &a,T b){if(b<a)a=b;}template<class T>inline void checkmax(T &a,T b){if(b>a)a=b;}const double eps=1e-15;const double pi=acos(-1.0);#define two(X) (1<<(X))#define twoL(X) (((int64)(1))<<(X))#define contain(S,X) (((S)&two(X))!=0)#define containm(S,X) (((S)&twoL(X))!=0)#define SIZE(A) ((int)A.size())#define LENGTH(A) ((int)A.length())#define MP(A,B) make_pair(A,B)#define MUL(a,b) (int)(((int64)a*b%MOD))#define FOR(i,a,b) for(i=(a);i<(b);i++)#define REP(i,n) for(int i=0;i<(n);i++)typedef pair<int,int>ipair;#define MOD (1000000007)using namespace std;const int size=203;#define ll doubleint n;struct two_matrix{ ll array[size][size];};two_matrix pw;two_matrix unit;two_matrix ret;int main(){int i,t,u,v,j,k;int x;double num[size],f;double ans;while(cin>>n>>x,n+x){for(i=0;i<n;i++){for(j=0;j<n;j++){if(i==j){pw.array[i][i]=1.0;unit.array[i][j]=1.0;}else{pw.array[i][i]=0.0;unit.array[i][j]=0.0;}}}for(i=0;i<n;i++){cin>>num[i];}cin>>t;for(i=0;i<t;i++){cin>>u>>v>>f;unit.array[u][u]-=f;unit.array[v][u]+=f;}while(x>0){if(x&1)for(i=0;i<n;i++)for(j=0;j<n;j++){ret.array[i][j]=0;for(k=0;k<n;k++)ret.array[i][j]=(ret.array[i][j]+(pw.array[i][k]*unit.array[k][j]));}for(i=0;i<n;i++)for(j=0;j<n;j++)pw.array[i][j]=ret.array[i][j];for(i=0;i<n;i++)for(j=0;j<n;j++){ret.array[i][j]=0;for(k=0;k<n;k++)ret.array[i][j]=(ret.array[i][j]+(unit.array[i][k]*unit.array[k][j]));}for(i=0;i<n;i++)for(j=0;j<n;j++)unit.array[i][j]=ret.array[i][j];x>>=1;}for(ans=0,j=0;j<n;j++){ans+=pw.array[n-1][j]*num[j];}printf("%.0lfn",ans);}return 0;}
转载请注明:文章转载自 www.mshxw.com
本文地址:https://www.mshxw.com/it/374531.html
我们一直用心在做
关于我们 文章归档 网站地图 联系我们

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

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