【组合数学】LGP3330（ZJOI2011）+UOJ#209

最新推荐文章于 2021-12-01 11:57:09 发布
最新推荐文章于 2021-12-01 11:57:09 发布
阅读量271 收藏
点赞数
分类专栏： 数论-组合数学
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Dream_Lolita/article/details/78782773
版权

LGP3330 看电影(MOVIE)

原题地址

【题目大意】
一共有n个人，k个格子（格子标号为1~k）
循环n轮，每轮随机一个k以内的整数（设为m），若m~k这些格子里有空的，就把一个人放到格子里，否则这个人将站着。

问全部人都坐着的概率是多少。

【题目分析】
一眼组合数学，实在不会可以考虑打表。
组合数学中经典的古典概型，概率等于合法的方案数除以总方案数。

【解题思路】
我们易得总方案数= KN 。
　
合法方案数的计算：考虑在最后面加一个座位，然后所有的座位连成一个环。
这样我们保证不论怎么放，最后都会有一个空的位置，那么我们从这个空的位置断开（将它丢到），必然是一个合法的情况。

现在那么总方案数等于 (K+1)N ，每种环都算了 (K+1) 次，所以除以 (K+1) ，最后抽调一个空白的座位来构造一个合法的方案，所以要乘 (N−K+1 )。

综上所述，答案就是：

(K+1)N−1×(N−K+1)KN

要用高精度 。

【代码】

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

const int MAXN=205;
const int MAXL=505;
int prinum,T,n,k,len;
int a[MAXL],pri[MAXN],ans[MAXN],ans2[MAXN];
int p[MAXN][MAXN];
bool ispri[MAXN];

inline void get_prime()
{
    memset(ispri,true,sizeof(ispri));
    ispri[1]=false;
    for(int i=2;i<MAXN;++i)
    {
        if(ispri[i])
        {
            ++prinum;
            pri[prinum]=i;
        }
        for(int j=1;j<=prinum && pri[j]*i<MAXN;++j)
        {
            ispri[pri[j]*i]=false;
            if(i%pri[j]==0)
                break;
        }
    }
}

inline void get_num()
{
    int x;
    for(int i=1;i<MAXN;++i)
    {
        x=i;
        for(int j=1;j<=prinum;++j)
            while(x%pri[j]==0)
            {
                p[i][pri[j]]++;
                x/=pri[j];
            }
    }
}

inline void _reset()
{
    memset(ans,0,sizeof(ans));
    memset(ans2,0,sizeof(ans2));
}

inline void mul(int x)
{
    if(len==0)
    {
        len=1;a[len]=x;
        for(len++;a[len];)
        {
            a[len+1]+=a[len]/10;
            a[len++]%=10;
        }
        --len;
    }
    else
    {
        for(int i=1;i<=len;++i)
            a[i]*=x;
        for(int i=1;i<=len;++i)
        {
            a[i+1]+=a[i]/10;
            a[i]%=10;   
        }
        for(len++;a[len];)
        {
            a[len+1]+=a[len]/10;
            a[len++]%=10;   
        }
        --len;
    }
}

inline void get_ans()
{   
/*  for(int i=1;i<10;++i)
        printf("%d ",ans[i]);
    printf("\n");
    for(int i=1;i<10;++i)
        printf("%d ",ans2[i]);
    printf("\n");*/

    len=0;
    memset(a,0,sizeof(a));
    for(int i=1;i<MAXN;++i)
        for(int j=1;j<=ans[i];++j)
            mul(i);
    if(!len)
        printf("1");
    for(int i=len;i;--i)
        printf("%d",a[i]);
    printf(" ");    

    len=0;  
    memset(a,0,sizeof(a));
    for(int i=1;i<MAXN;++i)
        for(int j=1;j<=ans2[i];++j)
            mul(i);
    if(!len)
        printf("1");
    for(int i=len;i;--i)
        printf("%d",a[i]);
    printf("\n");
}

inline void solve()
{   
    scanf("%d%d",&n,&k);
    if(n>k)
    {
        printf("0 1\n");
        return;
    }

    for(int i=1;i<=k+1;++i)
        ans[i]=p[k+1][i]*(n-1);

    if(k-n>0)
        for(int i=1;i<=k-n+1;++i)
            ans[i]+=p[k-n+1][i];
    if(k>1)
        for(int i=1;i<=k;++i)
        {
            ans2[i]=p[k][i]*n-ans[i];
            ans[i]-=p[k][i]*n;
        }

    get_ans();
}

int main()
{
    freopen("LGP3330.in","r",stdin);
    freopen("LGP3330.out","w",stdout);

    get_prime();
    get_num();

    scanf("%d",&T);
    while(T--)
    {
        _reset();
        solve();
    }

    return 0;
}

UOJ#209. 【UER #6】票数统计

原题地址

【题目大意】
一共有n个数字 ，每个数字要么是0，要么是1。
有m个限制 (x，y) ，表示前x个数字里有y个1，或后y个数字里有x个1
求满足所有限制的序列个数。

n≤5000，m≤1000，答案对998244353取模，一个点至多5组数据

【题目分析】
都是骗人的！

【解题思路】
当x!=y的时候，我们已经确定了这个限制是对前缀还是后缀的。
当x=y的时候，我们只需要保留最大的那个x就行了。

然后，一波组合数学。
答案是x=y是前缀限制的个数+后缀限制的个数-前后缀都限制的个数。

那么，现在就是有一堆前后缀的限制，让你求方案数。

如果只有前缀或只有后缀就是一个划分成m个区间的组合数学。

如果两个都有就把前缀转成后缀或者把后缀转成前缀。
枚举1的总和就行了。

【代码】

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

typedef long long LL;
const int MAXN=5010;
const int mod=998244353;

int n,m,ans,tot,lz,T;
int C[MAXN][MAXN];

struct Tnode
{
    int x,y,bz,ty;
    Tnode(){}
    Tnode(int xx,int yy,int bzz,int tyy)
    {
        x=xx;y=yy;
        bz=bzz;ty=tyy;
    }
};
Tnode a[MAXN],p[MAXN];

inline void get_C()
{
    C[0][0]=1;
    for(int i=1;i<MAXN-5;++i)
    {
        C[i][0]=1;
        for(int j=1;j<=i;++j)
            C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;
    }
}

inline bool cmp(Tnode A,Tnode B)
{
    return A.x<B.x;
}

inline void init()
{
    scanf("%d%d",&n,&m);

    tot=ans=lz=0;
    for(int i=1;i<=m;++i)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        if(x>y)
            a[++tot]=(Tnode){x,y,0,0};
        else
        if(x<y)
            a[++tot]=(Tnode){n-y,x,1,0};
        else
            lz=max(lz,x);
    }
    a[++tot]=(Tnode){lz,lz,0,1};
    a[++tot]=(Tnode){n-lz,lz,1,2};

    sort(a+1,a+tot+1,cmp);

//  for(int i=1;i<=tot;++i)
//      printf("%d %d %d %d\n",a[i].x,a[i].y,a[i].bz,a[i].ty);
}

inline void solve()
{
    for(int typ=1;typ<=3;++typ)
    {
        for(int i=1;i<=n;++i)
        {
            int cnt=0;bool flag=true;
            for(int j=1;j<=tot;++j)
            {
                if(a[j].ty==typ)
                    continue;
                p[++cnt].x=a[j].x;
                if(a[j].bz)
                    p[cnt].y=i-a[j].y;
                else
                    p[cnt].y=a[j].y;

                if(p[cnt].y<0 || p[cnt].y>i)
                {
                    flag=false;
                    break;
                }

                if(p[cnt].x<p[cnt].y)
                {
                    flag=false;
                    break;
                }

                if(p[cnt].x-p[cnt-1].x<p[cnt].y-p[cnt-1].y)
                {
                    flag=false;
                    break;
                }

                if(p[cnt].y-p[cnt-1].y<0)
                {
                    flag=false;
                    break;
                }
            }

            if(!flag)
                continue;

            int sum=1;
            for(int j=1;j<=cnt;++j)
                sum=1ll*sum*C[p[j].x-p[j-1].x][p[j].y-p[j-1].y]%mod;
            sum=1ll*sum*C[n-p[cnt].x][i-p[cnt].y]%mod;
//              printf("sum:%d\n",sum);

            if(typ<=2)
                ans=(ans+sum)%mod;
            else
                ans=(ans-sum+mod)%mod;
        }
    }
    printf("%d\n",ans);
}

int main()
{
//  freopen("UOJ209.in","r",stdin);
//  freopen("UOJ209.out","w",stdout);

    get_C();

    scanf("%d",&T);
    while(T--)
    {
        init();
        solve();
    }

    return 0;
}
Dream_Lolita
关注
专栏目录
UOJ#188. 【UR #13】Sanrd（min25筛）
m0_52639539的博客
05-02 355
UOJ#188. 【UR #13】Sanrd（min25筛） 题目链接：https://uoj.ac/problem/188 题目大意：定义数论函数 f(x)f(x)f(x) 的值为 xxx 的第二大质因子，其中质因子可以重复贡献，例如 f(36)=2f(36)=2f(36)=2 ，
【UOJ #209】【UER #6】票数统计
LYD729
07-06 644
Description 妹滋滋是一个善于编程的女孩子。 但是某一天，她一不小心把 UOJ 后台的票数统计程序写错了。 本来嘛在这种根本没有什么用的功能上出了 bug 也没有什么大关系，但是又有某一天，UOJ 突然就开始搞全民公投了。 这可怎么办呢？如果这个消息让别人知道的话自己肯定会被查表，更不要说让所有用户重新来投一次票了。 作为一个要强的女孩子，妹滋滋决定自力更生。 通过一些奥妙重
UOJ#209. 【UER #6】票数统计 组合+容斥
EM LGH
06-16 146
自己独立想出来的，开心. 首先，对于 $x$ 不等于 $y$ 的情况，显然只限制前缀/后缀. 然后如果没有 x 与 y 相等的情况的话我们完全可以枚举总的 1 的个数，然后限制后缀的就可以转化成限制前缀的. 如果引入 x 与 y 相等的情况，也就是说要求前缀或后缀都填满也按照上述方式处理即可. 但是要简单容斥一下，即 ans(前缀)+ans(后缀)-a...
[uoj#209][UER#6A]票数统计
时光真疯狂, 我一路执迷于匆忙.
07-04 890
Description给出n个数，每个数是0或1. 再给出m个限制，每个限制(x,y)表示“前x个数中有y个1”或“后y个数中有x个1” 求这样的序列的个数。 n<=5000,m<=1000Solution再一次被UER给虐了。 其实这道题劼鏼爷已经讲的很清楚了。（扑通扑通跪下来） 当x!=y的时候，很显然已经确定这个限制是限制前缀还是后缀的。 当x=y的时候，我们只需要保留最大的那个x
[我可能是个SB啊] UOJ #290 【ZJOI2017】仙人掌
雯舞
03-24 853
这好像不难啊 考场上智商掉线了吧 可能是脑袋被吃掉了吧先说正经的 这个题 要是不是棵仙人掌 那么就无解 要是是棵树 直接树形DP 那仙人掌呢 把所有环上的边拆掉 对剩下的森林中的每棵树DP 然后乘法原理就好了 然后 做完了…还是不知道考场上在干什么我可能忍不住分享下自己的悲惨经历啊 考场上想出了树形DP 然后就没有然后了？ 尼玛仙人掌这么好想 就被我弃了？！ 行吧 考挂不如狗 下午切这
《从Unknown谈一类支持末尾插入删除的区间信息维护方法》学习笔记+UOJ #191代码（及HACK5原理，雾）
litble的成(tui)长(fei)史
05-21 461
笔记 做带末尾插入删除的区间信息维护）的数据结构题的方法： 分块 思路：每次插删操作暴力重构最后一块。 支持插删操作，支持区间查询，支持在线 二进制分组 思路：若每次添加一个元素进数据结构里的复杂度比较高，则每次将这个元素单独放在最后一组，若最后一组与上一组的大小相同，就将这两组合并为同一组，不难发现最后得到的每个组大小都是2的次幂，并且互不相同，复杂度会是log⁡\loglog级的。 不支持删除...
[BZOJ4584/UOJ#204][APIO2016]Boat 赛艇（DP+组合数学+前缀和）
女装的你，如此好看！
04-14 270
容易得出这是一个关于上升子序列的DP。 一个直接的想法： f[i][j]f[i][j]f[i][j] 表示以 iii 为末尾，第 iii 个数的取值为 jjj 的方案数，但 aiaia_i 和 bibib_i 的规模是 10910910^9 ，显然过不去。 考虑离散化出一些区间 [li,ri][li,ri][l_i,r_i] ，使得每个 [li,ri][li,ri][l_i,r_i] 都满足：...
【uoj#209】[UER #6]票数统计 组合数+乱搞
weixin_30563319的博客
03-31 97
题目描述 一个长度为 $n$ 的序列，每个位置为 $0$ 或 $1$ 两种。现在给出 $m$ 个限制条件，第 $i$ 个限制条件给出 $x_i$ 、$y_i$ ，要求至少满足以下两个条件之一： 序列的前 $x_i$ 个位置中，恰好有 $y_i$ 个 $1$ ； 序列的后 $y_i$ 个位置中，恰好有 $x_i$ 个 $1$ ； 求有多少个序列满足所有限制条件。答案可能很大...
UOJ Round #20 T1 A. 【UR #20】跳蚤电话（组合数+树形DP）
ZSJZ_liuzian的博客
04-05 327
UOJ Round #20 T1 A. 【UR #20】跳蚤电话 题目大意 给出一棵树，求建出该树的不同操作方案数。建树方式如下：初始SSS集合只有111，操作111为取已连的边x,yx,yx,y和不在SSS的点zzz，删去边(x,y)(x,y)(x,y)，加入边(x,z),(y,z)(x,z),(y,z)(x,z),(y,z)，再把zzz放入SSS；操作222为取SSS内的点xxx和SSS外的点yyy，加入边(x,y)(x,y)(x,y)，再把yyy放入SSS。 2≤n≤1052\le n \le 10
UOJ#351-新年的叶子【树的直径,数学期望】
QuantAsk
12-01 1491
正题 题目链接:https://uoj.ac/problem/351 题目大意 给出nnn个点的一棵树，开始所有点都是白色，每次随机点黑一个叶子（可以重复点），求期望多少次能使得白色点构成的图直径发生变化。 答案对998244353998244353998244353取模 1≤n≤5×1051\leq n\leq 5\times 10^51≤n≤5×105 解题思路 考虑什么时候会直径会产生变化。 假设直径的长度LLL为偶数，那么所有的直径都有一个共同的中心点，设为xxx。此时我们需要在xxx的两棵子树
uoj #176. 新年的繁荣 字典树+Boruvka算法
beginend
03-06 1320
题意 有n个点，第i个点的权值为a[i]，在第i个点和第j个点之间连边的代价为a[i] and a[j]。问这个图的最大生成树。 n&lt;=100000,a[i]&lt;218n&lt;=100000,a[i]&lt;218nO(m2m)O(m2m)O(m2^m)。 我们可以在字典树每个节点维护其子树中连通块编号的最大值和最小值，这样就能判断子树中是否存在不同的连通块。 那么总的时间复杂...
【UOJ#37】 [清华集训2014] 主旋律
Yu Gi Oh!
05-16 287
题目链接 题目描述 给定一张强联通图,求有多少种边的存在情况满足图依然强联通。 n≤15n\leq15n≤15 Sol 首先正难则反,考虑用总数减去不强联通的。 考虑一张不强联通的图，缩点后一定是一个 DAG，好像可以对 DAG 进行计数。 诈一看这个做不了，因为缩点后计数是不可能在dp过程中实现的。 但我们按照 DAG 计数的思路的话其实并不需要真的知道 DAG 缩点后的形态。 我们类似 DAG...
UOJ #164. 【清华集训2015】V
Freopen的博客
01-21 313
题目 （居然）比想象中好写多了 操作： 区间赋值，区间加，区间减，区间和000取maxmaxmax，单点询问，单点求历史最大值。 此题所有的修改都可以用一个标记表达： (a,b)(a,b)(a,b)表示将xxx变为max(a+x,b)max(a+x,b)max(a+x,b) 那么标记的合并可以发现(a,b)∗(c,d)=(a+c,max(b+c,d))(a,b)*(c,d)=(a+c,max(b+...
UOJ#136. 【UER #3】开学前的作文
reflash的博客
12-09 558
到底是谁一直在错？
【组合数学+转化问题】BZOJ4005[JLOI2015]骗我呢
MEET YOU FIRST TIME
02-01 1340
一道不错的组合数学 转化问题的思路很优秀
【概率/组合数学】CF1060F Shrinking Tree
MEET YOU FIRST TIME
10-06 1134
【题目】 原题地址 给定一棵nnn个点的树，进行n−1n-1n−1轮操作，每轮操作随机选择一条边(u,v)(u,v)(u,v)，将u,vu,vu,v两点合成一个点，即删去这两个点后新建一个点xxx，将原来与uuu或vvv连边的点连向xxx，接着xxx的编号随机为u,vu,vu,v中的一个。 对于每个编号111~nnn，求它最终留下的概率。n≤50n\leq 50n≤50 【题目分析】 我的概率好弱...
【多项式转化+分治FFT+大整数处理】CC_LUCASTH Lucas Theorem
MEET YOU FIRST TIME
03-12 487
【题目】 Codechef 集合[n][n][n]表示由1∼n1\sim n1∼n所有整数构成的集合，定义： f(n,k)=∑S⊆[n]∣S∣=k∏x∈Sx f(n,k)=\sum_{S\subseteq [n]\\|S|=k} \prod_{x\in S}x f(n,k)=S⊆[n]∣S∣=k∑​x∈S∏​x 特别地，f(n,0)=1f(n,0)=1f(n,0)=1，求对于所有k∈[0,n]k\...
【数论-筛法/五边形数/组合数学】BZOJ4772 显而易见的数论
MEET YOU FIRST TIME
09-02 479
【题目】 原题地址 给你一个整数SSS，让你求所有整数划分的方案数的价值和，价值是个函数。 ∑划分方案p∑i=1m∑j=i+1mg(aF(pi,pj)mod&nbsp;k)∑划分方案p∑i=1m∑j=i+1mg(aF(pi,pj)mod&nbsp;k)\sum_{划分方案p}\sum_{i=1}^m\sum_{j=i+1}^mg(a_{F{(p_i,p_j)}mod\ k}) 其中F(pi...
【组合计数+NTT优化卷积】BZOJ5306 [HAOI2018] 染色
MEET YOU FIRST TIME
03-13 376
【题目】 lydsy 一个长度为nnn的序列，每个位置可以被染成mmm种颜色中的一种。若一种方案中出现次数恰好为SSS的颜色数有KKK种，则会有WKW_KWK​的愉悦值。问所有方案的愉悦值总和对100453580910045358091004535809取模的结果。 n≤107,m≤105,S≤150n\leq 10^7,m\leq 10^5,S\leq 150n≤107,m≤105,S≤150 ...
UOJ#435题目“Simple Tree”
最新发布
05-30
UOJ#435题目“Simple Tree”要求实现一种数据结构，支持以下两种操作： 1. 修改某个节点的权值； 2. 查询某个节点到根节点路径上权值的最大值。 因为需要查询某个节点到根节点的路径，所以可以考虑使用树链剖分...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值