[递推式求解、多指针、前缀和]XATU第七届算法大赛

文章目录

    • 赛后感想
    • 1、2题略
    • 3 乐华水上梦幻乐园之拖鞋失踪案[概率论、递推式求解]
      • 题意
      • 分析
      • 代码
    • 4、5、6、7、8题略
    • 9 Ugly Number[多指针]
      • 题意
      • 做法1[BFS with heap]
      • 做法2[多指针]
    • 10 怪兽大学
      • 题意
      • 分析
    • 总结

赛后感想

总结来看的感觉只有一个吧,就是痛心于我校的算法教学、竞赛水平的滞后。题目不难,压轴题更是我签完到之后直接A掉了,并不是我有多强,而是我校同仁总体水平确实有差距——这种差距不在于和985、211相比,而是和一些本地的、同类院校相比,差距也不可谓不大。不过幸而今年有小川、一鸣等同学的热心,希望XATU-PST能够越办越好吧。

以下过于简单的题目略过。

1、2题略

3 乐华水上梦幻乐园之拖鞋失踪案[概率论、递推式求解]

题意

n n n个人依次找拖鞋,第 1 1 1位会随机挑选一位幸运观众,此后每一位将:

  • 拿走自己的拖鞋,若自己的拖鞋还在
  • 随机挑选一位幸运观众,若自己的拖鞋不在了

问第 n n n位幸运儿拿到自己拖鞋的概率是?

分析

乍一看像是错位排列,但其实关系不大。本场比赛离AK就差这一道题。赛后看题解…好像确实不难。不过这里要提一句,网上的题解全都犯了一个不影响结论的错误,以下按照正确的讲。

f ( n ) f(n) f(n)代表第1个人坐在座位 C ( 1 ) C(1) C(1)后第 n n n个人坐到自己座位的概率,显然有以下公式:
f ( n ) = { 1 , if  C ( 1 ) = 1 0 , if  C ( 1 ) = n g ( C ( 1 ) ) , else f(n)= \begin{cases} 1, &\text{if $C(1)=1$}\\ 0, &\text{if $C(1)=n$}&\\ g(C(1)), &\text{else} \end{cases} f(n)=1,0,g(C(1)),if C(1)=1if C(1)=nelse
含义为:若1号坐到自己的位置,那么之后的人(包括 n n n)一定能坐到自己的位置;若坐到位置 n n n,则 n n n号一定坐不到自己的位置,否则问题转化为被1号玩家抢占位置的 C ( 1 ) C(1) C(1)号玩家找位置。网上的所有题解这里写的都是
g ( C ( 1 ) ) , else g(C(1)), \text{else} g(C(1)),else
但实际上除了第一位玩家,其他人只要需要开始找座位,就不存在找到自己座位的可能性,因此有:
g ( x ) = { 0 , if  C ( x ) = n g ( C ( x ) ) , else g(x)= \begin{cases} 0, &\text{if $C(x)=n$}&\\ g(C(x)), &\text{else} \end{cases} g(x)={0,g(C(x)),if C(x)=nelse
但是由于对递推式求解时这一项1能够消掉,所以这个错误并不影响答案的正确性。

由等概率原理,我们能写出 f ( x ) f(x) f(x) g ( x ) g(x) g(x)的显式公式:
f ( n ) = 1 + ∑ i = 2 n − 1 g ( i ) + 0 n f(n)=\frac{1+\sum_{i=2}^{n-1}g(i)+0}{n} f(n)=n1+i=2n1g(i)+0
g ( x ) = ∑ i = 2 x − 1 g ( i ) x g(x)=\frac{\sum_{i=2}^{x-1}g(i)}{x} g(x)=xi=2x1g(i)
于是有
g ( n ) = ∑ i = 2 n − 1 g ( i ) n g(n)=\frac{\sum_{i=2}^{n-1}g(i)}{n} g(n)=ni=2n1g(i)
g ( n − 1 ) = ∑ i = 2 n − 2 g ( i ) n − 1 g(n-1)=\frac{\sum_{i=2}^{n-2}g(i)}{n-1} g(n1)=n1i=2n2g(i)
两式做差得到
n g ( n ) − ( n − 1 ) g ( n ) = g ( n − 1 ) ng(n)-(n-1)g(n)=g(n-1) ng(n)(n1)g(n)=g(n1)
从而
g ( n ) = g ( n − 1 ) g(n)=g(n-1) g(n)=g(n1)
得到了这样一个意想不到的结果。
容易知道 g ( 2 ) = 0.5 g(2)=0.5 g(2)=0.5并且有特例 g ( 1 ) = 1 g(1)=1 g(1)=1,从而有答案
g ( x ) = { 1 , if  x = 1 0.5 , if  x > 1 g(x)= \begin{cases} 1, &\text{if $x=1$}&\\ 0.5, &\text{if $x>1$} \end{cases} g(x)={1,0.5,if x=1if x>1

代码

略。

4、5、6、7、8题略

9 Ugly Number[多指针]

题意

质因数只包含2、3、5的正整数称之为Ugly Number,求第 N N N个ugly number。

做法1[BFS with heap]

先交了一发暴力,果断TLE了。那么很快想到写BFS with heap,每次从堆里取出最小的数字为第++ i i i个ugly number,然后放入 i ∗ 2 i*2 i2 i ∗ 3 i*3 i3 i ∗ 5 i*5 i5。由于每次放入和取出heap的元素量同阶,因此复杂度是 O ( N l o g N ) O(NlogN) O(NlogN)

#include<iostream>
#include<algorithm>
#include<vector>
#include<cstdio>
#include<cstring>
#include<queue>
#include<set>
using namespace std;
typedef long long ll;
vector<ll> list;
set<ll> vis;
void gen(int cnt)
{
	priority_queue<ll, vector<ll>, greater<ll> > heap;
	heap.push(1);
	while (list.size() < cnt)
	{
		ll now = heap.top();
		heap.pop();
		if (vis.find(now) != vis.end())
			continue;
		vis.insert(now);
		list.push_back(now);
		heap.push(now*2);
		heap.push(now*3);
		heap.push(now*5);
	}
}
int n;
vector<ll> ask;
int maxx = 0;
int main()
{
	cin>>n;
	for (int i = 1; i <= n; i++)
	{
		ll x;
		cin>>x;
		ask.push_back(x);
		if (x>maxx)
			maxx = x;
	}
	gen(2000);
	for (auto v : ask)
		if (!v)
		cout<<0<<' ';
		else
		cout<<list[v-1]<<' ';
	return 0;
}

做法2[多指针]

赛后看到题解,了解到是来自力扣的一道原题,标准做法是多指针,时间复杂度 O ( N ) O(N) O(N),看完代码表示确实学到了。

刚听到可以用三指针的时候没太明白怎么保证每时刻找到最小值的时候还能跑遍各自乘法群的,但仔细看一下代码就明白了,以下是代码:

#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;

int main()
{
    int n;
    cin >> n;
    vector<int> a(n);
    for (int i = 0; i < n; i++)
        cin >> a[i];

    auto m = *max_element(a.begin(), a.end());
    vector<int> dp(m + 1);
    dp[1] = 1;
    int p2 = 1, p3 = 1, p5 = 1;
    for (int i = 2; i <= m; i++)
    {
        int num2 = dp[p2] * 2, num3 = dp[p3] * 3, num5 = dp[p5] * 5;
        dp[i] = min(min(num2, num3), num5);
        if (dp[i] == num2)
            p2++;
        if (dp[i] == num3)
            p3++;
        if (dp[i] == num5)
            p5++;
    }

    for (int i : a)
        cout << dp[i] << ' ';
    return 0;
}

与我的第一感不同, p 2 p_2 p2 p 3 p_3 p3 p 5 p_5 p5指向的是要被乘以2、3、5的原数字,而每次 p 2 p_2 p2 p 3 p_3 p3 p 5 p_5 p5只+1而非跳到新的数字上,所以跑遍正整数,从而 p 2 ∗ 2 p_2*2 p22 p 3 ∗ 3 p_3*3 p33 p 5 ∗ 5 p_5*5 p55跑遍丑数群。然后每时刻找到这三个产生数字中最小的即可。

这样看来,这种多指针问题适用于解决各种迭代产生某个群的问题,只要理解指针指原数字而非产生的数字,每次指针步进而非跳到新产生的数字即可。

10 怪兽大学

题意

给定一组数 A = { a 1 , a 2 , . . . , a n } A=\{a_1,a_2,...,a_n\} A={a1,a2,...,an},求 ∑ i ∑ j f l o o r ( a i a j ) \sum_i\sum_jfloor(\frac{a_i}{a_j}) ijfloor(ajai)

分析

据说又是力扣搬来的(都怪我不做力扣

做完带歪榜的不是签到题的签到题之后直接把这个所谓的压轴题A了。很裸的前缀和,有人说思路不好想,其实大概是这样:

向下取整 -> 某个区间对于结果的贡献相同 -> 离线区间计数 -> 前缀和

所以其实感觉还是很好想的,这个力扣困难题不要说放到ICPC了,NOIP普及组可能都是d1t1够不上的水准。代码就非常简单了:

#include<iostream>
#include<algorithm>
#include<vector>
#include<cstdio>
#include<cstring>
#include<set>
using namespace std;
int n;
vector<int> arr;
#define N 200005
int pre[N];
#define M 1000000007
int ans;
int range(int x, int y)
{
	return pre[y] - pre[x-1];
}
int main()
{
	cin>>n;
	for (int i  = 1; i <= n ;i++)
	{
		int tmp;
		cin>>tmp;
		arr.push_back(tmp);
	}
	for (auto v : arr)
		pre[v]++;
	for (int i = 1; i <= 200000; i++)
		pre[i] += pre[i-1];
	for (auto v : arr)
		for (int i = 1; i*v <= 100000; i++)
		{
			ans += i*range(i*v, (i+1)*v-1);
			ans %= M;
		}
	cout<<ans;
	return 0;
}

总结

一直以出题人的身份参与学校各项比赛,好久没有真正参赛了。一方面因为懒得出题,另一方面因为奖品比较好?(可惜被某人抢走了,啧),赛前还是很期待的。

不过要实话实说,赛中体验还是非常不好,因为遇到了非常多的题意不清、数据规模没给、测试数据与题给范围不符,甚至还有赛后才被我找出来的数据错误。挺拉跨的。

如果没有被错误的数据浪费的1个小时的话,唯一一道没有A掉的“阴间题”拖鞋失踪案应该有机会打表找到这个呆逼规律的。但总归还是自己太菜了,不会推公式,这一块还需要加强。

不过这道题还是非常不错的,总归是这次比赛里最像正经题目的了——虽然还是力扣扒的(你们怎么这么爱扒力扣啊喂!

跟红细胞一起进步了这么久,希望没有我在的日子里学弟们还能继续发扬光大吧。

热门文章

暂无图片
编程学习 ·

Java输出数组的内容

Java输出数组的内容_一万个小时-CSDN博客_java打印数组内容1. 输出内容最常见的方式// List<String>类型的列表List<String> list new ArrayList<String>();list.add("First");list.add("Second");list.add("Third");list.ad…
暂无图片
编程学习 ·

母螳螂的“魅惑之术”

在它们对大蝗虫发起进攻的时候&#xff0c;我认认真真地观察了一次&#xff0c;因为它们突然像触电一样浑身痉挛起来&#xff0c;警觉地面对限前这个大家伙&#xff0c;然后放下自己优雅的身段和祈祷的双手&#xff0c;摆出了一个可怕的姿势。我被眼前的一幕吓到了&#xff0c;…
暂无图片
编程学习 ·

疯狂填词 mad_libs 第9章9.9.2

#win7 python3.7.0 import os,reos.chdir(d:\documents\program_language) file1open(.\疯狂填词_d9z9d2_r.txt) file2open(.\疯狂填词_d9z9d2_w.txt,w) words[ADJECTIVE,NOUN,VERB,NOUN] str1file1.read()#方法1 for word in words :word_replaceinput(fEnter a {word} :)str1…
暂无图片
编程学习 ·

HBASE 高可用

为了保证HBASE是高可用的,所依赖的HDFS和zookeeper也要是高可用的. 通过参数hbase.rootdir指定了连接到Hadoop的地址,mycluster表示为Hadoop的集群. HBASE本身的高可用很简单,只要在一个健康的集群其他节点通过命令 hbase-daemon.sh start master启动一个Hmaster进程,这个Hmast…
暂无图片
编程学习 ·

js事件操作语法

一、事件的绑定语法 语法形式1 事件监听 标签对象.addEventListener(click,function(){}); 语法形式2 on语法绑定 标签对象.onclick function(){} on语法是通过 等于赋值绑定的事件处理函数 , 等于赋值本质上执行的是覆盖赋值,后赋值的数据会覆盖之前存储的数据,也就是on…
暂无图片
编程学习 ·

Photoshop插件--晕影动态--选区--脚本开发--PS插件

文章目录1.插件界面2.关键代码2.1 选区2.2 动态晕影3.作者寄语PS是一款栅格图像编辑软件&#xff0c;具有许多强大的功能&#xff0c;本文演示如何通过脚本实现晕影动态和选区相关功能&#xff0c;展示从互联网收集而来的一个小插件&#xff0c;供大家学习交流&#xff0c;请勿…
暂无图片
编程学习 ·

vs LNK1104 无法打开文件“xxx.obj”

写在前面&#xff1a; 向大家推荐两本新书&#xff0c;《深度学习计算机视觉实战》和《学习OpenCV4&#xff1a;基于Python的算法实战》。 《深度学习计算机视觉实战》讲了计算机视觉理论基础&#xff0c;讲了案例项目&#xff0c;讲了模型部署&#xff0c;这些项目学会之后可以…
暂无图片
编程学习 ·

工业元宇宙的定义与实施路线图

工业元宇宙的定义与实施路线图 李正海 1 工业元宇宙 给大家做一个关于工业元宇宙的定义。对于工业&#xff0c;从设计的角度来讲&#xff0c;现在的设计人员已经做到了普遍的三维设计&#xff0c;但是进入元宇宙时代&#xff0c;就不仅仅只是三维设计了&#xff0c;我们的目…
暂无图片
编程学习 ·

【leectode 2022.1.15】完成一半题目

有 N 位扣友参加了微软与力扣举办了「以扣会友」线下活动。主办方提供了 2*N 道题目&#xff0c;整型数组 questions 中每个数字对应了每道题目所涉及的知识点类型。 若每位扣友选择不同的一题&#xff0c;请返回被选的 N 道题目至少包含多少种知识点类型。 示例 1&#xff1a…
暂无图片
编程学习 ·

js 面试题总结

一、js原型与原型链 1. prototype 每个函数都有一个prototype属性&#xff0c;被称为显示原型 2._ _proto_ _ 每个实例对象都会有_ _proto_ _属性,其被称为隐式原型 每一个实例对象的隐式原型_ _proto_ _属性指向自身构造函数的显式原型prototype 3. constructor 每个prot…
暂无图片
编程学习 ·

java练习代码

打印自定义行数的空心菱形练习代码如下 import java.util.Scanner; public class daYinLengXing{public static void main(String[] args) {System.out.println("请输入行数");Scanner myScanner new Scanner(System.in);int g myScanner.nextInt();int num g%2;//…
暂无图片
编程学习 ·

RocketMQ-什么是死信队列?怎么解决

目录 什么是死信队列 死信队列的特征 死信消息的处理 什么是死信队列 当一条消息初次消费失败&#xff0c;消息队列会自动进行消费重试&#xff1b;达到最大重试次数后&#xff0c;若消费依然失败&#xff0c;则表明消费者在正常情况下无法正确地消费该消息&#xff0c;此时…
暂无图片
编程学习 ·

项目 cg day04

第4章 lua、Canal实现广告缓存 学习目标 Lua介绍 Lua语法 输出、变量定义、数据类型、流程控制(if..)、循环操作、函数、表(数组)、模块OpenResty介绍(理解配置) 封装了Nginx&#xff0c;并且提供了Lua扩展&#xff0c;大大提升了Nginx对并发处理的能&#xff0c;10K-1000K Lu…
暂无图片
编程学习 ·

输出三角形

#include <stdio.h> int main() { int i,j; for(i0;i<5;i) { for(j0;j<i;j) { printf("*"); } printf("\n"); } }
暂无图片
编程学习 ·

stm32的BOOTLOADER学习1

序言 最近计划学习stm32的BOOTLOADER学习,把学习过程记录下来 因为现在网上STM32C8T6还是比较贵的,根据我的需求flash空间小一些也可以,所以我决定使用stm32c6t6.这个芯片的空间是32kb的。 #熟悉芯片内部的空间地址 1、flash ROM&#xff1a; 大小32KB&#xff0c;范围&#xf…
暂无图片
编程学习 ·

通过awk和shell来限制IP多次访问之学不会你打死我

学不会你打死我 今天我们用shell脚本&#xff0c;awk工具来分析日志来判断是否存在扫描器来进行破解网站密码——限制访问次数过多的IP地址&#xff0c;通过Iptables来进行限制。代码在末尾 首先我们要先查看日志的格式&#xff0c;分析出我们需要筛选的内容&#xff0c;日志…
暂无图片
编程学习 ·

Python - 如何像程序员一样思考

在为计算机编写程序之前&#xff0c;您必须学会如何像程序员一样思考。学习像程序员一样思考对任何学生都很有价值。以下步骤可帮助任何人学习编码并了解计算机科学的价值——即使他们不打算成为计算机科学家。 顾名思义&#xff0c;Python经常被想要学习编程的人用作第一语言…
暂无图片
编程学习 ·

蓝桥杯python-数字三角形

问题描述 虽然我前后用了三种做法&#xff0c;但是我发现只有“优化思路_1”可以通过蓝桥杯官网中的测评&#xff0c;但是如果用c/c的话&#xff0c;每个都通得过&#xff0c;足以可见python的效率之低&#xff08;但耐不住人家好用啊&#xff08;哭笑&#xff09;&#xff09…