bzoj-3444 3444: 最后的晚餐(组合数学)

题目链接:

3444: 最后的晚餐

Time Limit: 5 Sec  Memory Limit: 128 MB

Description

【问题背景】

高三的学长们就要离开学校,各奔东西了。某班n人在举行最后的离别晚餐时,饭店老板觉得十分纠结。因为有m名学生偷偷找他,要求和自己暗恋的同学坐在一起。

【问题描述】

饭店给这些同学提供了一个很长的桌子,除了两头的同学,每一个同学都与两个同学相邻(即坐成一排)。给出所有信息,满足所有人的要求,求安排的方案总数(这个数字可能很大,请输出方案总数取余989381的值,也可能为0)。

Input

输入有m+1行,第一行有两个用空格隔开的正整数n、m,如题所示。接下来的m行,每一行有两个用空格隔开的正整数,第i行为Ai和Bi,表示Ai的暗恋对象为Bi,保证Ai互不相等。

Output

输出只有一行,这一行只有一个数字,如题所示。

Sample Input

4 2
1 2
4 3

Sample Output

8

题意:

思路:

只有单链和点才满足要求,分叉和环都不能有;链有两个方向放所以ans=(l+r)!*2^l;l为链的个数,r为点的个数;

AC代码:

/**************************************************************
    Problem: 3444
    User: LittlePointer
    Language: C++
    Result: Accepted
    Time:6100 ms
    Memory:67784 kb
****************************************************************/

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <bits/stdc++.h>
#include <stack>
#include <map>

using namespace std;

#define For(i,j,n) for(int i=j;i<=n;i++)
#define mst(ss,b) memset(ss,b,sizeof(ss));

typedef  long long LL;

template<class T> void read(T&num) {
    char CH; bool F=false;
    for(CH=getchar();CH<‘0‘||CH>‘9‘;F= CH==‘-‘,CH=getchar());
    for(num=0;CH>=‘0‘&&CH<=‘9‘;num=num*10+CH-‘0‘,CH=getchar());
    F && (num=-num);
}
int stk[70], tp;
template<class T> inline void print(T p) {
    if(!p) { puts("0"); return; }
    while(p) stk[++ tp] = p%10, p/=10;
    while(tp) putchar(stk[tp--] + ‘0‘);
    putchar(‘\n‘);
}

const LL mod=989381;
const double PI=acos(-1.0);
const int inf=1e9;
const int N=5e5+20;
const int maxn=1e4+220;
const double eps=1e-12;

int n,m,in[N],p[N],sum[N],num,vis[N];
vector<int>ve[N];
int findset(int x)
{
    if(x==p[x])return x;
    return p[x]=findset(p[x]);
}
inline void same(int x,int y)
{
    int fx=findset(x),fy=findset(y);
    if(fx!=fy)
    {
        p[fx]=fy;
        sum[fy]+=sum[fx];
        sum[fx]=0;
    }
}

map<int,int>mp[N];
void dfs(int cur,int fa)
{
    vis[cur]=1;
    num++;
    int len=ve[cur].size();
    for(int i=0;i<len;i++)
    {
        int temp=ve[cur][i];
        if(temp==fa)continue;
        dfs(temp,cur);
    }
}
int check()
{
    for(int i=1;i<=n;i++)
    {
        if(in[i]>2)return 0;
    }
    num=0;
    for(int i=1;i<=n;i++)
    {
        if(!vis[i])
        {
            if(in[i]==0)vis[i]=1,num++;
            else if(in[i]==1)dfs(i,0);
        }
    }
    if(num!=n)return 0;
    return 1;
}
LL pow_mod(int y)
{
    LL s=1,base=2;
    while(y)
    {
        if(y&1)s=s*base%mod;
        base=base*base%mod;
        y>>=1;
    }
    return s;
}
int main()
{
    //freopen("in.txt","r",stdin);
    read(n);read(m);
    for(int i=0;i<=n;i++)p[i]=i,sum[i]=1;
    int u,v;
    For(i,1,m)
    {
        read(u);read(v);
        if(mp[v][u])continue;
        ve[u].push_back(v);
        ve[v].push_back(u);
        same(u,v);
        in[u]++;
        in[v]++;
        mp[u][v]=1;
    }
    LL ans=0;
    if(check())
    {
        ans=1;
        int l=0,r=0;
        for(int i=1;i<=n;i++)
        {
            if(p[i]!=i)continue;
            if(sum[i]>1)l++;
            else r++;
        }
        for(int i=1;i<=l+r;i++)ans=ans*(LL)i%mod;
        ans=ans*pow_mod(l)%mod;
    }
    cout<<ans<<endl;

    return 0;
}

  

时间: 2024-10-24 15:13:40

bzoj-3444 3444: 最后的晚餐(组合数学)的相关文章

BZOJ 1211: [HNOI2004]树的计数( 组合数学 )

知道prufer序列就能写...就是求个可重集的排列...先判掉奇怪的情况, 然后答案是(N-2)!/π(d[i]-1)! --------------------------------------------------------------------------- #include<cstdio> #include<algorithm> #include<cstring> using namespace std; typedef long long ll; c

BZOJ 2111 ZJOI2010 Perm 排列计数 组合数学+Lucas定理

题目大意:求1~n的排列能组成多少种小根堆 考虑一个1~i的排列所构成的堆,l为左儿子大小,r为右儿子的大小 那么1一定是堆顶 左儿子和右儿子分别是一个堆 显然如果选出l个数给左儿子 那么左儿子的方案数显然是f[l],右儿子的方案数为f[r] 于是有f[i]=C(i-1,l)*f[l]*f[r] 于是我们线性筛处理出阶乘和阶乘的逆元 代入即可得到WA 原因是这题n可以大于p 此时要用到Lucas定理 坑死了 #include <cstdio> #include <cstring>

BZOJ 2467 中山市选2010 生成树 组合数学

题目大意:给定一个图,图的中心是一个n个点的多边形,每条边都外接一个五边形,求生成树个数 Matrix Tree定理?不会! 观察这个图 5n条边 4n个点 每个五边形都是一个环 必须拆掉一条边 拆掉之后发现4n个点 4n条边 是一个基环树 基环树的环上的边由中心多边形被拆掉的边所在的五边形的剩余边与中心多边形未被拆掉的边构成 容易发现这个环上任意拆掉一条边都会导致某个五边形被拆掉两条边 且一条边在中心多边形上 于是可知 这个图成为一棵树当且仅当一个五边形被拆掉两条边 剩余五边形被拆掉一条边 且

bzoj 2142: 礼物【中国剩余定理+组合数学】

参考:http://blog.csdn.net/wzq_qwq/article/details/46709471 首先推组合数,设sum为每个人礼物数的和,那么答案为 \[ ( C_{n}^{sum}C_{sum}^{w[1]}c_{sum-w[1]}^{w[2]}... \] 设w[0]=n-sum,然后化简成阶乘的形式: \[ \frac{n!}{w[0]!w[1]!...w[n]!} \] 注意到这里p不是质数,所以把p拆成质数的方相乘的形式,最后用中国剩余定理合并即可 然后现在的问题是怎

bzoj 3027: [Ceoi2004]Sweet【生成函数+组合数学】

首先根据生成函数的套路,这个可以写成: \[ \prod_{i=1}^{n}(1+x^1+x^2+...+x^{c[i]}) \] 然后化简 \[ =\prod_{i=1}^{n}\frac{1-x^{c[i]+1}}{1-x} \] \[ =\prod_{i=1}^{n}\frac{1}{1-x}*(1-x^{c[i]+1}) \] \[ =(1+x^1+x^2+...)^n*\prod_{i=1}^{n}(1-x^{c[i]+1}) \] 位数过多所以只考虑有常数项的位,后面那个式子可以df

BZOJ 3444: 最后的晚餐( )

把暗恋关系看成无向边, 那某个点度数超过2就无解.存在环也是无解.有解的话对连通分量进行排列就行了. ---------------------------------------------------------------------------------- #include<cstdio> #include<algorithm> #include<cstring> using namespace std; typedef long long ll; const

poj 3444 Wavelet Compression 模拟水题

水题,直接贴代码. //poj 3444 //sep9 #include <iostream> using namespace std; const int maxN=260; int a[maxN],b[maxN]; int main() { int i,n,m; while(scanf("%d",&n)==1&&n){ for(i=1;i<=n;++i) scanf("%d",&a[i]); m=1; while

BZOJ 1008 越狱 (组合数学)

题解:正难则反,从总数中减去全部相邻不相同的数目就是答案,n*(n-1)^(m-1):第一个房间有n中染色方案,剩下m-1个房间均只有n-1种染色方案,用总数减就是答案. #include <cstdio> const int mod=100003; typedef long long LL; LL n,m; LL power(LL a,LL b){ LL ans=1; a%=mod; while(b){ if(b&1)ans=ans*a%mod; a=a*a%mod,b>>

BZOJ 1005 明明的烦恼 (组合数学)

题解:n为树的节点数,d[ ]为各节点的度数,m为无限制度数的节点数. 则             所以要求在n-2大小的数组中插入tot各序号,共有种插法: 在tot各序号排列中,插第一个节点的方法有种插法: 插第二个节点的方法有种插法: ......... 另外还有m各节点无度数限制,所以它们可任意排列在剩余的n-2-tot的空间中,排列方法总数为 根据乘法原理: #include <cstdio> #include <cmath> int n,m,tot,i,j,d,down