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

;

void matrixmult(int a

,int b

){
int i,j,k;
int c

;
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

,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

;
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;
}