Description
给出一个二部图,左边
n
个点,标记分别为
Input
第一行输入一整数
(1≤T≤100,1≤n,m≤20,1≤q≤10,1≤ai,bi,ti≤109)
Output
对于每次查询,输出满足条件的子集数量
Sample Input
1
3 3 2
010
111
010
4 2 6
8 12 5
2 3
Sample Output
Case #1: 23 3
Solution
根据广义
Hall
引理,只需要从左边选取一个子集
X
使得其被一个覆盖盖住,从右边选取一个子集
Code
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
const int INF=0x3f3f3f3f,maxn=(1<<20)+5;
int T,n,m,q,a[22],b[22],L[maxn],R[maxn],dp[maxn];
char s[22];
int Solve(int *S,int n,int m,int A,int B)
{
int ans=0,N=1<<n;
for(int s=0;s<N;s++)
{
if((A&s)!=s)continue;
dp[s]=1;
int numa=0;
for(int i=0;i<n;i++)
if(s&(1<<i))
{
numa++,dp[s]&=dp[s^(1<<i)];
if(!dp[s])break;
}
if(!dp[s])continue;
int numb=0;
for(int i=0;i<m;i++)
if((B&(1<<i))&&(S[i]&s))numb++;
if(numb<numa)dp[s]=0;
if(dp[s])ans++;
}
return ans;
}
int main()
{
int Case=1;
scanf("%d",&T);
while(T--)
{
memset(L,0,sizeof(L));
memset(R,0,sizeof(R));
scanf("%d%d%d",&n,&m,&q);
for(int i=0;i<n;i++)
{
scanf("%s",s);
for(int j=0;j<m;j++)
if(s[j]=='1')L[i]|=(1<<j),R[j]|=(1<<i);
}
for(int i=0;i<n;i++)scanf("%d",&a[i]);
for(int i=0;i<m;i++)scanf("%d",&b[i]);
printf("Case #%d:",Case++);
while(q--)
{
int t,A=0,B=0;
scanf("%d",&t);
for(int i=0;i<n;i++)
if(a[i]%t==0)A|=(1<<i);
for(int i=0;i<m;i++)
if(b[i]%t==0)B|=(1<<i);
int ansl=Solve(R,n,m,A,B);
int ansr=Solve(L,m,n,B,A);
printf(" %I64d",(ll)ansl*ansr-1);
}
printf("\n");
}
return 0;
}