标签：幻想 include int while 解题 mul SHOI2016 mod getchar

「SHOI2016」黑暗前的幻想乡

sb题想不出来，应该去思考原因，而不是自暴自弃

一开始总是想着对子树做dp，但是状态压不起去，考虑用容斥消减一些条件变得好统计，结果越想越乱。

期间想过矩阵树定理，但没想清楚又被我忽略了。

其实非常简单

你对着所有的东西跑一遍生成树计数，然后你发现统计了同一个施工队的方案，然后发现可以枚举子集，就是个sb容斥了

Code:

#include <cstdio>
#include <cctype>
#include <algorithm>
const int mod=1e9+7;
const int N=20;
inline int add(int a,int b){return a+b>=mod?a+b-mod:a+b;}
#define mul(a,b) (1ll*(a)*(b)%mod)
int qp(int d,int k){int f=1;while(k){if(k&1)f=mul(f,d);d=mul(d,d),k>>=1;}return f;}
template <class T>
void read(T &x)
{
    x=0;char c=getchar();
    while(!isdigit(c)) c=getchar();
    while(isdigit(c)) x=x*10+c-'0',c=getchar();
}
struct node{int u,v;}E[N][N*N];
int n,ans,g[N][N],cnt[N];
int Gauss()
{
    int f=0,ret=1;
    for(int i=1;i<=n;i++)
    {
        int id=i;
        for(int j=i;j<=n;j++)
            if(g[j][i])
            {
                id=j;
                break;
            }
        if(g[id][i]==0) return 0;
        if(id!=i) std::swap(g[i],g[id]),f^=1;
        for(int j=i+1;j<=n;j++)
        {
            int inv=mul(g[j][i],qp(g[i][i],mod-2));
            for(int k=i;k<=n;k++)
                g[j][k]=add(g[j][k],mod-mul(g[i][k],inv));
        }
        ret=mul(ret,g[i][i]);
    }
    return f?add(mod,-ret):ret;
}
int main()
{
    read(n);
    --n;
    for(int i=1;i<=n;i++)
    {
        read(cnt[i]);
        for(int j=1;j<=cnt[i];j++)
            read(E[i][j].u),read(E[i][j].v);
    }
    for(int s=0;s<1<<n;s++)
    {
        int ct=0;
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                g[i][j]=0;
        for(int i=1;i<=n;i++)
            if(s>>i-1&1)
            {
                for(int j=1;j<=cnt[i];j++)
                {
                    int u=E[i][j].u,v=E[i][j].v;
                    g[u][v]=add(g[u][v],mod-1),g[v][u]=add(g[v][u],mod-1);
                    ++g[u][u],++g[v][v];
                }
                ++ct;
            }
        ct=(n-ct&1)?mod-1:1;
        ans=add(ans,mul(ct,Gauss()));
    }
    printf("%d\n",ans);
    return 0;
}

2019.3.12

标签：幻想,include,int,while,解题,mul,SHOI2016,mod,getchar
来源： https://www.cnblogs.com/butterflydew/p/10517161.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。