laok 发布的文章
模板---三分法
原题为luogu3382，传送门
给出一个N次函数$f$，保证在范围$[l,r]$内存在一点$x$，使得$[l,x]$上单调增，$[x,r]$上单调减。试求出$x$的值。
UOJ12 猜数
问题描述
这一天，小Y、小D、小C正在愉快地玩耍。
小Y是个数学家，他一拍脑袋冒出了一个神奇的完全平方数 $n$。
小D是个机灵鬼，很快从小Y嘴里套出了 $n$ 的值。然后在脑内把 $n$ 写成了 $a×b$的形式。其中 $a,b$都是正整数。
小C是个八卦狂，他发现小D从小Y那里获知了神奇的东西，于是死缠烂打追问小D。最后小D说道：“我可以告诉你正整数 $g$ 和 $l$ 的值，我保证 $ab=gl=n$ 且 $a,b$ 都是$g$的倍数。但是 $a,b$ 我可不能告诉你。”
这可急坏了小C。他决定退而求其次，找出$a+b$的最小值和最大值。请你帮帮他吧！
[NOI2014&UOJ#6]随机数生成器
这道题不要用C++11，否则luogu上T好几个点
题面见题目传送门
首先，因为$N,M \leq 5000$,所以$ O(nm)$ 的求$T_i$ 的纯模拟算法时限完全没有问题。
接着就是求路径序列了：千万注意要好好读题
“字典序最小” “从小到大排序”说明什么？要贪心做！
那么我们如何贪心呢？
直观地想到：从小到大考虑能不能走到它。
但是直接判断的话复杂度可是会超时的啊！
假设i在第$X_i$行，$Y_i$列。
那么显然：
当i和j都可以走，当$X_i < X_j$，$Y_i \leq Y_j$,利用这个性质，我们可以简化判断。
假设i选了，则对于所有j,当$X_j < X_i$，则$Y_j \geq Y_i$ ，反之同理。于是写出代码：
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define up(a,b,c) for (int a(b),end##a(c);a<=end##a;++a)
#define down(a,c,b) for (int a(b),end##a(c);a>=end##a;--a)
#define uup(a,c,b) for (int a(b),end##a(c);a<=end##a;++a)
#define udown(a,b,c) for (int a(b),end##a(c);a>=end##a;--a)
using namespace std;
int read(){
char s;
int k=0,base=1;
while((s=getchar())!='-'&&s!=EOF&&!(s>='0'&&s<='9'));
if(s==EOF)exit(0);
if(s=='-')base=-1,s=getchar();
while(s>='0'&&s<='9'){
k=k*10+(s-'0');
s=getchar();
}
return k*base;
}
void write(int x){
if(x<0){
putchar('-');
write(-x);
}
else{
if(x/10)write(x/10);
putchar(x%10+'0');
}
}
int t[25050000],ft[25050000],lj[5050],rj[5050],n,m,q;
void myswap(int &a,int &b){
int t=a;a=b;b=t;
}
int main(){
long long x=read(),a=read(),b=read(),c=read(),d=read();
n=read(),m=read(),q=read();
int nm=n*m;
up(i,1,nm){
x=(x*x*a+x*b+c)%d;
t[i]=i;
myswap(t[i],t[(int)x%i+1]);
}
up(i,1,q){
myswap(t[read()],t[read()]);
}
up(i,1,nm){
ft[t[i]]=i;
}
up(i,1,n)lj[i]=1,rj[i]=m;
int tt=n+m-1,h,l;
up(i,1,nm){
h=(ft[i]-1)/m+1;l=(ft[i]-1)%m+1;
if(lj[h]<=l&&rj[h]>=l){//可以选
printf("%d%c",i,--tt?' ':'\n');
if(!tt)return 0;
int k=h-1;
while(k){
rj[k]=rj[k]>l?l:rj[k];//根据性质
--k;
}
k=h+1;
while(k<=n){
lj[k]=lj[k]<l?l:lj[k];//根据性质
++k;
}
}
}
return 0;
}