标签：函数 二次 ll 组装 long rg 零件 CQOI2012 贪心
solution:
蒟蒻表示并不会模拟退火，所以用了差分数组加贪心吗。我们先来看题：
- 在数轴上的某个位置修建一个组装车间
- 到组装车间距离的平方的最小值。
- 1<=n<=20000
心路历程：
- 在一条直线上
- 距离的平方？（二次函数？）
- 1<=n<=10000？（nlogn(logn.....)?）
嗯？乍一看还真不知道怎么做啊！
可是在x轴上？还要距离的平方，好像几个二次函数加一起还是二次函数来着？
莫非？这整体就是一个单峰函数？
好吧，零件种类有多种，但（既然是平方）应该还是要从二次函数出手吧。
既然每种物资只能选一个的话那干脆写一个分段函数吧。|xi|<=200000，他给的又是整形，分界点最多只有400000个，好像可行
所以：
贪心：它的范围是-100000～100000，但如果我们将一个零件看成一个二次函数，数轴看成x轴，因为他们都是平方，每个零件所构成的函数都是\(y={(x+k)}^2\)的形式（ｋ表示零件在数轴上的位置）。
所以对于每一类零件，建立一个坐标轴，将零件画成二次函数的形式，这样对于每一个ｘ值，我们贪心的选ｙ值最小的那个零件。这样就能构成一个分段函数，而我们将每一种类的零件所构成的分段函数放在一起（用差分数组维护\(y=ax^2+bx+c\) 中的b与c），几个二次函数加一起还是二次函数，虽然分的段特别多，但绝对在我们的承受范围内，然后我们用二次函数求最小值的公式把最小值求出来即可（最小值可以在段外！！！）
零件位置给的是整形，分段的点的小数部分绝对为0或0.5,我们在差分数组前先处理一下（乘个二，加上200001，避免负数）（当然！你也可以用莫对思想维护！（比如鸽王发的那篇题解））
code：
#include<iostream>
#include<cstdio>
#include<algorithm>
#define ll long long
#define db double
#define rg register int
using namespace std;
struct su{
int x,y;
}t[200005];
int n,m,l,r;
ll s[800015];
ll v,k[800015];//请不要吝啬您的long long和double
db a,b,c,res=(ll)1<<50,ans,tot;//res初值赋大些
inline bool cmp(su x,su y){
return x.y==y.y?x.x<y.x:x.y<y.y;
}
int main(){
//freopen("battle.in","r",stdin);
//freopen("battle.out","w",stdout);
cin>>n>>m;a=n;
for(rg i=1;i<=m;++i)
cin>>t[i].x>>t[i].y;
sort(t+1,t+m+1,cmp);
for(rg i=1;i<=m;++i){
if(t[i].y!=t[i-1].y)l=1;
else l=t[i].x+t[i-1].x+400001;
if(t[i].y!=t[i+1].y)r=800002;
else r=t[i].x+t[i+1].x+400001;
v=-2*t[i].x;s[l]+=v;s[r]-=v;
v>>=1;v=v*v;k[l]+=v;k[r]-=v;
}
for(rg i=1;i<=800001;++i){
b=b+s[i];c=c+k[i];tot=c-(b*b)/(4*a);
if(tot<res)res=tot,ans=-b/(2*a);
}printf("%.4lf",ans);
return 0;
}
标签：函数,二次,ll,组装,long,rg,零件,CQOI2012,贪心 来源： https://www.cnblogs.com/812-xiao-wen/p/10339355.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。