题意
给你 n n n根木头,第 i i i根长 a i a_i ai
锯出 2 2 2段长度为 x x x的木头和 x x x段长度为 y y y的木头
木头锯完后不能再拼接且必须满足 x ≥ 2 , y ≥ 2 x\ge2,y\ge2 x≥2,y≥2
求 x × y x\times y x×y的最大值
n ≤ 5 × 1 0 5 , 2 ≤ a i ≤ 5 × 1 0 5 n\le5\times10^5,2\le a_i\le5\times10^5 n≤5×105,2≤ai≤5×105
题解
因为要锯 x x x段长度为 y y y的木头, y y y影响更大,所以考虑对于某一个 y y y怎么要求出最大的 x x x
对于一个 y y y能锯的的木头总数 c n t Y = ∑ i = 1 n ⌊ a i y ⌋ cntY=\sum_{i=1}^n\lfloor\frac{a_i}y\rfloor cntY=∑i=1n⌊yai⌋
剩下两根长度为 x x x的木头可以分两种情况讨论:①来自于同一根木头;②来自于不同木头
①:来自于同一根木头
假设限定 2 x ∈ [ k y , k y + y ) 2x\in[ky,ky+y) 2x∈[ky,ky+y),那么 2 x 2x 2x的最大值就是 k y + max { a i m o d y ∣ a i ≥ k y } ky+\max\{a_i\mod y|a_i\ge ky\} ky+max{aimody∣ai≥ky}
k k k从大到小枚举并记录 a i m o d y a_i\mod y aimody的最大值即可
剩余木头数量即为 c n t Y − k cntY-k cntY−k
②:来自于不同木头
假设限定 x ∈ [ k y , k y + y ) x\in[ky,ky+y) x∈[ky,ky+y),记录 a i m o d y a_i\mod y aimody的最大值 m 1 m_1 m1,次大值 m 2 m_2 m2
如果 x x x取 k y + m 2 ky+m_2 ky+m2,那么剩余木头数量就是 c n t Y − 2 k cntY-2k cntY−2k
如果 x x x取 k y + m 1 ky+m_1 ky+m1,这样就要把一段更长的木头(如果存在的话)多锯开一小段,那么剩余木头数量就是 c n t Y − 2 k − 1 cntY-2k-1 cntY−2k−1
再处理一些细节即可(详见代码)
记 A = max a i A=\max a_i A=maxai,时间复杂度 O ( n + ∑ y = 2 A ⌊ A y ⌋ ) = O ( n + A ln A ) O(n+\sum_{y=2}^A\lfloor\frac Ay\rfloor)=O(n+A\ln A) O(n+∑y=2A⌊yA⌋)=O(n+AlnA)
#include<bits/stdc++.h>
#define fp(i,a,b) for(register int i=a,I=b+1;i<I;++i)
#define fd(i,a,b) for(register int i=a,I=b-1;i>I;--i)
#define file(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)
template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
template<class T>inline bool cmin(T&a,const T&b){return a>b?a=b,1:0;}
using namespace std;
char ss[1<<17],*A=ss,*B=ss;
inline char gc(){return A==B&&(B=(A=ss)+fread(ss,1,1<<17,stdin),A==B)?-1:*A++;}
template<class T>inline void sd(T&x){
char c;T y=1;while(c=gc(),(c<48||57<c)&&c!=-1)if(c==45)y=-1;x=c-48;
while(c=gc(),47<c&&c<58)x=x*10+c-48;x*=y;
}
const int N=5e5+5,M=N<<1,inf=~0u>>1;
typedef long long ll;
typedef int arr[N];
struct Pair{
int Mod,Len;
Pair(int x=-1,int y=-1){Mod=x,Len=y;}
bool operator<(Pair b)const{return Mod==b.Mod?Len<b.Len:Mod<b.Mod;}
};
struct Rest{
Pair x,y;
void Up(Pair z){y=max(y,min(x,z)),x=max(x,z);}
};
int n,Mx;arr pr,cnt,sum;ll ans,cntY;
inline void Upd(int x,int y){x<2?0:cmax(ans,(ll)x*y);}
int main(){
#ifndef ONLINE_JUDGE
file("s");
#endif
sd(n);
int x;
fp(i,1,n)sd(x),++cnt[x],cmax(Mx,x);
fp(i,1,Mx){
sum[i]=sum[i-1]+cnt[i];
pr[i]=cnt[i]?i:pr[i-1];
}
fp(y,2,Mx){
cntY=0;Rest res;
for(int i=y;i<=Mx;i+=y)
cntY+=(ll)(sum[min(i+y-1,Mx)]-sum[i-1])*(i/y);
for(int k=Mx/y,m=Mx+1;k>=0;m=k*y,--k){
int val=pr[m-1],More=(res.x.Mod>=0)+(res.y.Mod>=0),ky=k*y;
if(val>=ky){
res.Up(Pair(val%y,val));
if(cnt[val]==1)val=pr[val-1];
if(val>=ky)res.Up(Pair(val%y,val));
}
if(res.x.Mod>=0)
Upd(min((ll)(ky+res.x.Mod)>>1,cntY-k),y);
if(res.y.Mod>=0){
Upd(min((ll)(ky+res.y.Mod),cntY-2*k),y);
if(More+(res.x.Len<m)>1)
Upd(min((ll)(ky+res.x.Mod),cntY-2*k-1),y);
}
}
}
printf("%lld\n",ans);
return 0;
}