51nod1806：wangyurzee的树（生成树计数 & 容斥）

最新推荐文章于 2021-05-18 18:39:07 发布
最新推荐文章于 2021-05-18 18:39:07 发布
阅读量270 收藏
点赞数
分类专栏： 数论
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/junior19/article/details/79815991
版权

1806 wangyurzee的树 
基准时间限制：1 秒 空间限制：131072 KB 分值: 80  难度：5级算法题
 收藏
 关注
wangyurzee有n个各不相同的节点，编号从1到n。wangyurzee想在它们之间连n-1条边，从而使它们成为一棵树。
可是wangyurzee发现方案数太多了，于是他又给出了m个限制条件，其中第i个限制条件限制了编号为u[i]的节点的度数不能为d[i]。
一个节点的度数，就是指和该节点相关联的边的条数。
这样一来，方案数就减少了，问题也就变得容易了，现在请你告诉wangyurzee连边的方案总数为多少。
答案请对1000000007取模。

样例解释
总方案共有3种，分别为{(1,2),(1,3)}，{(1,2),(2,3)}，{(2,3),(1,3)}。其中第二种方案节点1的度数为2，不符合要求，因此答案为2。
Input 
第一行输入2个整数n(1<=n<=1000000)，m(0<=m<=17)分别表示节点个数以及限制个数。
第2行到第m+1行描述m个限制条件，第i+1行为2个整数u[i]，d[i]，表示编号为u[i]的节点度数不能为d[i]。
为了方便起见，保证1<=ui<=m。同时保证1<=ui<=n，1<=di<=n-1，保证不会有两条完全相同的限制。
Output 
输出一行一个整数表示答案。
Input示例 
3 1
1 2
Output示例 
2

思路：根据Cayley公式，带标号无根生成树是pow(n,n-2)种；根据Prufer公式，每种无根树对应n-2长度的编码，序列中某编号出现的次数为度数-1，因此这题容斥一下就行，特判n=1。

# include <bits/stdc++.h>
# define A first
# define B second
using namespace std;
typedef long long LL;
const int maxn = 1e6;
const LL mod = 1e9+7;
pair<int,int>a[20];
LL fac[maxn+3]={1,1}, inv[maxn+3]={1,1}, fi[maxn+3]={1,1};
int vis[maxn+3];
LL qmod(LL x, LL y)
{
    LL res = 1;
    for(;y;y>>=1)
    {
        if(y&1) res=res*x%mod;
        x=x*x%mod;
    }
    return res;
}
int main()
{
    for(int i=2; i<=maxn; ++i)
    {
        fac[i]=fac[i-1]*i%mod;
        inv[i]=(mod-mod/i)*inv[mod%i]%mod;
        fi[i]=fi[i-1]*inv[i]%mod;
    }
    LL n, ans=0;
    int m;
    scanf("%lld%d",&n,&m);
    if(n==1) return 0*puts("1");
    for(int i=0,x,y; i<m; ++i)
    {
        scanf("%d%d",&x,&y);
        a[i]=make_pair(x,y);
    }
    sort(a,a+m);
    for(int i=0; i<1<<m; ++i)
    {

        int cnt=0, sum=0, flag=1, pre=-1;
        LL tmp = fac[n-2];
        for(int j=0; j<m; ++j)
        {
            if(i&(1<<j))
            {
                if(~pre && pre==a[j].A)
                {
                    flag = 0;
                    break;
                }
                pre = a[j].A;
                tmp = tmp*fi[a[j].B-1]%mod;
                ++cnt;
                sum += a[j].B-1;
            }
        }
        if(sum > n-2 || !flag) continue;
        tmp = tmp*fi[n-2-sum]%mod;
        tmp = tmp*qmod(n-cnt,n-2-sum)%mod;
        if(cnt&1) ans=(ans-tmp+mod)%mod;
        else ans=(ans+tmp)%mod;
    }
    printf("%lld\n",ans);
    return 0;
}


#include "Winsock.h" 
#include "windows.h" 
#include "stdio.h"  
#define RECV_PORT 3312  
#define SEND_PORT 4302  
#pragma   comment(lib, "wsock32.lib")
SOCKET sockclient, sockserver;
 struct sockaddr_in ServerAddr;//服务器地址
 struct sockaddr_in ClientAddr;//客户端地址 

/***********************全局变量***********************/
int Addrlen;//地址长度 
char filename[20];//文件名 
char order[10];//命令  
char rbuff[1024];//接收缓冲区  
char sbuff[1024];//发送缓冲区  


DWORD StartSock()    //初始化winsock   
{
    WSADATA WSAData;
    if (WSAStartup(MAKEWORD(2, 2), &WSAData) != 0)
    {
        printf("socket init fail!\n");
        return (-1);
    }
    return(1);
}

DWORD CreateSocket()
{
    sockclient = socket(AF_INET, SOCK_STREAM, 0);
    if (sockclient == SOCKET_ERROR)
    {
        printf("sockclient create fail ! \n");
        WSACleanup();
        return(-1);
    }
    ServerAddr.sin_family = AF_INET;
    ServerAddr.sin_addr.s_addr = htonl(INADDR_ANY);
    ServerAddr.sin_port = htons(RECV_PORT);
    if (bind(sockclient, (struct  sockaddr  FAR  *)&ServerAddr, sizeof(ServerAddr)) == SOCKET_ERROR)
    {                     //bind函数将套接字和地址结构绑定   
        printf("bind is the error");
        return(-1);
    }
    return (1);
}

int SendFileRecord(SOCKET datatcps, WIN32_FIND_DATA *pfd)     //用来发送当前文件记录 
{
    char filerecord[MAX_PATH + 32];
    FILETIME ft;         //文件建立时间   
    FileTimeToLocalFileTime(&pfd->ftLastWriteTime, &ft);
    SYSTEMTIME lastwtime;     //SYSTEMTIME系统时间数据结构   
    FileTimeToSystemTime(&ft, &lastwtime);
    char *dir = pfd->dwFileAttributes&FILE_ATTRIBUTE_DIRECTORY ? "<DIR>" : " ";
    sprintf(filerecord, "%04d-%02d-%02d %02d:%02d  %5s %10d   %-20s\n",
        lastwtime.wYear,
        lastwtime.wMonth,
        lastwtime.wDay,
        lastwtime.wHour,
        lastwtime.wMinute,
        dir,
        pfd->nFileSizeLow,
        pfd->cFileName);
    if (send(datatcps, filerecord, strlen(filerecord), 0) == SOCKET_ERROR)
    { //通过datatcps接口发送filerecord数据，成功返回发送的字节数   
        printf("Error occurs when sending file list!\n");
        return 0;
    }
    return 1;
}

int SendFileList(SOCKET datatcps)
{
    HANDLE hff;//建立一个线程  
    WIN32_FIND_DATA fd;   //搜索文件   
//可以通过FindFirstFile（）函数根据当前的文件存放路径查找该文件来把待操作文件的相关属性读取到WIN32_FIND_DATA结构中去  
    if (hff == INVALID_HANDLE_VALUE)//发生错误  
    {
        const char *errstr = "can't list files!\n";
        printf("list file error!\n");
        if (send(datatcps, errstr, strlen(errstr), 0) == SOCKET_ERROR)
        {
            printf("error occurs when senging file list!\n");
        }
        closesocket(datatcps);
        return 0;
    }

    BOOL fMoreFiles = TRUE;
    while (fMoreFiles)
    {//发送此项文件信息   
        if (!SendFileRecord(datatcps, &fd))
        {
            closesocket(datatcps);
            return 0;
        }
        //搜索下一个文件    
        fMoreFiles = FindNextFile(hff, &fd);
    }
    closesocket(datatcps);
    return 1;
}

int SendFile(SOCKET datatcps, FILE* file)
{
    printf(" sending file data..");
    for (;;)   //从文件中循环读取数据并发送客户端   
    {
        int r = fread(sbuff, 1, 1024, file);//把file里面的内容读到sbuff缓冲区   
        if (send(datatcps, sbuff, r, 0) == SOCKET_ERROR)
        {
            printf("lost the connection to client!\n");
            closesocket(datatcps);
            return 0;
        }
        if (r<1024)//文件传送结束    
            break;
    }
    closesocket(datatcps);
    printf("done\n");
    return 1;
}

//连接  
DWORD ConnectProcess()
{
    Addrlen = sizeof(ClientAddr);
    if (listen(sockclient, 5)<0)
    {
        printf("Listen error");
        return(-1);
    }
    printf("服务器监听中...\n");
    for (;;)
    {
        sockserver = accept(sockclient, (struct sockaddr FAR *)&ClientAddr, &Addrlen);
        //accept函数取出连接队列的第一个连接请求，sockclient是处于监听的套接字ClientAddr 是监听的对象地址，         
        //Addrlen是对象地址的长度 
        for (;;)
        {
            memset(rbuff, 0, 1024);
            memset(sbuff, 0, 1024);
            if (recv(sockserver, rbuff, 1024, 0) <= 0)
            {
                break;
            }
            printf("\n");
            printf("获取并执行的命令为：");
            printf(rbuff);
            if (strncmp(rbuff, "get", 3) == 0)
            {
                strcpy(filename, rbuff + 4); printf(filename);
                FILE *file; //定义一个文件访问指针    
                //处理下载文件请求    
                file = fopen(filename, "rb");//打开下载的文件，只允许读写   
                if (file)
                {
                    sprintf(sbuff, "get file %s\n", filename);
                    if (!send(sockserver, sbuff, 1024, 0))
                    {
                        fclose(file);      return 0;
                    }
                    else
                    {//创建额外数据连接传送数据     
                        if (!SendFile(sockserver, file))
                            return 0;
                        fclose(file);
                    }
                }//file   

                else//打开文件失败    
                {
                    strcpy(sbuff, "can't open file!\n");
                    if (send(sockserver, sbuff, 1024, 0))
                        return 0;
                } //lost 
            }//get

            if (strncmp(rbuff, "put", 3) == 0)
            {
                FILE *fd;
                int count;
                strcpy(filename, rbuff + 4);
                fd = fopen(filename, "wb");
                if (fd == NULL)
                {
                    printf("open file %s for weite failed!\n", filename);
                    return 0;
                }
                sprintf(sbuff, "put file %s", filename);
                if (!send(sockserver, sbuff, 1024, 0))
                {
                    fclose(fd);
                    return 0;
                }
                while ((count = recv(sockserver, rbuff, 1024, 0))>0)//recv函数返回接受的字节数赋给count          
                    fwrite(rbuff, sizeof(char), count, fd);
                //把count个数据长度为size0f（）的数据从rbuff输入到fd指向的目标文件
                printf(" get %s succed!\n", filename);
                fclose(fd);
            }//put 

            if (strncmp(rbuff, "pwd", 3) == 0){
                char   path[1000];
                GetCurrentDirectory(1000, path);//找到当前进程的当前目录  
                strcpy(sbuff, path);
                send(sockserver, sbuff, 1024, 0);
            }//pwd   

            if (strncmp(rbuff, "dir", 3) == 0){
                strcpy(sbuff, rbuff);
                send(sockserver, sbuff, 1024, 0);
                SendFileList(sockserver);//发送当前列表 
            }//dir 

            if (strncmp(rbuff, "cd", 2) == 0)
            {
                strcpy(filename, rbuff + 3);
                strcpy(sbuff, rbuff);
                send(sockserver, sbuff, 1024, 0);
                SetCurrentDirectory(filename);//设置当前目录 
            }//cd  
            closesocket(sockserver);
        }//for 2
    }//for 1 
}


int main()
{
    if (StartSock() == -1)
        return(-1);
    if (CreateSocket() == -1)
        return(-1);
    if (ConnectProcess() == -1)
        return(-1);
    return(1);
}

junior19
关注
专栏目录
51Nod1766: 树上的最远点对 题解
IcePrincess_1968's World
07-11 305
这题用到了树的直径的非常好的性质，有必要记录 树的直径满足这样一个性质： 设一棵树的一个点集SSS的直径的端点集合是{a1,b1a1,b1a_1,b_1}(如果有多条直径选取任意一条)，另一个点集TTT的直径的端点集合是{a2,b2a2,b2a_2,b_2}，则点集S⋃TS⋃TS\bigcup T的直径的端点集合{a3,b3a3,b3a_3,b_3}∈∈\in{a1,b1,a2,b2a1,b1...
51nod 1352：集合计数
07-05 1530
1352 集合计数 基准时间限制：1 秒 空间限制：131072 KB 分值: 20 难度：3级算法题  收藏  关注 给出N个固定集合{1，N},{2,N-1},{3,N-2},...,{N-1,2},{N,1}.求出有多少个集合满足：第一个元素是A的倍数且第二个元素是B的倍数。 提示： 对于第二组测试数据，集合分别是：{1,1
51nod 1601 完全图的最小生成树计数 字典树+最小生成树
beginend
03-06 542
题意 给定一个长度为n的数组a[1..n]，有一幅完全图，满足(u,v)的边权为a[u] xor a[v] 求边权和最小的生成树，你需要输出边权和还有方案数对1e9+7取模的值 1&amp;lt;=n&amp;lt;=10^5 0&amp;lt;=a[i]&amp;lt;2^30 分析 根据异或的性质，我们可以从高位到低位贪心。对于当前位，我们可以把点集分为当前位是1和是0两部分，显然最优的连边方案是两个集合自成...
[NOI2007]生成树计数
BeNoble_的博客
03-30 489
题意 一个nnn个点的树,点iii只能向[i−k,i−1][i−k,i−1][i-k,i-1]内的点连边,求有标号生成树的个数 n≤1015,k≤5n≤1015,k≤5n\le10^{15},k\le5 题解 题面误导向系列,显然不能像他说的那样做吧 因为发现kkk特别小,所以可以考虑用状压表示iii和前k−1k−1k-1个点的联通情况 考虑DP,f[i][j]DP,f[i][j]...
【51nod1806】wangyurzee的树（prufer序列）（容斥原理）
zxyoi_dreamer的博客（不定期诈尸）
07-14 181
传送门 题解： 关于度数的限制我们可以考虑放到prufer序列上。 一个点的度数就是它在prufer序列中出现次数+1。 所以问题就变成了一个长度为n−2n-2n−2的序列，每个位置的数可以是[1,n][1,n][1,n]中间的任意整数，限制某些数出现次数不能为某个值，求方案数。 很显然的容斥，考虑用总方案数减去所有不合法的方案，就是求满足0,1,2⋯m0,1,2\cdots m0,1,2⋯m个...
51NOD 1806 wangyurzee的树（容斥原理 + 组合数学）
代码改变世界
06-29 892
传送门 wangyurzee有n个各不相同的节点，编号从1到n。wangyurzee想在它们之间连n-1条边，从而使它们成为一棵树。 可是wangyurzee发现方案数太多了，于是他又给出了m个限制条件，其中第i个限制条件限制了编号为u[i]的节点的度数不能为d[i]。 一个节点的度数，就是指和该节点相关联的边的条数。 这样一来，方案数就减少了，问题也就变得容易了，现在请你告诉wangyur
51Nod-1806-wangyurzee的树
逐梦者
07-09 511
ACM模版描述题解尝试了这个题，百度到需要用到一种我没有接触过的序列，叫做 prufer 序列，十分强大的一个工具，网上查查能查到很多关于这个的讲解，这里有一个比较重要的是，prufer 序列的每一种序列对应一种生成树，序列中每个结点出现的次数比每个结点在该生成树中的度数少一，然后就是组合容斥的问题了。贴一下官方题解，和网上大神们讲的没有太大区别。这里注意要用到欧拉函数法求阶乘逆元，因为组合的公式里
51nod:51nod.com的代码
05-03
51点 51nod.com的代码
leetcodeoj和leetcode-My51Nod:个人题解（画风和leetcode不一样的OJ网站）
07-07
【标题】 "LeetCode OJ和LeetCode-My51Nod:个人题解与开源系统解析" 在编程世界中，LeetCode是一个广受欢迎的在线编程挑战平台，它提供了丰富的算法题目供程序员们练习和提升自己的技能。LeetCode OJ（Online Judge...
算法-树的直径（51Nod-2602）.rar
09-16
标题中的“算法-树的直径（51Nod-2602）”指的是一个算法问题，来源于在线编程挑战平台51Nod的第2602题。这道题目通常涉及图论中的树结构以及如何计算树的直径。树的直径是指在所有可能的路径中，最长的一条路径的...
BZOJ 1494: [NOI2007]生成树计数
frog的小池塘
07-25 2501
这题确实比较难， 所以我写解题报告也尽量从浅显处着手，希望能够帮到大家。 题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1494 题目大意：有n(n 算法： 从这道题的数据范围，不难看出要使用状态压缩。 首先，我们来确定前k个点的状态。
利用原始socket简单实现FTP的客户端和服务器端程序
热门推荐
ikunkun博客
07-06 2万+
1.设计目的 本设计旨在利用原始socket简单实现FTP(File Transfer Protocol,文件传输协议)的客户端和服务器端程序，能够实现get、put、pwd、dir、cd等基本交互命令。2.具体要求 用socket 编程接口编写两个程序，分别为客户程序(client.c)和服务器程序（server.c），该程序应能实现下述命令功能： get：取远方的一个文件 put：传给
计算机网络socket搭建实验（付源码）
淡风wisdon－大大的博客
12-26 2429
大二期间学习计算机网络时，做的实验：       本课程设计主要说明了在windows下如何用c语言编写设计ftp服务器和客户端的基本方法及其实现。包含了TCP/IP协议基本通信原理，Winsock编程基本原理以及ftp工作原理。       关键词：TCP/IP   Winsock   ftp   文件传输 文件传输源程序代码 （1）客户端client程序代码 ---- cl
02文件传输程序
MarkLee0419的博客
06-19 1020
1. FTP传输协议具体内容 1.1 实验目的 文件传送是各种计算机网络都实现的基本功能，文件传送协议是一种最基本的应用层协议按照客户/服务器的模式进行工作，提供交互式的访问，是INTERNET使用最广泛的协议之一。 本实验的目的是，学会利用已有网络环境设计并实现简单应用层协议，掌握TCP/IP 网络应用程序基本的设计方法和实现技巧。 1.2 实验内容和要求 1、实验内...
生成树的计数Matrix-Tree定理
07-19 1万+
摘要          在信息学竞赛中，有关生成树的最优化问题如最小生成树等是我们经常遇到的，而对生成树的计数及其相关问题则少有涉及。事实上，生成树的计数是十分有意义的，在许多方面都有着广泛的应用。本文从一道信息学竞赛中出现的例题谈起，首先介绍了一种指数级的动态规划算法，然后介绍了行列式的基本概念、性质，并在此基础上引入Matrix-Tree定理，同时通过与一道数学问题的对比，揭示了该定理所包
C语言按两个字节读写二进制文件,C语言 读写二进制文件(示例代码)
weixin_39742568的博客
05-18 576
查找了比较多的资源， 发现没有办法把text 文件转成binary文件仅作为记录，不过这个例子可以去除换行符。#include#include#defineN255intmain(){chara[N];FILE*fp1,*fp2;fp1=fopen("test_seq.fa","r");fp2=fopen("testSeq.dat","wb");/*inputtextfile...
生成树计数问题——矩阵树定理及其证明
不来也不去的一只失忆蝴蝶
02-07 1万+
生成树计数问题给一副n个节点的无向图G，求一个包含n-1条边的边集使得边集的边构成一颗树，问这样的边集的数量。矩阵树定理以下我们都不对重边与自环进行讨论。 先定义度数矩阵D，是一个n*n的矩阵。 Di,i=节点i的度数，对于i不等于j，Di,j=0。 再定义邻接矩阵A，也是一个n*n的矩阵。 i与j有边相连就有Ai,j=1否则Ai,j=0。 最后定义基尔霍夫矩阵C=D-A。 那么，Ci,
CF692D2：The Wall (medium)（组合数）
junior19的博客
08-23 1977
D2. The Wall (medium) time limit per test 2 seconds memory limit per test 256 megabytes input standard input output standard output Heidi the Cow is aghast: cracks in th
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值