NKOJ 1901 喜欢奇数的面包师
时间限制 : 20000 MS 空间限制 : 65536 KB
问题描述
有n(n<=100)个面包师,每个面包师门上有粉笔画的横线,每个面包师有若干个喜欢的的面包师。
每个月门上粉笔横线数目为奇数的面包师将获奖,他获奖后会在他喜欢的面包师的门上加上一条粉笔横线,问第t(t<=1,000,000,000)个月获奖的面包师有多少人?
输入格式
第一行,一个整数k,表示接下来有k组测试数据
对于每组测试数据,格式如下:
第一行两个整数n和t
接下来n行,每行描述一个面包师,每行首先是面包师的名字(不超过20个小写字母构成);接着是一个整数,表示该面包师的门上已经画上的横线数目;接下来是一个整数,表示该面包师喜欢的其他面包师的数目,最后是该面包师喜欢的其他面包师的名字。
输出格式
对于每组测试数据,输出一个整数,表示第t个月获奖的面包师人数。
样例输入
2
3 2
bessie 2 3 bessie linda mary
mary 1 1 linda
linda 0 1 bessie
2 2
siavosh 1 2 siavosh mohammad
mohammad 1 0
样例输出
2
0
来源 Tehran Sharif
思路:
用一维数组f1表示第一个月各个面包师的得分,f[i]为1表示该面包师得分为奇数,为0表示偶数
用二维数组a表示各个面包师的“复杂关系”a[i][j]表示i喜欢j
然后发现f1*a后所得的g数组表示第二个月各面包师增加得分的增加量,(g+f1)%2即为第二个月的f2,即f(i)=f(i-1) *a+f(i-1)
即 f(i)=f(i-1) * (a+E),E为单位矩阵,迭代后得f(i)=f1*(a+E)^(i-1);
优化:因为该题中矩阵运算全是模2意义下,所以用&代替*,^代替+。
代码:
#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
const int need=103;
typedef int int_[need][need];
string name[need];
int n,t,f[need];
int_ a,ans,c;
string bb;int tot=0;
int get_id()
{
cin>>bb;
for(int i=1;i<=tot;i++)
if(bb==name[i]) return i;
name[++tot]=bb;
return tot;
}
void matrix_multi(int_ a,int_ b)
{
memset(c,0,sizeof(c));
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
for(int k=1;k<=n;k++)
c[i][j]^=a[i][k]&b[k][j];
memcpy(a,c,sizeof(c));
}
void matrix_power(int b)
{
for(int i=1;i<=n;i++) ans[i][i]=1;
while(b)
{
if(b&1) matrix_multi(ans,a);
b>>=1;
matrix_multi(a,a);
}
}
void clean()
{
tot=0;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++) ans[i][j]=a[i][j]=0;
f[i]=0;
name[i]="";
}
}
int main()
{
int k;scanf("%d",&k);
while(k)
{
k--;
scanf("%d%d",&n,&t);
for(int i=1,b,p,q;i<=n;i++)
{
b=get_id();
scanf("%d%d",&p,&q);
f[b]=p&1;
for(int j=1;j<=q;j++)a[b][get_id()]=1;
}
for(int i=1;i<=n;i++) a[i][i]^=1;
matrix_power(t-1);
int cnt=0;
for(int i=1,now;i<=n;i++)
{
now=0;
for(int j=1;j<=n;j++) now^=f[j]&ans[j][i];
if(now) cnt++;
}
printf("%d\n",cnt);
clean();
}
}