## CHN_JZ的博客

while(!die) is_a_oier();

# [二维树状数组]BZOJ 1452—— [JSOI2009]Count

## 解题思路

#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=25,maxm=125,tt=45989;
struct jz{
int x[maxm][maxm],n,m;
jz operator*(const jz &b)const{
jz c;memset(c.x,0,sizeof(c.x));
c.n=n,c.m=b.m;
for (int i=1;i<=c.n;i++)
for (int j=1;j<=c.m;j++)
for (int k=1;k<=b.n;k++)
c.x[i][j]=(c.x[i][j]+x[i][k]*b.x[k][j])%tt;
return c;
}
}a,b;
int n,m,t,A,B,tot=1,nxt[maxm],son[maxm],lnk[maxn],ans;
nxt[++tot]=lnk[x];lnk[x]=tot;son[tot]=y;
}
void work(){
for (int i=0;i<n;i++)
for (int j=lnk[i];j;j=nxt[j])
for (int k=lnk[son[j]];k;k=nxt[k])
if (k!=(j^1)) b.x[j][k]++;
}
jz qsm(jz w,int b){
jz c;c.n=w.n;c.m=w.m;
for (int i=1;i<=c.n;i++)
for (int j=1;j<=c.m;j++) c.x[i][j]=(i==j);
while (b>0){
if (b%2==1) c=c*w;
w=w*w;
b=b>>1;
}
return c;
}
int main(){
freopen("exam.in","r",stdin);
freopen("exam.out","w",stdout);
scanf("%d%d%d%d%d",&n,&m,&t,&A,&B);
for (int i=1;i<=m;i++){
int x,y;scanf("%d%d",&x,&y);
}
work();
for (int j=lnk[A];j;j=nxt[j]) a.x[1][j]++;
a.n=1;b.n=a.m=b.m=tot;
a=a*qsm(b,t-1);
for (int j=lnk[B];j;j=nxt[j]) ans=(ans+a.x[1][j^1])%tt;
printf("%d\n",ans);
}

#### BZOJ1452: [JSOI2009]Count 二维树状数组

2016-10-13 19:38:10

#### 【二维树状数组】BZOJ1452 [JSOI2009]Count

2017-09-29 20:02:38

#### bzoj1452: [JSOI2009]Count(二维树状数组)

2017-09-13 22:21:29

#### BZOJ 1452: [JSOI2009]Count 二维树状数组

2017-05-02 21:03:50

#### 二维树状数组（bzoj 1452: [JSOI2009]Count）

2017-08-07 23:44:29

#### BZOJ 1452 [JSOI2009]Count 二维树状数组

2015-03-19 23:07:08

#### BZOJ 1452: [JSOI2009]Count （二维树状数组）

2014-08-21 23:49:07

#### 【bzoj 1452】[JSOI2009]Count（二维树状数组）

2016-10-09 20:09:01

#### 1452: [JSOI2009]Count 二维树状数组

2017-07-28 10:54:51

#### 【BZOJ】1452[JSOI2009]Count

2018-03-16 21:26:13