Codevs 1507 酒厂选址

1507 酒厂选址

时间限制: 1 s

空间限制: 128000 KB

题目等级 : 黄金 Gold

传送门

题目描述 Description

Abstinence（戒酒）岛的居民们酷爱一种无酒精啤酒。以前这种啤酒都是从波兰进口，但今年居民们想建一个自己的啤酒厂。岛上所有的城市都坐落在海边，并且由一条沿海岸线的环岛高速路连接。酒厂的投资者收集了关于啤酒需求量的信息，即每天各城市消费的啤酒桶数。另外还知道相邻城市之间的距离。每桶啤酒每英里的运费是1元。日运费是将所需要的啤酒从酒厂运到所有城市所必需的运费之和。日运费的多少和酒厂的选址有关。投资者想找到一个合适的城市来修建酒厂，以使得日运费最小。

请设计一个程序：从文件bre.in 读入城市的数目、相邻两城市间的距离以及每个城市消费的啤酒桶数，计算最小的日运费，将结果写到输出文件bre.out中。

输入描述 Input Description

第一行是一个整数n（5 <= n <= 10000） ，表示城市的数目。 城市沿高速路编号，使得相邻的城市的编号也相邻（城市1和n也被认为是相邻）。 以下的n行，每行有两个非负整数。第I+1行的数 zi、di分别是城市I每日的啤酒消费量（桶）和从城市I沿高速路到下一个城市的距离（英里）。高速路的总长不会超过65535 英里。每座城市的日消费量不会超过255桶。

输出描述 Output Description

一个整数，表示所需的最小日运费（元）。

样例输入 Sample Input

6

1 2

2 3

1 2

5 2

1 10

2 3

样例输出 Sample Output

41

数据范围及提示 Data Size & Hint

分类标签 Tags

启发式搜索 模拟 搜索

/*
n^2做法.
一开始搞了一个前缀和和一个后缀和.
然后跑的时候傻眼了（可能有负的然后各种特判）
又看了看图
处理所有权值
反向跑的时候减去正向跑的权然后取小
一个前缀和就搞定了.
ans 可能会特别大（巨坑）.
*/
#include<iostream>
#include<cstdio>
#define MAXN 10001
#define LL long long
using namespace std;
LL tot=1e18;
int n,s[MAXN],t[MAXN],max1;
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d%d",&t[i],&s[i]);
        max1+=s[i];
        s[i]+=s[i-1];
    }
    for(int k=1;k<=n;k++)
    {
        LL ans=0;
        for(int i=1;i<=n;i++)
        {
            int x=k,y=i;
            if(x>y) swap(x,y);
            ans+=min(s[y-1]-s[x-1],max1-(s[y-1]-s[x-1]))*t

[i];
        }
        tot=min(tot,ans);
    }
    printf("%lld",tot);
    return 0;
}

/*
o(n)做法.
二分答案类似.
闲的我比着Pascal题解打了一遍.
枚举酒厂位置，
可以找到从一个点开始为分界，
一部分点跑左边，一部分跑右边，
该点记为mid，
mid会随着酒厂的右移而右移，具有二分性。
每次移动时计算答案贡献.
*/
#include<iostream>
#include<cstdio>
#define MAXN 10001
using namespace std;
long long tot,best,ans,ansm,mid,ans1mid,ans2mid;
long long l,r,min1,s[MAXN],t[MAXN],n,m,tots,tott;
int main()
{
    scanf("%d",&n);
    for(int i=0;i<n;i++)
    {
        scanf("%d%d",&t[i],&s[i]);
        tots+=s[i];tott+=t[i];
    }
    for(int i=0;i<n;i++)//确定第一次mid的位置.
    {
        m+=s[i];
        min1+=t[i];
        if(m>(tots>>1))
        {
            mid=i;
            m-=s[i];
            l=min1;
            r=tott-min1;
            break;
        }
    }
    for(int i=1;i<=mid;i++)//正向累加.
    {
        tot+=s[i-1];
        ans+=t[i]*tot;
    }
    tot=0;
    for(int i=n-1;i>=mid+1;i--)//反向累加.
    {
        tot+=s[i];
        ans+=t[i]*tot;
    }
    best=ans;//best最优解.
    for(int i=1;i<n;i++)
    {
        m-=s[i-1];
        l-=t[i-1];
        ans1mid=mid;
        while(m<=(tots>>1))//求新的中间位置.
        {
            m+=s[mid];
            mid=(mid+1)%n;
            l+=t[mid];
        }
        l-=t[mid];
        mid=(mid+n-1)%n;
        m-=s[mid];
        r=tott-l;
        ans2mid=mid;
        ansm=m;
        while(ans2mid!=ans1mid)
        {
            ans+=t[ans2mid]*(ansm+s[i-1]-(tots-ansm-s[i-

1]));
            ans2mid=(ans2mid+n-1)%n;
            ansm-=s[ans2mid];
        }
        ans+=(r-l)*s[i-1];
        best=min(best,ans);
        if(i==mid)//特殊情况.
        {
            m=tots;
            l=tott;
            ans=0;
            int tj,j=(i+1)%n;
            tj=s[i];
            while(j!=i)
            {
                ans+=t[j]*tj;
                tj+=s[j];
                j=(j+1)%n;
            }
        }
    }
    printf("%lld",best);
    return 0;
}