2016.7.13 刷题打卡

到22:29 才写完一道题：
bzoj 1898
主要是 搞了搞矩阵乘法在 信息学中的应用 没错就是 国家集训队论文
另外翻出来 lyd讲的 沼泽鳄鱼，写了些：

这是一道 dp题 但是 光看k就知道这个要 优化，怎么优化呢？

首先我们考虑 没有 鳄鱼的情况那问题就变成 了 从s到e路径为k 的 路径。
但是现在 我们有一群烦人的 小生物 鳄鱼。
没关系，不虚他！ 我们只要建立以下 Gi 即删掉 i时刻的不能经过的 点的相关的边。

那么此问题变成： π（i=1-k） Gi

那我如何快速的乘出这些东西呢？

这个时候题目里的重要信息 T=2，3,4 就要用上了！

这些的最小公倍数 是 12
所以 我们有 字丑勿怪。。。

那么上代码：

#include<cstdio>  
#include<cstdlib>  
#include<cstring>  
// by mar_ch
#define maxn 60  
#define mod 10000  

typedef int pp[maxn][maxn];  

int n,m,start,end,k,Nfish,w[21][5],ans=0,len[maxn];  
pp t,f[12],guy,all,a,apple;  

//矩乘；   
inline void mul(pp x,pp &y)  
{  
    memset(t,0,sizeof(t));  
    for(int i=0;i<n;i++)  
     for(int j=0;j<n;j++)  
      for(int k=0;k<n;k++)  
       t[i][j]=(t[i][j]+x[i][k]*y[k][j]%mod)%mod;  
    memcpy(y,t,sizeof(y));  
}  

inline void fast(int y)  
{   
    for(;y;y>>=1,mul(all,all))  
    if(y&1)mul(all,guy);  
}  

int main()  
{  
    scanf("%d%d%d%d%d",&n,&m,&start,&end,&k);  
    for(int i=1,x,y;i<=m;i++){  
        scanf("%d%d",&x,&y);  
        a[x][y]=a[y][x]=1;  
    }  
    scanf("%d",&Nfish);  
    for(int i=1,t;i<=Nfish;i++){  
        scanf("%d",&len[i]);  
        for(int j=0;j<len[i];j++)  
        scanf("%d",&w[i][j]);  
    }  
    for(int i=0;i<n;i++)all[i][i]=1;  

    for(int i=0;i<12;i++){  
        memcpy(f[i],a,sizeof(f[i]));  
        for(int j=1;j<=Nfish;j++){  
            for(int k=0;k<n;k++)  
             f[i][w[j][(i+1)%len[j]]][k]=f[i][k][w[j][i%len[j]]]=0;  
        }  
        mul(f[i],all);  
    }  
    for(int i=0;i<n;i++)guy[i][i]=1;  
    fast(k/12);  
    for(int i=0;i<k%12;i++)mul(f[i],guy);  
    printf("%d\n",guy[end][start]);  
    return 0;  
}  

嗯嗯。

然后 容我洗个澡，回来再大战 poj 3613

嗯哼，题解明天。。

#include<cstdio>  
#include<cstring>  
#include<climits>  
#define min(a,b) ((a)<(b))?(a):(b)  

const int N = 205;  
const int INF = 1000000100;  

int m,p[1001];  
int ans[N][N];  

void matrixmult(int a[N][N],int b[N][N]){  
    int i,j,k;  
    int c[N][N];  
    for(i=0;i<m;i++){  
        for(j=0;j<m;j++){  
            c[i][j]=INF;  
            for(k=0;k<m;k++)  
                c[i][j]=min( a[i][k]+b[k][j] , c[i][j] );  
        }  
    }  
    for(i=0;i<m;i++){  
        for(j=0;j<m;j++)a[i][j]=c[i][j];  
    }  
}  

void mult(int s[N][N],int k){  
    int i,j;  
    for(i=0;i<m;i++){  
        for(j=0;j<m;j++)  
            ans[i][j]=s[i][j];  
    }  
    k--;  
    while(k>0){  
        if(k&1) matrixmult(ans,s);  
        k=k>>1;  
        matrixmult(s,s);  
    }  
}  

int main()  
{  
    int mat[N][N];  
    int i,j,n,t,s,e,from,to,val;  
    while(~scanf("%d%d%d%d",&n,&t,&s,&e))  
    {  
        for(i=0;i<N;i++){  
            for(j=0;j<N;j++)mat[i][j]=INF;  
        }  
        memset(p,-1,sizeof(p));  
        m=0;  
        for(i=0;i<t;i++){  
            scanf("%d%d%d",&val,&from,&to);  
            if(p[from]<0) p[from] = m++;  
            if(p[to]<0) p[to] = m++;  
            if(val<mat[ p[from] ][ p[to] ])  
                mat[ p[from] ][ p[to] ] = mat[ p[to] ][ p[from] ] = val;  
        }  
        s=p[s]; e=p[e];  
        mult(mat,n);  
        printf("%d\n",ans[s][e]);  
    }  
    return 0;  
}