2. 程式人生 > 實用技巧 >隨機化的正確開啟方式

隨機化的正確開啟方式

阿新 • • 發佈：2020-11-23

P4168 [Violet]蒲公英

目錄

題目

傳送門

思路

經典的線上求眾數問題

預處理

首先,離散化時絕對跑不掉的,設a為離散化後的序列,c為原序列,b為離散化輔助陣列

struct node {
	int dat , id;
}b[nn];
bool cmp(node a , node b){return a.dat < b.dat;}

	//離散化 
	for(int i = 1 ; i <= n ; i++)
		b[i].dat = a[i] , b[i].id = i;
	sort(b + 1 , b + n + 1 , cmp);
	int cnt_ = 0;
	for(int i = 1 ; i <= n ; i++) {
		if(b[i].dat != b[i - 1].dat) cnt_++;
		a[b[i].id] = cnt_;
	}

把序列a分成t塊(提前說明t約為3次根號下n,具體原因在時間複雜度中講),每段長度為len=n/t,然後是分塊的常規操作:設L,R表示每一塊的左右端點,pos表示每個點所屬分塊

	t = 0;
	while(t * t * t < n)++t;
	len = n / t;
	
	L[1] = 1 , R[1] = len;
	for(int i = 2 ; i <= t ; i++)
		L[i] = R[i - 1] + 1,
		R[i] = len * i;
	if(R[t] < n)
		++t , L[t] = R[t - 1] + 1 , R[t] = n;
	
	for(int i = 1 ; i <= t ; i++)
		for(int j = L[i] ; j <= R[i] ; j++)
			pos[j] = i;

此外,我們設cnt[i][j][k]表示數字k在第i個塊到第j個塊出現的次數,zs[i][j]表示第i個塊到到j個塊的眾數的下標(原諒我不知道眾數的英文就直接上拼音了)

	//預處理cnt 和 眾數 
	for(int i = 1 ; i <= t ; i++)
		for(int j = i ; j <= t ; j++)
			for(int k = L[i] ; k <= R[j] ; k++) {
				++cnt[i][j][a[k]];
				if(cnt[i][j][a[k]] > cnt[i][j][a[zs[i][j]]] || (cnt[i][j][a[k]] == cnt[i][j][a[zs[i][j]]] && a[k] < a[zs[i][j]]))
					zs[i][j] = k;
			}

對於每一個詢問

若l,r屬於同一塊,則直接暴力

		for(int i = l ; i <= r ; i++) {
			++tmp_cnt[a[i]];
			if(tmp_cnt[a[i]] > tmp_cnt[a[ans]] || (tmp_cnt[a[i]] == tmp_cnt[a[ans]] && a[i] < a[ans]))
				ans = i;
		}
		for(int i = l ; i <= r ; i++)
			--tmp_cnt[a[i]];//這裡直接減應該比memset快(我沒試過),memset是針對整個陣列(就是O(m*n)了),而此時r-l不超過len,是根號級別
		return c[ans];

對於其他情況:

p=pos[l],q=pos[r]

和分塊模板一樣,我們把[l,r]分為:開頭:[l,R[p]),中間:[L[p+1],R[q-1]],結尾:(L[q],r](這裡注意下括號的意義,有的是下標,有的是區間)

顯然,最終眾數出為塊p+1~q-1的眾數,或開頭,結尾兩段的數之中

因此,我們令ans=zs[p+1][q-1],然後在cnt[p+1][q-1][]的基礎上加上開頭,結尾兩段的數,直接統計答案即可

	++p , --q;//方便起見
	ans = zs[p][q];
	for(int i = l ; i <= R[p - 1] ; i++) {
		++cnt[p][q][a[i]];
		if(cnt[p][q][a[i]] > cnt[p][q][a[ans]] || (cnt[p][q][a[i]] == cnt[p][q][a[ans]] && a[i] < a[ans]))
			ans = i;
	}
	for(int i = L[q + 1] ; i <= r ; i++) {
		++cnt[p][q][a[i]];
		if(cnt[p][q][a[i]] > cnt[p][q][a[ans]] || (cnt[p][q][a[i]] == cnt[p][q][a[ans]] && a[i] < a[ans]))
			ans = i;
	}
	for(int i = l ; i <= R[p - 1] ; i++)--cnt[p][q][a[i]];//這兩行用於清除開頭,結尾的影響
	for(int i = L[q + 1] ; i <= r ; i++)--cnt[p][q][a[i]];
	
	return c[ans];

時空複雜度

時間複雜度(前半段為預處理,後半段為每一次查詢):

\[O(nt^2+m\frac{n}{t}) \]

題目已經給出m,n在一個數量級(最大資料),因此,我們考慮讓兩邊儘量平均,則有方程組

\[nt^2=m\frac{n}{t} \]

解得t等於三次根號下m,約等於三次根號下n

空間複雜度:

\[O(nt^2) \]

由於t才去到30~40,所以是可以接受的

程式碼(含對拍)

宣告:使用對拍檔案時需要關掉強制線上,即直接輸入l,r(詳情參考std.cpp)

第一次交時忘記了強制線上(聽取WA聲一片),曾一度懷疑人生:難道我暴力都寫錯了?

不然就一次AC了

話說樣例真水啊,沒開強制線上都能過

AC程式碼(tested.cpp)

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#define nn 40010
#define max_t 64+10
using namespace std;
int read() {
	int re = 0;
	bool sig = false;
	char c = getchar();
	while(c < '0' || c > '9') {
		if(c == '-')sig = true;
		c = getchar();
	}
	while(c >= '0' && c <= '9')
		re = (re << 1) + (re << 3) + c - '0',
		c = getchar();
	return sig ? -re : re;
}
int n , m;
int t , len;
int a[nn];
int c[nn]; 
int zs[max_t][max_t];
int L[max_t] , R[max_t];
int pos[nn];
unsigned short cnt[max_t][max_t][nn];

struct node {
	int dat , id;
}b[nn];
bool cmp(node a , node b){return a.dat < b.dat;}
void Init() {
	//離散化 
	for(int i = 1 ; i <= n ; i++)
		b[i].dat = a[i] , b[i].id = i;
	sort(b + 1 , b + n + 1 , cmp);
	int cnt_ = 0;
	for(int i = 1 ; i <= n ; i++) {
		if(b[i].dat != b[i - 1].dat) cnt_++;
		a[b[i].id] = cnt_;
	}
	//分塊 
	t = 0;
	while(t * t * t < n)++t;
	len = n / t;
	
	L[1] = 1 , R[1] = len;
	for(int i = 2 ; i <= t ; i++)
		L[i] = R[i - 1] + 1,
		R[i] = len * i;
	if(R[t] < n)
		++t , L[t] = R[t - 1] + 1 , R[t] = n;
	
	for(int i = 1 ; i <= t ; i++)
		for(int j = L[i] ; j <= R[i] ; j++)
			pos[j] = i;
	//預處理cnt 和 眾數 
	for(int i = 1 ; i <= t ; i++)
		for(int j = i ; j <= t ; j++)
			for(int k = L[i] ; k <= R[j] ; k++) {
				++cnt[i][j][a[k]];
				if(cnt[i][j][a[k]] > cnt[i][j][a[zs[i][j]]] || (cnt[i][j][a[k]] == cnt[i][j][a[zs[i][j]]] && a[k] < a[zs[i][j]]))
					zs[i][j] = k;
			}
}
int tmp_cnt[nn];
int query(int l , int r) {
	int p = pos[l] , q = pos[r];
	int ans = 0;
	if(p == q) {
		for(int i = l ; i <= r ; i++) {
			++tmp_cnt[a[i]];
			if(tmp_cnt[a[i]] > tmp_cnt[a[ans]] || (tmp_cnt[a[i]] == tmp_cnt[a[ans]] && a[i] < a[ans]))
				ans = i;
		}
		for(int i = l ; i <= r ; i++)
			--tmp_cnt[a[i]];
		
		return c[ans];
	} 
	++p , --q;
	ans = zs[p][q];
	for(int i = l ; i <= R[p - 1] ; i++) {
		++cnt[p][q][a[i]];
		if(cnt[p][q][a[i]] > cnt[p][q][a[ans]] || (cnt[p][q][a[i]] == cnt[p][q][a[ans]] && a[i] < a[ans]))
			ans = i;
	}
	for(int i = L[q + 1] ; i <= r ; i++) {
		++cnt[p][q][a[i]];
		if(cnt[p][q][a[i]] > cnt[p][q][a[ans]] || (cnt[p][q][a[i]] == cnt[p][q][a[ans]] && a[i] < a[ans]))
			ans = i;
	}
	for(int i = l ; i <= R[p - 1] ; i++)--cnt[p][q][a[i]];
	for(int i = L[q + 1] ; i <= r ; i++)--cnt[p][q][a[i]];
	
	return c[ans];
} 
int main() {
	n = read(),	m = read();
	for(int i = 1 ; i <= n ; i++)
		c[i] = a[i] = read();
	Init();
	
	int lastans = 0;
	for(int i = 1 ; i <= m ; i++) {
		int l , r;
		l = read(),	r = read();
		l = ((l + lastans - 1) % n) + 1;
		r = ((r + lastans - 1) % n) + 1;
		if(l > r) {int tmp = l ; l = r ; r = tmp;}
//		query(l , r);
		printf("%d\n" , lastans = query(l , r));
	}
	return 0;
}

暴力(std.cpp)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define nn 40010
using namespace std;
int read() {
	int re = 0;
	bool sig = false;
	char c = getchar();
	while(c < '0' || c > '9') {
		if(c == '-')sig = true;
		c = getchar();
	}
	while(c >= '0' && c <= '9')
		re = (re << 1) + (re << 3) + c - '0',
		c = getchar();
	return sig ? -re : re;
}

int n , m ;
int a[nn] , c[nn];
struct node {
	int dat , id;
}b[nn];
bool cmp(node a , node b){return a.dat < b.dat;}

int cnt[nn];
int main() {
	n = read(),	m = read();
	for(int i = 1 ; i <= n ; i++)
		a[i] = b[i].dat = c[i] = read(),	b[i].id = i;
	
	sort(b + 1 , b + n + 1 , cmp);
	int cnt_ = 0;
	for(int i = 1 ; i <= n ; i++) {
		if(b[i].dat != b[i - 1].dat) cnt_++;
		a[b[i].id] = cnt_;
	} 
	
	for(int i = 1 ; i <= m ; i++) {
		int l =read() , r = read();
		memset(cnt , 0 , sizeof(cnt));
		int ans = 0;
		for(int j = l ; j <= r ; j++) {
			cnt[a[j]]++;
			if(cnt[a[j]] > cnt[a[ans]] || (cnt[a[j]] == cnt[a[ans]] && a[j] < a[ans]))
				ans = j;
		}
		printf("%d\n" , c[ans]);
	}
	return 0;
}

隨機資料(random.cpp)

#include <bits/stdc++.h>
using namespace std;
int random(int r , int l = 1) {
	return (l == r ? l : ((long long)rand() * rand() % (r - l) + l ));
}
int main() {
	srand((unsigned)time(0));
	
	int n = 40000 , m = 50000;
	printf("%d %d\n" , n , m);
	for(int i = 1 ; i <= n ; i++)
		printf("%d " , random(1000000000));
	putchar('\n');
	
	for(int i = 1 ; i <= m ; i++) {
		int l = random(n) , r = random(n , l);
		printf("%d %d\n" , l , r); 
	} 
		
	return 0;
}

對拍控制(compare.cpp)

將continue刪掉,std的雙斜槓去掉,即開啟std和tested的對拍

#include <bits/stdc++.h>
using namespace std;
int main() {
	
	while(true) {
		system("random.exe > input.txt");
		puts("random");
		
//		system("std.exe < input.txt > output1.txt");
//		puts("std");
		int t = clock();
		if(system("tested.exe < input.txt > output2.txt") != 0) {//檢驗執行時錯誤,記得return 0
			cout << "RE";
			return 0;
		}
		puts("tested");
		printf(">time:%d\n" , clock() - t);
		
		continue;
		if(system("fc output1.txt output2.txt")) {
			cout << "WA";
			system("start input.txt");
			return 0;
		}
		
	}
	return 0;
}

相關推薦

隨機化的正確開啟方式

前言 最近在考場上發現有的題目暴力搜尋並不是很好打 而且很耗時間 效率低下

詳解B+樹及其正確開啟方式

前言 hello，小可愛們，繼上篇長文的更新，我又又又來了。 前面我們長大了InnoDB資料頁的7個組成部分，各個資料頁組成了一個雙向連結串列，而每個資料頁中的記錄按照主鍵從小到大的順序組成一個單連結串列，每個資料

雲開發 VSCode 外掛 Cloudbase Toolkit 的正確開啟方式

什麼是 Cloudbase Toolkit Tencent CloudBase Toolkit 是雲開發的 VS Code（Visual Studio Code）外掛。該外掛可以讓您更好地在本地進行雲開發專案開發和程式碼除錯，並且輕鬆將專案部署到雲端。

雲開發 VSCode 外掛 Cloudbase Toolkit 的正確開啟方式及應用場景分析

什麼是 Cloudbase Toolkit Tencent CloudBase Toolkit 是雲開發的 VS Code（Visual Studio Code）外掛。該外掛可以讓您更好地在本地進行雲開發專案開發和程式碼除錯，並且輕鬆將專案部署到雲端。

基於Redis分散式鎖的正確開啟方式

分散式鎖是在分散式環境下（多個JVM程式）控制多個客戶端對某一資源的同步訪問的一種實現，與之相對應的是執行緒鎖，執行緒鎖控制的是同一個JVM程式內多個執行緒之間的同步。分散式鎖的一般實現方法是在應用伺服器之

rclone 雲盤同步工具的正確開啟方式

Rclone 是一款的命令列工具，支援在不同物件儲存、網盤間同步、上傳、下載資料。

idea克隆非Maven專案的正確開啟方式

下載 clone專案，可以直接下載zip，也可以使用idea的新建方式 輸入URL 或者下載zip，然後用idea開啟

python中四捨五入的正確開啟方式

round()函式 （注意：下面的我也不清楚是否正確，我只是發表一下我的觀點） 對於簡單的舍入運算，使用內建的 round(value,ndigits) 函式即可

golang的正確開啟方式【第四章 golang特性】

技術標籤：go 一，map 定義 result := make(map[string]interface{}) 使用 result["total"] = total

《生化危機8》二週目正確開啟方式及建議

《生化危機8》中很多小夥伴覺得開啟二週目就是重新開一個新檔進行第二遍的遊戲，其實不然，本作是有正確的二週目進入方法的，那麼現在為大家帶來“奧塞羅特1981”分享的《生化危機8》二週目正確開啟方式及建議，希望

vue列印正確開啟方式

1.首先下載安裝包 npm install vue-print-nb --save 2.在main.js中引用 import Print from \'vue-print-nb\'

Vue配置環境變數的正確開啟方式

目錄第一 配置package.on第二 配置介紹 簡單說明 看詳情點選上面tps第三 在根目錄新建檔案第四第五總結:第一 配置package.json

跟著麓安一起上分，新年廣場正確開啟方式

不要在匹配到新年廣場時躡手躡腳，如何拒絕在新年廣場翻車，快點跟著麓安小課堂一起學習這些上分小技巧吧！在我們觀看麓安精彩的上分表演時，你會發現麓安有很多點位都值得去學習，接下來我就帶大家一起來看一看！

B站學習（存檔）-python學習中Anaconda和Pycharm的正確開啟方式

虛擬環境 建立虛擬環境conda create -n name python=3.8 啟用虛擬環境conda activate name conda 安裝

Git 工作流的正確開啟方式

前言 一直在使用git做版本控制，也一直工作很順利，直到和別人發生衝突的時候。這才注意到git 工作流並不是那麼簡單。比如，之前遇到的清理歷史。百度到的資料很多，重複性也很多，但實踐性操作很少，我很難直接理解

【Note】 tqdm 簡介及正確的開啟方式

tqdm簡介及正確開啟方式 什麼是tqdm？ tqdm是一個快速的，易擴充套件的進度條提示模組，官方網站

開啟Python列表解析式的正確使用方式

Python 是一種極其多樣化和強大的程式語言！當需要解決一個問題時，它有著不同的方法。在本文中，將會展示列表解析式

透徹！這才是正確的分散式事務的開啟方式

尊重原創版權: https://www.gewuweb.com/hot/18048.html 透徹！這才是正確的分散式事務的開啟方式

vue 相鄰自定義元件渲染錯誤正確的開啟方式

話不多說看問題： 當封裝自定義元件時例如（自定義下拉列表）兩個相同的元件在多次v-if變化時偶爾會發生渲染錯誤，明明賦值正確但是組建中的ajax方法可能返回的資料亂掉，或者其他神邏輯錯誤。

一篇文章帶你解讀redis分散式鎖的發展史和正確實現方式

來源：點我達技術 http://tech.dianwoda.com/ 前言 近兩年來微服務變得越來越熱門，越來越多的應用部署在分散式環境中，在分散式環境中，資料一致性是一直以來需要關注並且去解決的問題，分散式鎖也就成為了一種廣泛