描述

Ssoier在紧张的学习中,杜老师每天给他们传授精妙的知识。

杜老师为了活跃气氛,设计了一个点名器,这个点名器包含一个长度为M的数组(下标1开始),每个元素是一个oier的名字,每次点名的时候,点名器会等概论随机生成一个1到M的整数,对应的人就要回答问题。当然杜老师有喜欢的人,他会故意让一些名字重复出现以增加某些人被点中的概率。

LDX感觉不是很好,他不希望被点中。他找到杜老师,不过杜老师表示他确实喜欢ldx,不过杜老师还是想给ldx一点好处,他被点名器的数组给ldx看,然后让ldx对数组进行T此修改,每次修改LDX选择的某个位置的名字,然后把他改成其他oier的名字(一定是换成另外人的名字,不允许不修改)。LDX希望让自己在点名器的出现次数不超过K。

请你帮帮LDX计算有多少种满足条件的修改方案。两个方案不同当且仅当存在i,两方案中的第i次修改操作不同。答案对1e9+7取模。

输入

为了简化问题,为每个OIer分配一个1 到N 的整数。zgs是1号。

第一行三个整数N,M,T,K,含义如上。

接下来一行M 个整数表示数组元素对应的选手编号

输出

输出一行一个整数表示答案。

样例输入

2 3 1 0

2 2 2

样例输出

0

提示

对于20% 的数据,N,M,T <= 10。

对于50% 的数据,N,M,T<= 50。

对于70% 的数据,N,M,T<=200

对于100% 的数据,1<=M,T<=2000; 2<=N<=10^9; 0<=K<=M。


这是一个简单的dp。

直接f[i][j]" role="presentation" style="position: relative;">f[i][j]f[i][j]表示前i轮剩下j个数值为1的方案数量。

显然f[i][j]只与f[i-1][j],f[i-1][j-1],f[i-1][j+1]有关。

随便转移就能过了。

代码:

#include<bits/stdc++.h>
#define ll long long
#define mod 1000000007
#define N 2005
using namespace std;
inline ll read(){
    ll ans=0;
    char ch=getchar();
    while(!isdigit(ch))ch=getchar();
    while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
    return ans;
}
ll n,m,t,k,a[N],f[N][N],ans=0,cnt=0;
int main(){
    n=read(),m=read(),t=read(),k=read();
    for(ll i=1;i<=m;++i){a[i]=read();if(a[i]==1)++cnt;}
    f[0][cnt]=1;
    for(ll i=1;i<=t;++i){
        for(ll j=m;~j;--j){
            f[i][j]=(f[i-1][j+1]*(j+1)%mod*(n-1)%mod+f[i-1][j]*(m-j)%mod*(n-2)%mod)%mod;
            if(j)(f[i][j]+=f[i-1][j-1]*(m-j+1))%=mod;
        }
    }
    for(int i=k;~i;--i)(ans+=f[t][i])%=mod;
    cout<<ans;
    return 0;
}

2018.09.15点名器(简单dp)的更多相关文章

  1. Lean Data Innovation Sharing Salon(2018.09.15)

    时间:2018.09.15地点:北京国华投资大厦

  2. 4.15 每周作业 —— 简单DP

    免费馅饼 Time Limit : 2000/1000ms (Java/Other)   Memory Limit : 65536/32768K (Java/Other) Total Submissi ...

  3. 2018.09.01 poj3071Football(概率dp+二进制找规律)

    传送门 概率dp简单题. 设f[i][j]表示前i轮j获胜的概率. 如果j,k能够刚好在第i轮相遇,找规律可以发现j,k满足: (j−1)>>(i−1)" role=" ...

  4. 2018.09.01 独立集(树形dp)

    描述 给定一颗树(边权为1),选取一个节点子集,使得该集合中任意两个节点之间的距离都大于K.求这个集合节点最多是多少 输入 第一行是两个整数N,K 接下来是N-1行,每行2个整数x,y,表示x与y有一 ...

  5. 2018.09.15 hdu3018Ant Trip(欧拉路)

    传送门 显然答案等于各个连通分量的笔画数之和. 因此我们dfs每个连通分量计算对答案的贡献. 对于一个连通分量,如果本来就有欧拉回路那么只需要一笔. 否则需要寄点数/2那么多笔才能画完. 知道这个结论 ...

  6. 2018.09.15模拟总结(T1,T3)

    过了一周,终于迎来了第二次模拟(这不是期待的语气),看第一周毒瘤程度,我就觉得接下来的模拟只能更毒瘤. 花了10多分钟读完了三道题,觉得暴力还是挺好写的,然后在每一道题都思索那么几分钟后,觉得还是写暴 ...

  7. hdu 2018 母牛的故事(简单dp)

    母牛的故事 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Submi ...

  8. 2018.09.15 poj1734Sightseeing trip(floyd求最小环)

    跟hdu1599差不多.. 只是需要输出方案. 这个可以递归求解. 代码: #include<iostream> #include<cstdio> #include<cs ...

  9. 2018.09.15 hdu1599find the mincost route(floyd求最小环)

    传送门 floyd求最小环的板子题目. 就是枚举两个相邻的点求最小环就行了. 代码: #include<bits/stdc++.h> #define inf 0x3f3f3f3f3f3f ...

随机推荐

  1. leetcode476

    public class Solution { public int FindComplement(int num) { //计算数字二进制的反码 var list = new List<int ...

  2. 免费json API

    免费json API http://www.bejson.com/knownjson/webInterface/

  3. TTreeView.OnCustomDrawItem

    TTreeNode *node; node = , "AAAA"); TreeView1->Items->AddChild(node, "aaa1" ...

  4. VB6 如何创建一个标准控制台程序

    打开 VB6 并新建一个标准EXE程序,把窗口删掉,然后再加入一个模块. 在模块中加入AllocConsole.FreeConsole.SetConsoleTitle.Sleep的API声明: Pub ...

  5. Tocmat 启动错误 Port 8005 required by tomcat v7.0 server at localhost is already in use

    1: netstat -ano!findstr 8005 2: taskkill /pid 6476 /f

  6. os.popen与os.system区别

    os.system 调用系统命令,完成后退出,返回结果是命令执行状态,一般是0 os.popen 可以实现一个“管道”,从这个命令获取的值可以在python 中继续被使用 #该方法不但执行命令还返回执 ...

  7. 分享 - 27 个机器学习、数学、Python 速查表

      转载自:伯乐在线 - iPytLab,原文链接,侵删 机器学习涉及到的方面非常多.当我开始准备复习这些内容的时候,我找到了许多不同的”速查表”, 这些速查表针对某一主题都罗列出了所有我需要知道的知 ...

  8. object-c语法

    Objective-C:C的超集 Objective-Objective-C是C语言的严格超集--任何C语言程序不经修改就可以直接通过Objective-C编译器,在Objective-C中使用C语言 ...

  9. delphi 图片加水印源代码

    unit UWaterMark; interface uses {$IFNDEF DELPHIXE2ANDUP} windows,SysUtils,classes,graphics,Gdiplus; ...

  10. spring 中的断言的作用

    org.springframework.util.AssertAssert翻译为中文为"断言".用过JUNIT的应该都知道这个概念了.就是断定某一个实际的值就为自己预期想得到的,如 ...