제출 #93486

#	제출 시각	아이디	문제	언어	결과	실행 시간	메모리
93486		jufreire	수열 (APIO14_sequence)	C++11	100 / 100	1292 ms	85744 KiB

sequence

#include <stdio.h>
#include <stdlib.h>

long long xf;

class solution
{
public:
  int lastcut; //"this" cut //**
  //int ncuts;
  //int cuts[200];
  //int prevcut;
  long long b; //valor anterior
  long long a; //soma do inicio ate o lastcut
               //  long long preva; //soma do zero ate o lastcut anterior
               //  long long melhor; //a partir de onde esta eh a melhor
  long long valueat(long long x) //valor no final se este fosse o ultimo corte
  {
    if (x < a) return -1;
    return a * (x - a) + b;
  }
  long long acima(solution nova)
  {
    long long num;
    long long den;
    num = a*a - nova.a*nova.a + nova.b - b;
    den = a - nova.a;

    if (den == 0) //retas paralelas
    {
      if (valueat(xf) >= nova.valueat(xf)) return xf+1;
      else return nova.a;
    }

    long long x = num/den;
    if (x<nova.a) return nova.a;
    if (valueat(x) < nova.valueat(x) ) return x;
    return x+1; //talvez tivesse arredondado para baixo
  }

  long long acima_bin(solution nova) // retorna A PARTIR DE ONDE a nova esta acima desta. Pode ser xf se nunca.
  {
    long long x0 = a;
    if (nova.a > x0) x0 = nova.a; //nao adianta comparar antes de comecar
    long long x1 = xf;
    long long m;
    if (valueat(xf)>nova.valueat(xf)) return xf + 1;
    if (valueat(x0)<nova.valueat(x0)) return x0; //solucao nova ja surge acima da antiga. o que retornar?
    do
    {
      m = (x0 + x1) / 2;
      if (valueat(m) > nova.valueat(m)) x0 = m;
      else x1 = m;
    } while (x1>x0 + 1);
    return x1;
  } //eh, aqui pode melhorar bastante. sao duas retas, nao preciso dessa busca binaria toda para saber onde se encontram.
};

int k;

/*
int compar(const void* p1, const void* p2) //dessa vez vou querer em ordem da menos para mais inclinada
{
solution * m1 = (solution *)p1;
solution * m2 = (solution *)p2;
long long a = m1->a - m2->a;
if (a>0) return 1;
if (a<0) return -1;
return 0;
}
*/

int tamanho[201];
solution  **s;

typedef struct anterior
{
  int cut;
  int prox;
} anterior;

//anterior ** pc; //isto vai guardar apenas o prevcut do solution

int ** pc2;

void tam(int vv)
{
  tamanho[vv] += 1000;
  s[vv] = (solution*)realloc(s[vv], tamanho[vv] * sizeof(solution));
//  pc[vv] = (anterior*)realloc(pc[vv], tamanho[vv] * sizeof(anterior));
}
void copysol(int vv, int to, int from)
{
  s[vv][to].lastcut = s[vv][from].lastcut; //**
//  s[vv][to].ncuts = s[vv][from].ncuts;
  s[vv][to].a = s[vv][from].a;
  s[vv][to].b = s[vv][from].b;
//  s[vv][to].prevcut = s[vv][from].prevcut;
//  pc[vv][to].cut = pc[vv][from].cut;
//  pc[vv][to].prox = pc[vv][from].prox;
}


int main()
{
  int n, i;
  long long * a;
  long long * b;

//  freopen("60","r",stdin);

  scanf("%d", &n);
  scanf("%d", &k);
  a = (long long *)malloc(n * sizeof(long long));
  for (i = 0; i<n; i++) scanf("%lld", &a[i]);
  b = (long long *)malloc(n * sizeof(long long));
  b[n - 1] = a[n - 1];
  for (i = n - 2; i >= 0; i--) b[i] = b[i + 1] + a[i];
  //b[i] guarda a soma de i ate o fim
  xf = b[0];

  s = (solution**)malloc((k + 1) * sizeof(solution*));
//  pc = (anterior**)malloc((k + 1) * sizeof(anterior*));
  pc2 = (int**)malloc((k + 1) * sizeof(int*));

  for (i = 0; i <= k; i++)
  {
//    pc[i] = 0;
    pc2[i] = (int*) malloc(n*sizeof(int));
    s[i] = 0;
    tamanho[i] = 1000;
    //tam(i);
  }
  tam(0);

  int nsols[201];
  nsols[0] = 1;
  for (i = 1; i <= k; i++) nsols[i] = 0;

//  s[0][0].lastcut = 0;
//  s[0][0].ncuts = 0;
  s[0][0].a = 0;
  s[0][0].b = 0;
//  s[0][0].prevcut = -1;
//  pc[0][0].prox = -1;
//  pc[0][0].cut = 0;

  int ki, si, ss, ci;
  long long x, t1, t2;
  for (ki = 1; ki <= k; ki++) //fazer o proximo corte
  {
    tam (ki);
    if (ki>1) free(s[ki-2]);
//    printf("debug k=%d\n",ki);
    si = 0; //primeira solucao anterior
    ss = nsols[ki - 1]; //quantidade delas
    for (ci = ki; ci<n; ci++) //tem que ter espaco no inicio para pelo menos ki-1 cortes
    {
      x = b[0] - b[ci]; //ponto onde tenho que analisar as solucoes anteriores
      while (si<ss - 1 && s[ki - 1][si + 1].lastcut < ci && s[ki - 1][si].valueat(x) <= s[ki - 1][si + 1].valueat(x)) si++; //se tudo der certo bastava um if, mas por via das duvidas
      //while (si<ss - 1 && pc[ki - 1][si + 1].cut < ci && s[ki - 1][si].valueat(x) <= s[ki - 1][si + 1].valueat(x)) si++; //se tudo der certo bastava um if, mas por via das duvidas
      s[ki][nsols[ki]].lastcut = ci; //**
//      s[ki][nsols[ki]].ncuts = ki;
      s[ki][nsols[ki]].a = x;
      s[ki][nsols[ki]].b = s[ki - 1][si].valueat(x);
      //for (i=0;i<ki-1;i++) s[ki][nsols[ki]].cuts[i] = s[ki - 1][si].cuts[i];
      //s[ki][nsols[ki]].cuts[i] = ci;
//      s[ki][nsols[ki]].prevcut = si; //na verdade nao eh o cut, eh a solucao que tem o cut, pra ler recursivamente.
//      pc[ki][nsols[ki]].prox = si;
//      pc[ki][nsols[ki]].cut = ci;
      pc2[ki][ci] = s[ki-1][si].lastcut;

      //esta eh a nova solucao proposta. agora tenho que decidir se ela sera incluida ou nao.
//      printf("debug: %lld x + %lld - ", s[ki][nsols[ki]].a, s[ki][nsols[ki]].b);
      if (nsols[ki] == 0 || s[ki][nsols[ki] - 1].acima(s[ki][nsols[ki]]) <= xf) //ta errado
      {
        //        printf("accepted\n");
        nsols[ki]++;
        if (nsols[ki] == tamanho[ki]) tam(ki);
        //foi inserida na ordem certa de coeficiente angular, nem preciso ordenar.
        //so preciso ver se tem alguem antes dela que foi toda coberta. 
        //e elas nao estao espalhadas. se alguem foi coberto, foram as ultimas.
        i = nsols[ki] - 1;
        do
        {
          if (i<2) break;
          t1 = s[ki][i - 1].acima(s[ki][i]); //i passa de i-1
          t2 = s[ki][i - 2].acima(s[ki][i - 1]); //i-1 passa de i-2
          if (t1>t2) break; //porque esse eh o certo, cruzar em ordem.
                            //          printf("  removing old solution %lld %lld\n", s[ki][i - 1].a, s[ki][i - 1].b);
          copysol(ki, i - 1, i);
          i--;
          nsols[ki]--;
        } while (1);
      }
      else
      {
        //        printf("rejected\n");
      }

    }
//    printf("%d\n",nsols[ki]);
  }

  //a que ganha eh a mais inclinada com k cortes no valor de xf
  printf("%lld\n", s[k][nsols[k] - 1].valueat(xf));
  //for (i=0;i<k;i++) printf("%d ",s[k][nsols[k] - 1].cuts[i]);
  int pr[201];
  int pr2[201];
  int p = nsols[k] - 1;
  
  /*
  for (i = k; i >= 1; i--)
  {
    pr[i] = s[i][p].lastcut;
    p = s[i][p].prevcut;
  }
  p = nsols[k] - 1;
  for (i = k; i >= 1; i--)
  {
    pr2[i] = pc[i][p].cut;
    p = pc[i][p].prox;
  }

  for (i = 1; i <= k; i++) printf("%d ", pr2[i]);
  printf("\n");
  */


  //tentando o novo
  p = s[k][nsols[k] - 1].lastcut;

  for (i = k; i >= 1; i--)
  {
    pr2[i] = p;
    p = pc2[i][p];
  }

  for (i = 1; i <= k; i++) printf("%d ", pr2[i]);
  printf("\n");



}

컴파일 시 표준 에러 (stderr) 메시지

sequence.cpp: In function 'int main()':
sequence.cpp:208:7: warning: unused variable 'pr' [-Wunused-variable]
   int pr[201];
       ^~
sequence.cpp:112:8: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
   scanf("%d", &n);
   ~~~~~^~~~~~~~~~
sequence.cpp:113:8: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
   scanf("%d", &k);
   ~~~~~^~~~~~~~~~
sequence.cpp:115:30: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
   for (i = 0; i<n; i++) scanf("%lld", &a[i]);
                         ~~~~~^~~~~~~~~~~~~~~

• Subtask 1: 11 / 11
• Subtask 2: 11 / 11
• Subtask 3: 11 / 11
• Subtask 4: 17 / 17
• Subtask 5: 21 / 21
• Subtask 6: 29 / 29