标签:scan ... mat 多个 ons ota sam tle 12px
Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)
Total Submission(s): 2990 Accepted Submission(s): 761
HDOJ 2008 Summer Exercise(4)- Buffet Dinner
代码:
1 /* 2 构造邻接矩阵A,f[i][j]表示i到j的道路,A^n中的各个元素f[i][j]表示从i到j走n步有多少种方法,因此就是矩阵乘法 3 A^t1,A^(t1+1)......A^t2;注意点数太大的用map离散化一下。 4 */ 5 #include<iostream> 6 #include<cstring> 7 #include<cstdio> 8 #include<map> 9 #include<cmath> 10 using namespace std; 11 const int mod=2008; 12 const int maxn=35; 13 map<string,int>m; 14 struct Lu 15 { 16 int f[maxn][maxn]; 17 }e,g[10004]; 18 Lu solve(Lu x,Lu y,int n) 19 { 20 Lu z; 21 memset(z.f,0,sizeof(z.f)); 22 for(int k=1;k<=n;k++) 23 { 24 for(int i=1;i<=n;i++) 25 { 26 if(!x.f[i][k]) continue; 27 for(int j=1;j<=n;j++) 28 { 29 z.f[i][j]+=x.f[i][k]*y.f[k][j]; 30 z.f[i][j]%=mod; 31 } 32 } 33 } 34 return z; 35 } 36 int main() 37 { 38 int n; 39 while(~scanf("%d",&n)) 40 { 41 m.clear(); 42 char a[30],b[30],v1[30],v2[30]; 43 int cnt=0; 44 memset(e.f,0,sizeof(e.f)); 45 for(int i=0;i<n;i++) 46 { 47 scanf("%s%s",&a,&b); 48 if(m[a]==0) m[a]=++cnt; 49 if(m[b]==0) m[b]=++cnt; 50 e.f[m[a]][m[b]]++; 51 } 52 memset(g[0].f,0,sizeof(g[0].f)); 53 for(int i=1;i<=cnt;i++) 54 g[0].f[i][i]=1; 55 for(int i=1;i<=10001;i++) 56 g[i]=solve(g[i-1],e,cnt); 57 int k,t1,t2; 58 scanf("%d",&k); 59 while(k--) 60 { 61 scanf("%s%s%d%d",v1,v2,&t1,&t2); 62 if(m[v1]==0||m[v2]==0||t2==0) 63 printf("0\n"); 64 else 65 { 66 // if(t1>t2) swap(t1,t2); 67 int ans=0; 68 for(int i=t1;i<=t2;i++){ 69 if(i==0) continue; 70 ans=(ans+g[i].f[m[v1]][m[v2]])%mod; 71 } 72 printf("%d\n",ans); 73 } 74 } 75 } 76 return 0; 77 }
标签:scan ... mat 多个 ons ota sam tle 12px
原文地址:http://www.cnblogs.com/--ZHIYUAN/p/6193208.html