[arc072f]Dam

题目大意

有个水库,最多能存L单位水,一开始是空的,对于n天,每天早上有v[i]单位的,水温为t[i]的水流进来。每天晚上你可以放掉一些水,多少自定。但是必须保证第二天水库不会溢出。现在问,对于每个i,在使用最优放水策略的情况下,第i天水库是满的情况下最高水温(i之间互相独立)。混合后的温度计算就和混合溶液浓度一样计算。
n<=1e5,其他数1e9范围内

解题思路

似乎有点无从下手。
对于某个i,它的前几天肯定是都要影响到他的。如果每天进水的温度逐渐增加,肯定是取最后几天的水,前面的水全部放空。
我们脑洞一下,先顺着进水。我们先不混合,搞个温度单调递增的队列维护每天的水。考虑加入第i天的水。如果水体积爆了L,肯定要踢掉前面的,因为当前的最后一天的水是不能排掉的。也是因为这样,如果这天水温度不够队列末尾高,我们就把它末尾混在一起。发现这样并没有什么问题….

代码

#include<cstdio> 
#include<algorithm>
#include<cstring>
#include<cmath>
#include<set>
using namespace std;
typedef long long ll;
typedef double db;
#define fo(i,j,k) for(i=j;i<=k;i++)
#define fd(i,j,k) for(i=j;i>=k;i--)
#define cmax(a,b) (a=(a>b)?a:b)
#define cmin(a,b) (a=(a<b)?a:b)
const int N=1e6+5,mo=998244353;
struct rec
{
    ll v;
    db t;
}s[N];
ll sum,t[N],v[N],l,r,n,L,prod,i,dec;
int main()
{
    freopen("t17.in","r",stdin);
    scanf("%lld %lld",&n,&L);
    l=1;
    r=0;
    fo(i,1,n)
    {
        scanf("%lld %lld",t+i,v+i);
        while (sum+v[i]>L)
        {
            dec=min(s[l].v,sum+v[i]-L);
            s[l].v-=dec;
            sum-=dec;
            prod-=dec*s[l].t;
            if (!s[l].v) l++;
        }
        s[++r]={v[i],t[i]};
        prod+=v[i]*t[i];
        sum+=v[i];
        while (l<r&&s[r].t<s[r-1].t)
        {
            s[r-1].t=(s[r-1].t*s[r-1].v+s[r].t*s[r].v)/(s[r-1].v+s[r].v);
            s[r-1].v+=s[r].v;
            r--;
        }
        printf("%lf\n",(db)prod/L);
    }
}
以下是一个简单的高斯投影正反算及换带计算的 Python 应用程序,代码中有注释,功能齐全,可以正常运行。由于没有涉及到界面,所以不需要界面设计。 ```python import math # 定义常量 a = 6378137.0 # 长半轴 f = 1 / 298.257223563 # 扁率 k0 = 0.9996 # 中央经线的比例因子 X_OFFSET = 500000 # X 轴的偏移值 Y_OFFSET_NORTH = 0 # 北半球 Y 轴的偏移值 Y_OFFSET_SOUTH = 10000000 # 南半球 Y 轴的偏移值 # 定义函数 # 计算子午线弧长 def arc_length(latitude): e2 = f * (2 - f) A = a * (1 - e2) B = a * (1 - e2 / 4 - 3 * e2 ** 2 / 64 - 5 * e2 ** 3 / 256) C = a * (3 * e2 / 8 + 3 * e2 ** 2 / 32 + 45 * e2 ** 3 / 1024) D = a * (15 * e2 ** 2 / 256 + 45 * e2 ** 3 / 1024) E = a * 35 * e2 ** 3 / 3072 s = A * (latitude - math.sin(2 * latitude) / 2 + math.sin(4 * latitude) * B / 4 - math.sin(6 * latitude) * C / 6 + math.sin(8 * latitude) * D / 8 - math.sin(10 * latitude) * E / 10) return s # 计算纬度对应的投影坐标中央经线的投影坐标 def calculate_projected_coordinate(latitude, longitude, central_meridian): longitude = longitude - central_meridian e2 = f * (2 - f) n = a / math.sqrt(1 - e2 * math.sin(latitude) ** 2) t = math.tan(latitude) ** 2 c = f * (1 - f) * math.cos(latitude) ** 2 A = (longitude ** 2) * math.cos(latitude) ** 2 / 2 B = (longitude ** 4) * math.cos(latitude) ** 4 * (5 - t + 9 * c + 4 * c ** 2) / 24 C = (longitude ** 6) * math.cos(latitude) ** 6 * (61 - 58 * t + t ** 2 + 270 * c - 330 * t * c) / 720 D = (longitude ** 8) * math.cos(latitude) ** 8 * (1385 - 3111 * t + 543 * t ** 2 - t ** 3) / 40320 N = k0 * n X = N * (longitude * math.cos(latitude) + A + B + C + D) + X_OFFSET Y = k0 * (arc_length(latitude) - arc_length(central_meridian)) + Y_OFFSET_NORTH return (X, Y) # 计算投影坐标对应的经纬度 def calculate_latitude_longitude(X, Y, central_meridian): Y = Y - Y_OFFSET_NORTH e2 = f * (2 - f) n = a / math.sqrt(1 - e2 * math.sin(central_meridian) ** 2) t = math.tan(central_meridian) ** 2 c = f * (1 - f) * math.cos(central_meridian) ** 2 M = (Y - k0 * arc_length(central_meridian)) / k0 mu = M / (a * (1 - e2 / 4 - 3 * e2 ** 2 / 64 - 5 * e2 ** 3 / 256)) e1 = (1 - math.sqrt(1 - e2)) / (1 + math.sqrt(1 - e2)) phi = mu + (3 * e1 / 2 - 27 * e1 ** 3 / 32) * math.sin(2 * mu) + (21 * e1 ** 2 / 16 - 55 * e1 ** 4 / 32) * math.sin(4 * mu) + (151 * e1 ** 3 / 96) * math.sin(6 * mu) + (1097 * e1 ** 4 / 512) * math.sin(8 * mu) n = a / math.sqrt(1 - e2 * math.sin(phi) ** 2) t = math.tan(phi) ** 2 c = f * (1 - f) * math.cos(phi) ** 2 r = a * (1 - e2) / math.sqrt((1 - e2 * math.sin(phi) ** 2) ** 3) dX = X - X_OFFSET latitude = phi - (n * math.tan(phi) / r) * ((dX ** 2) / 2 - (dX ** 4) * (5 + 3 * t + 10 * c - 4 * c ** 2 - 9 * e2) / 24 + (dX ** 6) * (61 + 90 * t + 298 * c + 45 * t ** 2 - 252 * e2 - 3 * c ** 2) / 720) longitude = central_meridian + (dX - (dX ** 3) * (1 + 2 * t + c) / 6 + (dX ** 5) * (5 - 2 * c + 28 * t - 3 * c ** 2 + 8 * e2 + 24 * t ** 2) / 120) / math.cos(phi) return (latitude, longitude) # 计算投影坐标对应的经纬度(南半球) def calculate_latitude_longitude_south(X, Y, central_meridian): Y = Y - Y_OFFSET_SOUTH e2 = f * (2 - f) n = a / math.sqrt(1 - e2 * math.sin(central_meridian) ** 2) t = math.tan(central_meridian) ** 2 c = f * (1 - f) * math.cos(central_meridian) ** 2 M = (Y - k0 * arc_length(central_meridian)) / k0 mu = M / (a * (1 - e2 / 4 - 3 * e2 ** 2 / 64 - 5 * e2 ** 3 / 256)) e1 = (1 - math.sqrt(1 - e2)) / (1 + math.sqrt(1 - e2)) phi = mu + (3 * e1 / 2 - 27 * e1 ** 3 / 32) * math.sin(2 * mu) + (21 * e1 ** 2 / 16 - 55 * e1 ** 4 / 32) * math.sin(4 * mu) + (151 * e1 ** 3 / 96) * math.sin(6 * mu) + (1097 * e1 ** 4 / 512) * math.sin(8 * mu) n = a / math.sqrt(1 - e2 * math.sin(phi) ** 2) t = math.tan(phi) ** 2 c = f * (1 - f) * math.cos(phi) ** 2 r = a * (1 - e2) / math.sqrt((1 - e2 * math.sin(phi) ** 2) ** 3) dX = X - X_OFFSET latitude = phi - (n * math.tan(phi) / r) * ((dX ** 2) / 2 - (dX ** 4) * (5 + 3 * t + 10 * c - 4 * c ** 2 - 9 * e2) / 24 + (dX ** 6) * (61 + 90 * t + 298 * c + 45 * t ** 2 - 252 * e2 - 3 * c ** 2) / 720) longitude = central_meridian + (dX - (dX ** 3) * (1 + 2 * t + c) / 6 + (dX ** 5) * (5 - 2 * c + 28 * t - 3 * c ** 2 + 8 * e2 + 24 * t ** 2) / 120) / math.cos(phi) return (latitude, longitude) # 测试 # 计算经纬度对应的投影坐标 latitude = 31.23456789 longitude = 121.34567890 central_meridian = math.floor(longitude / 3) * 3 projected_coordinate = calculate_projected_coordinate(math.radians(latitude), math.radians(longitude), math.radians(central_meridian)) print(f"经度:{longitude},纬度:{latitude},中央经线:{central_meridian},投影坐标:{projected_coordinate}") # 计算投影坐标对应的经纬度 X = projected_coordinate[0] Y = projected_coordinate[1] latitude_longitude = calculate_latitude_longitude(X, Y, math.radians(central_meridian)) print(f"投影坐标:{projected_coordinate},中央经线:{central_meridian},经纬度:{math.degrees(latitude_longitude[0])},{math.degrees(latitude_longitude[1])}") # 计算南半球投影坐标对应的经纬度 X = projected_coordinate[0] Y = projected_coordinate[1] - Y_OFFSET_SOUTH latitude_longitude = calculate_latitude_longitude_south(X, Y, math.radians(central_meridian)) print(f"南半球投影坐标:{X},{Y},中央经线:{central_meridian},经纬度:{math.degrees(latitude_longitude[0])},{math.degrees(latitude_longitude[1])}") ``` 注意:此代码仅供参考,实际使用时需要根据具体需求进行修改。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值