C.Discrete Acceleration
发现自己傻了,这题明明可以直接O(n)模拟,我却写了个二分,但其实我的这个二分还是O(n)的复杂度,因为之前O(n)处理了左右端点到每个点的时间。
我的傻逼写法:
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<map>
#include<queue>
#include<cmath>
using namespace std;
template<class T>inline void read(T &x){x=0;char o,f=1;while(o=getchar(),o<48)if(o==45)f=-f;do x=(x<<3)+(x<<1)+(o^48);while(o=getchar(),o>47);x*=f;}
#define ll long long
#define ull unsigned long long
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define repb(i,a,b) for(int i=(a);i>=b;i--)
#define log(x) (31-__builtin_clz(x))
#define INF 0x3f3f3f3f
ll gcd(ll a,ll b){ while(b^=a^=b^=a%=b); return a; }
//#define INF 0x7fffffff
const int MAXN = 1e5+5;
int a[MAXN];
double ts[MAXN],tn[MAXN];// 记录正着反着到达每个点的时间
int n,len;
void solve(){
scanf("%d%d",&n,&len);
rep(i,1,n) cin>>a[i];// 保证输入升序
int v = 1;//速度
double t = 0;
a[0] = 0;a[n+1] = len;
ts[0] = tn[n+1] = 0;
rep(i,1,n){
t += (double)(a[i]-a[i-1])/v;
ts[i] = t;
v++;
}
v=1,t=0;
repb(i,n,1){
t += (double)(a[i+1]-a[i])/v;
tn[i] = t;
v++;
}
int mid,l = 0,r = n+1;
while(l<r){//找第一个正向会后到的点 ts>tn的
mid = (l+r)>>1;
if(ts[mid]<=tn[mid]) l = mid+1;
else r = mid;
}
double dist = a[l]-a[l-1];
double res = max(ts[l-1],tn[l]);
int vs = l-1+1,vn = n+1-l+1;
if(ts[l-1]>tn[l]){
swap(vs,vn);
}
dist -= vs*fabs(ts[l-1]-tn[l]);//vs先到多走一段
res+=(double)dist/(vs+vn);
printf("%llf\n",res);
}
int main(){
int z;
cin>>z;
while(z--) solve();
}
这题绕了很多弯子,但是最后的收获是发现cout输出double的时候是输出六位有效数字,printf输出double的时候默认是六位小数,我之前就因为用了cout输出所以出现了精度问题。
双指针模拟的写法:
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<map>
#include<queue>
#include<cmath>
using namespace std;
template<class T>inline void read(T &x){x=0;char o,f=1;while(o=getchar(),o<48)if(o==45)f=-f;do x=(x<<3)+(x<<1)+(o^48);while(o=getchar(),o>47);x*=f;}
int cansel_sync=(ios::sync_with_stdio(0),cin.tie(0),0);
#define ll long long
#define ull unsigned long long
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define repb(i,a,b) for(int i=(a);i>=b;i--)
#define log(x) (31-__builtin_clz(x))
#define INF 0x3f3f3f3f
ll gcd(ll a,ll b){ while(b^=a^=b^=a%=b); return a; }
//#define INF 0x7fffffff
const int MAXN = 1e5+5;
int n,len;
int a[MAXN];
void solve(){
scanf("%d%d",&n,&len);
rep(i,1,n) scanf("%d",&a[i]);
double t1=0,t2=0;//记录正反两人的时间
int l = 0,r = n+1;//双指针
a[0] = 0;a[n+1] = len;
int v1=1,v2=1;
while(l+1<r){//直到两人在相邻的加速点上
if(t1+(double)(a[l+1]-a[l])/v1<t2+(double)(a[r]-a[r-1])/v2)
t1+=(double)(a[l+1]-a[l])/v1,v1++,l++;
else t2+=(double)(a[r]-a[r-1])/v2,v2++,r--;
}
if(t1>t2) swap(v1,v2);
printf("%llf\n", max(t1,t2)+(a[r]-a[l]-v1*fabs(t1-t2))/(v1+v2));
}
int main(){
int z;
scanf("%d",&z);
while(z--) solve();
}
压了一下行,少了很多不必要的变量,并且因为少掉了正反两次遍历,效率比之前更高
D.Searchlights
暴力,计算往右走了x步时,至少需要向上走多少步才能全部脱离探照灯范围。
通过枚举n和m中的点,维护函数f的值
f[b[j].x-a[i].x]=max(f[b[j].x-a[i].x],b[j].y-a[i].y+1);
但是这里的f(x)并不是“往右走了x步时,至少需要向上走多少步才能全部脱离探照灯范围”
这个问题的实际答案应该是
g
(
x
)
=
m
a
x
{
f
[
i
]
,
i
>
=
x
}
g(x)=max\left\{f[i],i>=x\right\}
g(x)=max{f[i],i>=x}
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<map>
#include<queue>
#include<cmath>
using namespace std;
template<class T>inline void read(T &x){x=0;char o,f=1;while(o=getchar(),o<48)if(o==45)f=-f;do x=(x<<3)+(x<<1)+(o^48);while(o=getchar(),o>47);x*=f;}
#define ll long long
#define ull unsigned long long
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define repb(i,a,b) for(int i=(a);i>=b;i--)
#define log(x) (31-__builtin_clz(x))
#define INF 0x3f3f3f3f
ll gcd(ll a,ll b){ while(b^=a^=b^=a%=b); return a; }
//#define INF 0x7fffffff
const int MAXN = 2020;
int n,m;
struct Node{
int x,y;
};
int f[1000010]={0};
Node a[MAXN],b[MAXN];
int main(){
cin>>n>>m;
rep(i,1,n) cin>>a[i].x>>a[i].y;
rep(i,1,m) cin>>b[i].x>>b[i].y;
rep(i,1,n)
rep(j,1,m)
if(b[j].x>=a[i].x)
f[b[j].x-a[i].x] = max(f[b[j].x-a[i].x],b[j].y-a[i].y+1);
//f[x]记录往右走了x步之后至少需要向上走多少步才能脱离,但这里的f[x]并步完全准确
//实际应该往上走的值应该是max{f[i],i>=x}
int ans = INF;
int max_back = 0;//记录当前i到末尾位置的f最大值
repb(i,1000009,0){
max_back = max(max_back,f[i]);
ans = min(ans,max_back+i);
}
cout<<ans<<endl;
}
之后应该多练习一下稍难一点的题,不能老是签个到就走人呜呜
E.Avoid Rainbow Cycles
这题难点在于想到怎么建图。
建立一个二分图,左边每个集合建一个点Li,右边给1到n的每个数字建一个点Rj。
如果集合i中包含数字j,则建立Li到Rj的一条边。
在这个图上跑最大生成树即可
可以观察到如果出现了彩虹圈圈,在二分图上也会出现环。比如A1连了123三个点,也会有12,13,23三条边,如果删去1这条边,就只剩下23一条边,就相当于删去了集合A1中的这个元素。
这题建图还挺奇妙的。
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<map>
#include<queue>
#include<cmath>
using namespace std;
template<class T>inline void read(T &x){x=0;char o,f=1;while(o=getchar(),o<48)if(o==45)f=-f;do x=(x<<3)+(x<<1)+(o^48);while(o=getchar(),o>47);x*=f;}
int cansel_sync=(ios::sync_with_stdio(0),cin.tie(0),0);
#define ll long long
#define ull unsigned long long
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
#define repb(i,a,b) for(int i=(a);i>=b;i--)
#define log(x) (31-__builtin_clz(x))
#define INF 0x3f3f3f3f
ll gcd(ll a,ll b){ while(b^=a^=b^=a%=b); return a; }
//#define INF 0x7fffffff
const int MAXN = 1e5+5;
ll a[MAXN],b[MAXN];
int n,m;
struct Edge{
int u,v; ll k;
friend bool operator < (Edge a,Edge b){//边权大的排在前面
return a.k>b.k;
}
}edges[MAXN<<1];//题目保证元素数少于2e5;
//把1到m编给set,m+1到m+n编给数字节点
//并查集部分
int fa[MAXN<<1];
int rankk[MAXN<<1];//并查集树高优化
int find(int x){
if(fa[x]==x) return x;
else return fa[x] = find(fa[x]);//路径压缩
}
void merge(int x,int y){//树高优化
x = find(x);y = find(y);
if(x==y) return ;
if(rankk[x]>rankk[y]) fa[y] = x;
else{
fa[x] = y;
if(rankk[x]==rankk[y]) rankk[y]++;
}
}
int main(){
cin>>m>>n;
rep(i,1,m) cin>>a[i];
rep(i,1,n) cin>>b[i];
int si,hc;
int cnt = 0;//记录边数
rep(i,1,m){
cin>>si;
rep(j,1,si){
cnt++;
edges[cnt].u=i;
cin>>hc;
edges[cnt].v=m+hc;//x对应节点m+x;
edges[cnt].k=a[i]+b[hc];
}
}
sort(edges+1,edges+1+cnt);//边权大的往前排
rep(i,1,m+n)//初始化并查集
fa[i] = i,rankk[i] = 0;
ll ans = 0;
rep(i,1,cnt){
if(find(edges[i].u)!=find(edges[i].v)){
//cout<<"check"<<edges[i].u<<' '<<edges[i].v<<' '<<edges[i].k<<endl;
//这些不划掉的边从总边权中删去,最终会得到划掉的边的花费总和
merge(edges[i].u,edges[i].v);
}
else ans += edges[i].k;//这个边划掉所支付的代价
}
cout<<ans<<endl;
}