P6160 [Cnoi2020] 三分 / 向量

题意

传送门 P6160 [Cnoi2020]向量

题解
三分

⟨ a ⃗ , b ⃗ ⟩ = α , ⟨ b ⃗ , c ⃗ ⟩ = β \langle\vec{a},\vec{b}\rangle=\alpha,\langle\vec{b},\vec{c}\rangle=\beta a ,b =α,b ,c =β,那么有 ⟨ a ⃗ , c ⃗ ⟩ = 2 π − ( α + β ) 或 ∣ α − β ∣ \langle\vec{a},\vec{c}\rangle=2\pi-(\alpha+\beta)或\lvert \alpha-\beta\rvert a ,c =2π(α+β)αβ 猜测凸函数性质,三分套三分即可。由于 l o n g   d o u b l e long\ double long double 有效数位与测试数据极大值接近,此方法卡精度较为困难。

#include <bits/stdc++.h>
using namespace std;
#define rep(i, l, r) for (int i = l, _ = r; i < _; ++i)
#define db long double
const db pi = acos((db)-1.0), eps = 2e-19;
db R[3];

db F(db a, db b)
{
    db res = R[0] * R[1] * cos(a);
    res += R[1] * R[2] * cos(b);
    res += R[2] * R[0] * min(cos(a + b), cos(b - a));
    return res;
}

db calc(db a)
{
    db lb = 0, ub = pi;
    rep(i, 0, 100)
    {
        db lmid = lb + (ub - lb) / 3, umid = ub - (ub - lb) / 3;
        if (F(a, lmid) + eps < F(a, umid))
            ub = umid;
        else
            lb = lmid;
    }
    return F(a, ub);
}

int main()
{
    rep(i, 0, 3) scanf("%Lf", R + i);
    db lb = 0, ub = pi;
    rep(i, 0, 100)
    {
        db lmid = lb + (ub - lb) / 3, umid = ub - (ub - lb) / 3;
        if (calc(lmid) + eps < calc(umid))
            ub = umid;
        else
            lb = lmid;
    }
    db res = calc(ub);
    printf("%.1Lf\n", res + res * eps);
    return 0;
}
向量

已知 ∣ a ⃗ + b ⃗ + c ⃗ ∣ 2 = ∣ a ⃗ ∣ 2 + ∣ b ⃗ ∣ 2 + ∣ c ⃗ ∣ 2 + 2 a ⃗ ⋅ b ⃗ + 2 b ⃗ ⋅ c ⃗ + 2 c ⃗ ⋅ a ⃗ \lvert\vec{a}+\vec{b}+\vec{c}\rvert^2=\lvert\vec{a}\rvert^2+\lvert\vec{b}\rvert^2+\lvert\vec{c}\rvert^2+2\vec{a}\cdot\vec{b}+2\vec{b}\cdot\vec{c}+2\vec{c}\cdot\vec{a} a +b +c 2=a 2+b 2+c 2+2a b +2b c +2c a 关注 ∣ a ⃗ + b ⃗ + c ⃗ ∣ 2 \lvert\vec{a}+\vec{b}+\vec{c}\rvert^2 a +b +c 2,若 ∣ a ⃗ ∣ + ∣ b ⃗ ∣ > ∣ c ⃗ ∣ \lvert\vec{a}\rvert+\lvert\vec{b}\rvert> \lvert\vec{c}\rvert a +b >c ,则可以构造三角形,此时向量和模值为 0 0 0;反之,取 a ⃗ , b ⃗ \vec{a},\vec{b} a ,b c ⃗ \vec{c} c 方向相反即可,此时向量和模值最小,为 ∣ c ⃗ ∣ − ∣ a ⃗ ∣ − ∣ b ⃗ ∣ \lvert\vec{c}\rvert-\lvert\vec{a}\rvert-\lvert\vec{b}\rvert c a b

#include <bits/stdc++.h>
using namespace std;
typedef long double db;
#define rep(i, l, r) for (int i = l, _ = r; i < _; ++i)
db R[3];

int main()
{
    rep(i, 0, 3) cin >> R[i];
    db res = 0;
    rep(i, 0, 3) res -= R[i] * R[i];
    db d = R[0] + R[1] > R[2] ? 0 : R[2] - R[1] - R[0];
    cout << fixed << setprecision(1) << (res + d * d) / 2 << '\n';
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值