[NOIP普及组]税收与补贴问题

描述

每样商品的价格越低，其销量就会相应增大。现已知某种商品的成本及其在若干价位上的销量（产品不会低于成本销售），并假设相邻价位间销量的变化是线性的且在价格高于给定的最高价位后，销量以某固定数值递减。（我们假设价格及销售量都是整数）

对于某些特殊商品，不可能完全由市场去调节其价格。这时候就需要政府以税收或补贴的方式来控制。（所谓税收或补贴就是对于每个产品收取或给予生产厂家固定金额的货币）

你是某家咨询公司的项目经理，现在你已经知道政府对某种商品的预期价格，以及在各种价位上的销售情况。要求你确定政府对此商品是应收税还是补贴的最少金额（也为整数），才能使商家在这样一种政府预期的价格上，获取相对其他价位上的最大总利润。

总利润=单位商品利润*销量

单位商品利润 = 单位商品价格 – 单位商品成本 （– 税金 or  + 补贴）

输入

每个测试文件只包含一组测试数据，每组输入的第一行为政府对某种商品的预期价，第二行有两个整数，第一个整数为商品成本，第二个整数为以成本价销售时的销量售，以下若干行每行都有两个整数，第一个为某价位时的单价，第二个为此时的销量，以一行-1，-1表示所有已知价位及对应的销量输入完毕，输入的最后一行为一个单独的整数表示在已知的最高单价外每升高一块钱将减少的销量。

输出

输出有两种情况：若在政府预期价上能得到最大总利润，则输出一个单独的整数，数的正负表示是补贴还是收税，数的大小表示补贴或收税的金额最小值。若有多解，取绝对值最小的输出。如在政府预期价上不能得到最大总利润，则输出"NO SOLUTION"。

输入样例 1 

31
28 130
30 120
31 110
-1 –1
15

输出样例 1

4

来源

sources: NOIP全国联赛普及组 2000年NOIP全国联赛普及组

#include<cstdio>
#include<algorithm>
#define maxn 100000
using namespace std;
struct tnode{int p,s;}d[maxn];
int p1,q=0,r0=0,r1,r2;
int p2,sum=1,ans=0;
void redirect()
{
  freopen("1.in","r",stdin);
  freopen("1.out","w",stdout);
}
bool cmp1(tnode a,tnode b)
{
  return a.p<b.p;
}
bool cmp2(tnode a,tnode b)
{
  if(a.s<b.s)return 1;
  if(a.s>b.s)return 0;
  return a.p>b.p;
}
void readdata()
{
  int i,j,k,l;
  scanf("%d%d%d",&p1,&d[sum].p,&d[sum].s);
  p2=d[sum].p,p1-=p2;
  while(d[sum].p!=-1 && d[sum].s!=-1)
    {
      d[sum].p-=p2,sum++;
      scanf("%d%d",&d[sum].p,&d[sum].s);
    }
  sum--,l=sum;
  sort(d+1,d+sum+1,cmp1);
  
  for(i=1;i<sum;i++)
    {
      if(d[i+1].p==d[i].p+1)continue;
      k=(d[i+1].s-d[i].s)/(d[i+1].p-d[i].p);
      d[++l].p=d[i].p+1,d[l].s=d[i].s+k;
      for(j=2;j<d[i+1].p-d[i].p;j++)
        d[++l].p=d[l-1].p+1,d[l].s=d[l-1].s+k;
    }
  
  scanf("%d",&k);
  d[++l].p=d[sum].p+1,d[l].s=d[sum].s-k;
  while(d[l].s>=k)d[++l].p=d[l-1].p+1,d[l].s=d[l-1].s-k;  
  
  sort(d+1,d+l+1,cmp2);
  for(sum=1,i=2;i<=l;i++)if(d[i].s!=d[i-1].s)d[++sum]=d[i];  
}
void work()
{
  int i,j,k,x;
  for(i=1;i<=sum;i++)if(d[i].p==p1){q=i;break;}
  if(q==0){printf("NO SOLUTION\n");return;}
  
  k=d[q].p*d[q].s;
  for(i=1;i<=sum;i++)if(d[i].p*d[i].s>k)r0++;
  if(r0==0){printf("0\n");return;}
  
  r1=r0;
  for(x=1;x<=maxn;x++)
    {
      k=(d[q].p+x)*d[q].s,r2=0;
      for(i=1;i<=sum;i++)if((d[i].p+x)*d[i].s>k)r2++;
      if(r2==0){ans=x;break;}
      if(r2>r1)break;
      r1=r2;
    }
  r1=r0;
  for(x=1;x<=d[q].p;x++)
    {
      k=(d[q].p-x)*d[q].s,r2=0;
      for(i=1;i<=sum;i++)if((d[i].p-x)*d[i].s>k)r2++;
      if(r2==0)
        {
          if(ans==0 || ans>x)ans=-x;
          break;        
        }
      if(r2>r1)break;
      r1=r2;
    }
  if(ans==0)printf("NO SOLUTION\n");
  else printf("%d\n",ans);  
}
int main()
{
  //redirect();
  readdata();
  work();
  return 0;
}