矩阵取数游戏c语言,NOIP2002pj产生数[floyd 高精度]

noip2002普及组第三题

－－－－－－－－－－－－－－－－

一个数可以变换多次，floyd求传递闭包(初始化d[i][i]=1)，乘法原理更新答案

要用高精度，注意输出

//

// main.cpp

// noip2002产生数

//

// Created by Candy on 9/10/16.

// Copyright © 2016 Candy. All rights reserved.

//

#include

#include

#include

#include

using namespace std;

typedef unsigned long long ll;

const int N=,B=1e4;

inline int read(){

int x=,f=;char ch=getchar();

while(ch''){if(ch=='-')f=-;ch=getchar();}

while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}

return x*f;

}

int k,x,y,d[N][N],f[N];

char s[N];

void floyd(){

for(int i=;i<=;i++) d[i][i]=;

for(int k=;k<=;k++)

for(int i=;i<=;i++)

for(int j=;j<=;j++)

d[i][j]=d[i][j]||(d[i][k]&&d[k][j]);

for(int i=;i<=;i++)

for(int j=;j<=;j++) if(d[i][j]) f[i]++;

}

struct big{

int d[],size;

big(){size=;}

} ans;

void chengInt(big &a,int k){

int g=,i;

for(i=;i<=a.size;i++){

int tmp=a.d[i]*k;

a.d[i]=(tmp+g)%B;

g=(tmp+g)/B;

}

while(g){

a.d[i++]=g%B; a.size++;

g/=B;

}

}

int main(int argc, const char * argv[]) {

scanf("%s%d",s,&k);

for(int i=;i<=k;i++) scanf("%d%d",&x,&y),d[x][y]=;

floyd();

ans.d[]=;

int len=strlen(s);

for(int i=;i

int a=s[i]-'';

chengInt(ans,f[a]);

//printf("f %d %d\n",a,f[a]);

}

for(int i=ans.size;i>=;i--){

if(i!=ans.size){

if(ans.d[i]

else if(ans.d[i]

else if(ans.d[i]

}

cout<

}

return ;

}