【每日一题】6月12日题目精讲 数论

题号 NC13884
名称 Paint Box
来源 西南交通大学第十三届ACM决赛
戳我进入往期每日一题汇总贴~
往期每日一题二期题单

图片说明

感谢@p3ng推荐题目,100牛币奖励已发放~ 如果你在题库做题时遇到了喜欢的题目,欢迎推荐给邓老师~ 点击查看详情

题解

定义用不超过 种颜色涂完 个盒子的方案数为

意为:第一个盒子有 种选法,第二个之后的每一个盒子都需要和前一个不一样即为 种选法,共选了次。

于是,根据容斥原理,用正好 种颜色涂完 个盒子的方案数为:
题中给了 种颜色,从 种颜色选 k 种的方案数为

题目中要求的 n 个盒子,用 m 种颜色中选出 正好 k 种颜色,且相邻盒子颜色不同的方案数为:
于是,设
则有:

注意,在这里,并不代表数学上的含义,而是找到了一个利于计算的通式
先看特殊值
n=m=k=4,ans=24
n=m=k=5,ans=120
n=m=k 时,ans = k!
m=k>n时,ans=0
这些可以作为编程时的边界条件
考虑到n,m,k的数值,需要用到快速幂与逆元与组合数公式运算

注意到中的 的值,是不能直接用 O(1) 时间进行代入计算的。于是,可以考虑打表

但是,我们退回来想,利用式子 呢? 可以利用公式
即第一部分可以在O(k)的时间复杂度,利用阶乘和逆元的阶乘得到
第二部分可以利用阶乘、逆元和逆元的阶乘得到

代码如下:

#include<stdio.h>
#include<cstdio>
#include<iostream>
using namespace std;

#define ll long long

const ll mod=1e9+7;
const int maxn=1e6+10;

int T;
ll p[maxn], invp[maxn], inv[maxn];

//(a^b)%p
ll qp(ll a, ll b){
    a %= mod;
    ll ret = 1;
    while(b){
        if (b & 1)
            ret = ret * a % mod;
        a = (a * a) % mod;
        b >>= 1;
    }
    return ret;
}

void init(){
    p[0] = 1;
    for(int i = 1; i <= 1000000; i++)
        p[i] = p[i - 1] * i % mod;
    invp[1000000] = qp(p[1000000], mod - 2);
    for(int i = 1000000 - 1; i>=0; i--)
        invp[i] = invp[i + 1] * (i + 1) % mod;
    inv[1] = 1;
    for(int i = 2; i <= 1000000; i++)
        inv[i] = 1ll * (mod - mod / i) * inv[mod % i] % mod;
}

ll C(ll n, ll m){
    return p[n] * invp[m] % mod * invp[n - m] % mod;
}

int main()
{
    init();
    scanf("%d",&T);
    while(T--)
    {
        ll n,m,k;
        scanf("%lld%lld%lld",&n,&m,&k);
        ll Cmk = 1;
        for(int i=1; i<=k; i++)
            Cmk = Cmk * (m - k + i) % mod * inv[i] % mod;
        ll tmp = 0;
        for(int i=0; i<k; i++)
            if (i % 2 == 0)
                tmp = (tmp + C(k, k-i) * (k - i) % mod * qp(k - i - 1, n - 1) % mod + mod) % mod;
            else
                tmp = (tmp - C(k, k-i) * (k - i) % mod * qp(k - i - 1, n - 1) % mod + mod) % mod;
        ll ans = Cmk * tmp % mod;
        printf("%lld\n", ans);
    }
    return 0;
}

欢迎各位大佬来做题写题解,也欢迎大家踊跃在当日讨论贴中提问!

活动奖励:

在牛客博客中写出题解,并回复地址
审核通过可获得(依据题目难度和题解的内容而定)

本道题目6月19日中午12:00之前写的题解有获得牛币资格~

.牛币兑换中心

牛客博客开通方式

  1. 如何开通牛客博客:https://www.nowcoder.com/discuss/202952
  2. 如何使用博客搬家功能:进入博客--->设置--->底部博客搬家
  3. 如果你对牛客博客有任何意见或建议:牛客博客意见反馈专贴
全部评论
大家也这么怕数学题么? 那以后要循序渐进多来点
1 回复 分享
发布于 2020-06-11 18:49
https://blog.nowcoder.net/n/d938236cec1048f58229c85ee2a8c014
1 回复 分享
发布于 2022-11-02 14:09 山东
https://blog.nowcoder.net/n/03dde93fab5c4cbd8c5a0cc024f7d72f
点赞 回复 分享
发布于 2020-06-16 15:26
https://blog.nowcoder.net/n/5def76ea197748f284c250b9ebe6ef2d
点赞 回复 分享
发布于 2020-06-12 17:36
点赞 回复 分享
发布于 2020-06-11 18:29
点赞 回复 分享
发布于 2020-06-11 14:31
https://blog.nowcoder.net/n/56522415c491401e88539f49911bc388
点赞 回复 分享
发布于 2020-06-11 13:44
点赞 回复 分享
发布于 2020-06-11 12:47

相关推荐

点赞 评论 收藏
分享
07-11 15:12
门头沟学院 Java
别人在上班,我就在工位上看看视频啥的,这正常吗?
程序员小白条:实习就是摸鱼,只是公司指标,把你进来了,可能那时候客户很多,但等你进来的时候,已经是淡季了,根本没多少需求,或者说根本不适合实习生去完成,因此你就每天干坐着就行,可能1,2个月都没需求
实习生的蛐蛐区
点赞 评论 收藏
分享
避坑恶心到我了大家好,今天我想跟大家聊聊我在成都千子成智能科技有限公司(以下简称千子成)的求职经历,希望能给大家一些参考。千子成的母公司是“同创主悦”,主要经营各种产品,比如菜刀、POS机、电话卡等等。听起来是不是有点像地推销售公司?没错,就是那种类型的公司。我当时刚毕业,急需一份临时工作,所以在BOSS上看到了千子成的招聘信息。他们承诺无责底薪5000元,还包住宿,这吸引了我。面试的时候,HR也说了同样的话,感觉挺靠谱的。于是,我满怀期待地等待结果。结果出来后,我通过了面试,第二天就收到了试岗通知。试岗的内容就是地推销售,公司划定一个区域,然后你就得见人就问,问店铺、问路人,一直问到他们有意向为止。如果他们有兴趣,你就得摇同事帮忙推动,促进成交。说说一天的工作安排吧。工作时间是从早上8:30到晚上18:30。早上7点有人叫你起床,收拾后去公司,然后唱歌跳舞(销售公司都这样),7:55早课(类似宣誓),8:05同事间联系销售话术,8:15分享销售技巧,8:30经理训话。9:20左右从公司下市场,公交、地铁、自行车自费。到了市场大概10点左右,开始地推工作。中午吃饭时间大约是12:00,公司附近的路边盖饭面馆店自费AA,吃饭时间大约40分钟左右。吃完饭后继续地推工作,没有所谓的固定中午午休时间。下午6点下班后返回公司,不能直接下班,需要与同事交流话术,经理讲话洗脑。正常情况下9点下班。整个上班的一天中,早上到公司就是站着的,到晚上下班前都是站着。每天步数2万步以上。公司员工没有自己的工位,百来号人挤在一个20平方米的空间里听经理洗脑。白天就在市场上奔波,公司的投入成本几乎只有租金和工资,没有中央空调。早上2小时,晚上加班2小时,纯蒸桑拿。没有任何福利,节假日也没有3倍工资之类的。偶尔会有冲的酸梅汤和西瓜什么的。公司的晋升路径也很有意思:新人—组长—领队—主管—副经理—经理。要求是业绩和团队人数,类似传销模式,把人留下来。新人不能加微信、不能吐槽公司、不能有负面情绪、不能谈恋爱、不能说累。在公司没有任何坐的地方,不能依墙而坐。早上吃早饭在公司外面的安全通道,未到上班时间还会让你吃快些不能磨蹭。总之就是想榨干你。复试的时候,带你的师傅会给你营造一个钱多事少离家近的工作氛围,吹嘘工资有多高、还能吹自己毕业于好大学。然后让你早点来公司、无偿加班、抓住你可能不会走的心思进一步压榨你。总之,大家在找工作的时候一定要擦亮眼睛,避免踩坑!———来自网友
qq乃乃好喝到咩噗茶:不要做没有专业门槛的工作
点赞 评论 收藏
分享
星辰再现:裁员给校招生腾地方
点赞 评论 收藏
分享
评论
点赞
2
分享

创作者周榜

更多
牛客网
牛客网在线编程
牛客网题解
牛客企业服务