7_11模拟赛

最新推荐文章于 2023-03-12 16:52:22 发布
最新推荐文章于 2023-03-12 16:52:22 发布
阅读量303 收藏 1
点赞数
分类专栏： gfoj 比赛 dp 数论-数学 文章标签： 动态规划
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/mc_dl/article/details/75007201
版权
gfoj 同时被 3 个专栏收录
26 篇文章 0 订阅
订阅专栏
dp
10 篇文章 0 订阅
订阅专栏
数论-数学
3 篇文章 0 订阅
订阅专栏

gfoj

Problem 720: Chess

Time Limit: 1000 ms Memory Limit: 255998 KB

Problem Description
在一个n*m的棋盘上，要求每一行每一列都至少有一个棋子的总方案数是多少？由于答案有可能很大，所以输出答案对1000000007取模。

Input
两个整数n和m。

   30%                1<=n,m<=5

   100%              1<=n,m<=50

Output
一个整数表示总方案数模1000000007的值。

Sample Input
2 3

Sample Output
25

本题各种做法

dp(n^3)

dp[ i ] [ j ] 表示前i行中有j列放有棋子的方案数

#include<cstdio>
#include<iostream>
#define MOD 1000000007
#define N 60
using namespace std;
typedef long long LL;
LL a[N],dp[N][N],p,c[N][N];
int n,m;
int main(){
    scanf("%d%d",&n,&m);
    c[0][0]=1;
    for(int i=1;i<=m;i++)
        for(int j=0;j<=i;j++)
            if(j==0)c[i][j]=1;
            else c[i][j]=(c[i-1][j-1]+c[i-1][j])%MOD;
    a[0]=1;
    for(int i=1;i<=50;i++)a[i]=a[i-1]*2%MOD;
    dp[0][0]=1;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            for(int k=0;k<j;k++){
                p=(dp[i-1][k]*a[k]%MOD)*c[m-k][j-k]%MOD;
                dp[i][j]=(dp[i][j]+p)%MOD;
            }
            p=dp[i-1][j]*(a[j]-1)%MOD;
            dp[i][j]=(dp[i][j]+p)%MOD;
        }
    }
    cout<<dp[n][m];
}

lzk(n^4)

#include<cstdio>
#define ll long long
using namespace std;

const int N=50+5;
const ll mod=1e9+7;

ll c[N][N];
ll f[N][N];

ll qpow(int x) {
    if(x==0) return 1ll;
    ll ans=qpow(x/2);
    ans=ans*ans%mod;
    if(x&1) ans=ans*2%mod;
    return ans%mod;
}

int main() {
    c[0][0]=1;
    for(int i=1; i<=50; i++) {
        c[i][0]=1;
        for(int j=1; j<=i; j++) c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;
    }

    int n, m;
    scanf("%d%d", &n, &m);
    f[1][1]=1;
    for(int x=1; x<=n; x++)
        for(int y=1; y<=m; y++) {
            if(x==1&&y==1) continue;
            f[x][y]=qpow(x*y);
            for(int i=0; i<x; i++)
                for(int j=0; j<y; j++) {
                    if(i==0&&j==0) continue;
                    f[x][y]=(f[x][y]+mod-c[x][i]*c[y][j]%mod*f[x-i][y-j]%mod)%mod;
                }
            f[x][y]=(f[x][y]-1)%mod;
        }
    printf("%lld\n", f[n][m]);
    /*
    printf("\n");
    for(int i=1; i<=n; i++) {
        for(int j=1; j<=m; j++) printf("%lld ", f[i][j]);
        printf("\n");
    }
    */
    return 0;
}

hz 容斥原理

∑i=0n(−1)i×Cin×(2n−i−1)m

(2n−1)m−C(n,1)∗(2n−1−1)m+C(n,2)∗(2n−2−1)m−....+...

最后是加或减要看n的奇偶 

#define mod 1000000007
using namespace std;

typedef long long ll;
int m,n,i;
ll ans;
ll pow(ll t,int k) {
    ll res = 1;
    for (int i=0;(k>>i);i++) {
        if ((k>>i)&1) res = (res * t) % mod;
        t = (t * t) % mod;
    } return res;
}

ll C(int a,int b) {
    ll res = 1;
    if (b*2 > a) b = a - b;
    for (int i=0;i<b;i++) res = ((res * (a-i)) % mod) * pow(i+1,mod-2) % mod;
    return res;
}

int main() {
    cin >> n >> m;
    for (i=n;i;i--) {
        if ((n-i) & 1) ans = (ans + mod - (C(n,n-i) * pow(pow(2,i)-1,m) % mod)) % mod;
        else ans = (ans + C(n,n-i) * pow(pow(2,i)-1,m)) % mod;
    }
    cout << ans;
}

ssn&lsb 递归

#include<iostream>
#include<cstring>
#define ll long long
using namespace std;
const int maxn=60, MOD=1000000007;
ll f[maxn], ys[maxn], n, m;
ll mi(ll a, ll k) {
    if(k == 0) return 1;
    if(k == 1) return a%MOD;
    ll l = mi(a, k/2)%MOD;
    if(k & 1) return ((l*l%MOD)*a)%MOD;
    else return (l*l)%MOD;
}
ll C(ll a, ll k) {
    memset(ys, 0, sizeof(ys));
    ll ret = 1;
    for(int i = a; i >= a-k+1; i--) {
        int b = i;
        for(int j = 2; j <= b; j++) while(b%j == 0)
            {ys[j]++; b /= j; }
    }
    for(int i = k; i >= 1; i--) {
        int b = i;
        for(int j = 2; j <= b; j++) while(b%j == 0) 
            {ys[j]--; b /= j; }
    }
    for(int i = 1; i <= 50; i++) while(ys[i]--) 
        ret = (ret*i)%MOD;
    return ret%MOD;
}
ll solve(ll x) {
    if(f[x]) return f[x]; 
    ll ret;
    ret = mi(mi((ll)2, x)-1, m);
    for(int i = 1; i <= x; i++) 
        ret = (ret-C(x, i) * solve(x-(ll)i)+MOD)%MOD;
    f[x] = ret;
    return ret%MOD;
}
int main() {
    cin >> n >> m;
    cout << solve(n) << endl;
    return 0;
}

Problem 721: INTSUB

Time Limit: 1000 ms Memory Limit: 255999 KB

Problem Description
给出一个大小为2n的集合X=｛1，2，3，…，2n-1，2n｝，问该集合中有趣子集的数目，答案mod1e9+7。

x的有趣子集定义为，该子集中至少有两个数，a和b，b是a的倍数且a是集合中最小的元素。

Input
输入一个整数n

40% 1<=n<=50

100% 1<=n<=1000

Output
输出子集数目mod1e9+7的结果

Sample Input
3
Sample Output
47

分析

b是a的倍数
b有【n/a】个
所以答案为

∑i=1n2n/i−1∗2n−i−n/i+1

完整AC代码

其实不用快速幂

#include<cstdio>
#include<iostream>
#define MOD 1000000007
#define N 1010
using namespace std;
typedef long long LL;
int n,p;
LL a[2*N],ans,q;
LL Qpow(int q){
    LL p=2,ans=1;
    while(q!=0){
        if(q%2)ans=(ans*p)%MOD;
        p=(p*p)%MOD;
        q/=2;
    }
    return ans;
}
int main(){
    scanf("%d",&n);
    n*=2;
//  for(int i=0;i<=n;i++)a[i]=Qpow(i);
//  for(int i=0;i<=n;i++)cout<<a[i]<<'\n';
    for(int i=1;i<=n/2;i++){
        p=(n/i)-1;
        if(!a[p])a[p]=Qpow(p);
        if(!a[n-i-p])a[n-i-p]=Qpow(n-i-p);
        q=((a[p]-1)*a[n-i-p])%MOD;
        ans=(ans+q)%MOD;
//      cout<<ans<<'\n';
    }
    cout<<ans;
}

Problem 722: MINSUB

Time Limit: 3000 ms Memory Limit: 256000 KB

Problem Description
给一个n*m矩阵M，和一个整数k，对于矩阵M的子矩阵M’，定义min(M’)为M’矩阵中元素的最小值。

我们需要找出这样一个子矩阵，该矩阵的面积至少为K，且min(M’)最大化。在min(M’)相同的情况下，面积越大越好。面积的定义为该矩阵的行数*列数。

Input
第一行两个整数n和m，表示矩阵的大小。

   30%                1<=n,m<=50

   100%              1<=n,m<=1000     1<=k<=n*m

Output
输出两个值，第一个是最大的min(M’)，第二个是该矩阵的最大面积。

Sample Input
3 3 2
1 2 3
4 5 6
7 8 9
Sample Output
8 2

分析

二分矩阵元素大小最小值
把矩阵转换为0-1矩阵
扫描线用单调队列（栈）统计比当前小的元素的位置
反过来再扫一边
统计答案
继续二分
时限3s

完整AC程序(3s,stl栈)

#include<cstdio>
#include<iostream>
#include<stack>
#define N 1010
using namespace std;
struct Node{
    int v,w;
};
stack<Node>s;
int n,m,a[N][N],f[N][N],b[N],l,r,mid,c[N],k,ans,num,w;
int Solve(int x){
    int ans=0;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            if(a[i][j]>=x)f[i][j]=f[i-1][j]+1;
            else f[i][j]=0;
    for(int i=1;i<=n;i++){
        while(!s.empty()) s.pop();
        for(int j=1;j<=m;j++){
            while(!s.empty()&&s.top().v>f[i][j]) s.pop();
            if(!s.empty()){
                if(s.top().v==f[i][j])b[j]=b[s.top().w];
                else b[j]=s.top().w+1;
            }
            else b[j]=j;
            s.push((Node){f[i][j],j});
        }
        while(!s.empty()) s.pop();
        for(int j=m;j>=1;j--){
            while(!s.empty()&&s.top().v>f[i][j]) s.pop();
            if(!s.empty()){
                if(s.top().v==f[i][j])c[j]=b[s.top().w];
                else c[j]=s.top().w-1;
            }
            else c[j]=j;
            s.push((Node){f[i][j],j});
        }
        for(int j=1;j<=m;j++)
            ans=max(ans,(c[j]-b[j]+1)*f[i][j]);
    }
    return ans;
}
int main(){
    freopen("data.txt","r",stdin);
    scanf("%d%d%d",&n,&m,&k);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            scanf("%d",&a[i][j]);
            r=max(r,a[i][j]);
        }
    l=0;
    while(l<r){
        int mid=(l+r)/2;
        if((w=Solve(mid))>=k){
            if(mid>ans){
                ans=mid;
                num=w;
            }
            l=mid+1;
        }
        else r=mid;
    }
    if((w=Solve(l))>=k){
        if(l>ans){
            ans=l;
            num=w;
        }
    }
    printf("%d %d\n",ans,num);
}

标程（二分经我优化，可以跑进2s，类似kmp自配函数的单调队列）

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <ctime>
#include <iostream>
#include <algorithm>

using namespace std;

int n, m, a[1001][1001], test, k, f[1001][1001], U[1001], D[1001],ma;

int calc(int now) {
    memset(f, 0, sizeof(f));
    for (int i = 1; i <= n; i++)
        for (int j = 1; j <= m; j++)
            if (a[i][j] >= now) f[i][j] = f[i][j - 1] + 1;
            else f[i][j] = 0;
    int area = 0;
    for (int i = 1; i <= m; i++)
    {
        for (int j = 1; j <= n; j++)
        {
            int now = j - 1;
            for (; now && f[j][i] <= f[now][i]; now = U[now]);
            U[j] = now;
        }
        for (int j = n; j; --j)
        {
            int now = j + 1;
            for (; now != n + 1 && f[j][i] <= f[now][i]; now = D[now]);
            D[j] = now;
        }
        for (int j = 1; j <= n; j++)
            area = max(area, (D[j] - U[j] - 1) * f[j][i]);
    }
    return area;
}

int main() {
        freopen("data.txt","r",stdin);
        scanf("%d%d%d", &n, &m, &k);
        for (int i = 1; i <= n; i++)
            for (int j = 1; j <= m; j++){
                scanf("%d", &a[i][j]);
                ma=max(ma,a[i][j]);
            }

        int Left = 0, Right = ma, Mid = (Left + Right) >> 1;
        for (; Left + 1 < Right; Mid = (Left + Right) >> 1)
            if (calc(Mid) >= k)
                Left = Mid;
            else
                Right = Mid;
        printf("%d %d\n", Left, calc(Left));
}
mc_dl
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
2020 第十一届蓝桥杯校内模拟赛c/c++ B组
12-22
7.单词hello 思路:第一次不是元音字母的时候，找到结束位置记录，并break掉，后面从上次结束的下一个位置开始，遍历找到第一次不是辅音的字母…如此四次，注意中间判断是否达到字符串长度 8.正整数序列 思路：这个题...
8.21 模拟赛
lunch__的博客
08-23 159
T1 这个题看起来比较简单但是模型转化还是比较难想的 考场上我大概想到把答案记为总和减去选定的nnn个互不相同的宽 但是没有跳出这个模型 题解的思想巧妙多了 把每个权值离散后看作一个点，每对限制关系记为一条边 题目就可以转化为在每一条边上选择一个端点，不可以重复选择 求最小权值和 由于答案保证有解 所以最后得到的多个联通块要么是树，要么是基环树 如果是基环树的话，环上的权值是一...
第十四届蓝桥杯第三期模拟赛 C/C++ B组 原题与详解
热门推荐
weixin_67596609的博客
03-07 2万+
本篇文章对第十四届蓝桥杯第三期模拟赛所有的题目进行了详细解析。原题+代码+详解，一定会让你理解到位，赶快进来学习一下吧！！！
python乒乓球比赛规则_python模拟体育竞技分析--采用乒乓球规则
weixin_39593247的博客
11-27 2495
采用乒乓球比赛规则‪‬‪‬‪‬‪‬‪‬‮‬‫‬‪‬‪‬‪‬‪‬‪‬‪‬‮‬‭‬‫‬‪‬‪‬‪‬‪‬‪‬‮‬‫‬‭‬‪‬‪‬‪‬‪‬‪‬‮‬‫‬‫‬‪‬‪‬‪‬‪‬‪‬‮‬‫‬‫‬‪‬‪‬‪‬‪‬‪‬‮‬‭‬‪‬a.一局比赛：‪‬‪‬‪‬‪‬‪‬‮‬‫‬‪‬‪‬‪‬‪‬‪‬‪‬‮‬‭‬‫‬‪‬‪‬‪‬‪‬‪‬‮‬‫‬‭‬‪‬‪‬‪‬‪‬‪‬‮‬‫‬‫‬‪‬‪‬‪‬‪‬‪‬‮‬‫‬‫‬‪‬‪‬‪‬‪‬...
2020蓝桥杯省赛模拟赛
王跃坤的博客
10-11 2993
目录第一题：第二题：第三题：第四题：第五题：第六题：第七题：第八题：第九题：第十题： 第一题： 全排列问题+hashSet去重 package A; import java.util.*; public class Main { /** * 问题描述 　　将LANQIAO中的字母重新排列，可以得到不同的单词，如LANQIAO、AAILNOQ等，注意这7个字母都要被用上，单词不一定有具体的英文意义。 　　请问，总共能排列如多少个不同的单词。 答案提交 　　这是一道结果填空的题，你只需要算出结果后提
蓝桥杯第十四届蓝桥杯模拟赛第三期考场应对攻略（C/C++）
m0_71934846的博客
03-12 2万+
蓝桥杯第十四届蓝桥杯模拟赛第三期考场应对攻略（C/C++）
Java 第十一届 蓝桥杯 省模拟赛 螺旋矩阵
01-20
　13 20 19 18 7 　12 11 10 9 8 输入格式 输入的第一行包含两个整数 n, m，分别表示螺旋矩阵的行数和列数。 　第二行包含两个整数 r, c，表示要求的行号和列号。 输出格式 输出一个整数，表示螺旋矩阵中第 r 行第 ...
Java 第十一届 蓝桥杯 省模拟赛 反倍数
12-21
反倍数 题目 问题描述 给定三个整数 a, b, c，如果一个整数既不是 a ...以下这些数满足要求：1, 5, 7, 11, 13, 17, 19, 23, 25, 29。 评测用例规模与约定 对于 40% 的评测用例，1 <= n <= 10000。 　对于 80% 的
乒乓球比赛
09-19
乒乓球单打比赛程序模拟。比赛的双方简称A和B，他们的比赛规则采用： 1， 比赛胜负采用7局4胜制：谁先取得4局胜利就赢得比赛； 2， 每局比分采用11分制： A、如果比赛双方谁获得11分并且领先对方2分或2分以上，则该...
oj_468 出去玩
ljj的博客
03-20 1314
题目描述Problem 468: 出去玩
【gfoj】Problem#787: 密室逃脱
ljj的博客
08-27 933
Problem 787: 密室逃脱 Time Limit: 2000 ms Memory Limit: 524288 KBProblem Description 小J被关在密室里！密室的构造是一个H×W的棋盘，即有H行W列，第i行第j列的房间用Ai,j 表示。若Ai,j =’#’，则该房间的门是锁上的；若Ai,j =’.’，则该房间是可以自由进出的；小J被困在Ai,j =’S’的地方，该房间
【三分套三分】传送带
ljj的博客
10-25 472
题目在一个2维平面上有两条传送带，每一条传送带可以看成是一条线段。两条传送带分别为线段AB和线段CD。FTD在AB上的移动速度为P，在CD上的移动速度为Q，在平面上的移动速度R。现在FTD想从A点走到D点，他想知道最少需要走多长时间在一个2维平面上有两条传送带，每一条传送带可以看成是一条线段。两条传送带分别为线段AB和线段CD。FTD在AB上的移动速度为P，在CD上的移动速度为Q，在平面上的移动速度
oj 作业7 运输装备
ljj的博客
03-23 435
运输装备
中转站
ljj的博客
04-17 367
中转站
Problem 802: 递增数列
ljj的博客
08-25 364
Problem 802: 递增数列 Time Limit: 1000 ms Memory Limit: 262144 KBProblem Description 给定一个包含n个元素的数组，每次你可以把一个元素+1或-1，操作后的数可以为负数，零，正数。你需要操作最小次数，使得该数组严格递增。Input第一行一个整数n表示数组大小 接下来一行n个数字表示数组元素40% 1 <= n <= 1
【线段树】【CDQ分治】回转寿司
ljj的博客
10-23 360
题目大意给你一个序列，求连续子序列和在【L，R】之间的方案数 N≤100000，|Ai|≤100000，0≤L, R≤109. 分析这是我打的第一题CDQ（太菜了） 我对这题印象很深刻 当时大家有各种做法 好像都是线段树？ 然后这时出现了一股清流 dhr的CDQ分治（orz dhr 好短啊） 然后愉悦的改完后就没管了很久以后yzx讲CDQ分治时想到我曾经打过一道先递归 再更新答案 然后
oj_10 蚯蚓
ljj的博客
04-17 342
luogu原题：蚯蚓
提高组联赛训练07.19
ljj的博客
07-19 322
题目链接点这里王国道路（图论）Spfa(883ms)#include<cstdio> #include<iostream> #include<queue> #define INF 2147483647 #define N 2010 using namespace std; typedef long long LL; LL ans; queue<int>Q; struct Edge{ int
【dp】【优化】Problem 791: 硬币
ljj的博客
08-27 310
Problem 791: 硬币Time Limit: 2000 ms Memory Limit: 524288 KBProblem Description 有X+Y+Z个人，从1到X+Y+Z编号。第i个人有Ai个金币，Bi个银币，Ci个铜币。小J想从X个人手中得到他们的全部金币，Y个人手中得到他们的全部银币，Z个人手中得到他们的全部铜币。一个人只会将他的一种硬币全部给小J。小J想知道她最多能得到多
matlab 模拟22人足球比赛
最新发布
05-30
下面是一个简单的 MATLAB 代码，用于模拟22人足球比赛： ``` % 定义球队和比赛时间 teamA = {'A1', 'A2', 'A3', 'A4', 'A5', 'A6', 'A7', 'A8', 'A9', 'A10', 'A11'}; teamB = {'B1', 'B2', 'B3', 'B4', 'B5', 'B6'...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值