#!/usr/bin/perl -w

# mié 01 sep 2021 03:16:24 CEST  ->     Esta versión 5 no debe ser publicada hasta revisión de mínima verosimilitud
# lun 04 oct 2021 10:31:18 CEST  ->     Como no hacen comentarios, se publica

use strict;
# use Memoize;

# memoize("rollpob");

my $DEBUG = 1;

my $INGRESO_FISCAL_2020 = 164_630; # millones €  -> es el más bajo desde 2011

my @pob;
my @tmp;
open IN, '<', "pob_edad.csv" or die $!;
@tmp = <IN>;
for (@tmp){
    chomp;
    chop; # windows
    next if (/^#/);
    m/^(\d+) \s (.+) \; (\d+)$/x;
    if (defined $1 && defined $3) {
	$pob[$1] = $3;
	print "$1 $3\n" if $DEBUG;
    }else{
	warn "DATO CHUNGO";
    }
}
close IN;

my $SUPUESTO_FIN_TRABAJO = 67;    # pone rentas de mayores de 65 luego hay duda

my $SUPUESTA_COTIZACION = 0.279;  # by hand

my @renta = (0) x 101;
open IN2, '<', "renta_edad.csv" or die $!; # problema: es renta NETA media "estadística" y lo que interesa es INGRESO INTEGRO TRABAJO FISCAL
@tmp = <IN2>;
for (@tmp){
    chomp;
    chop;
    next if (/^#/);
    if (/^\+/ && /2020\;/){
	m/\;(\d+)$/;
	my $data = $1;                     # problema: son muy pocos datos
	warn "Dato chungo" unless defined $data; 
	$data /= (1 - $SUPUESTA_COTIZACION) ;   # así "aumento" la renta neta a la renta bruta
	print "$data\n" if $DEBUG;
	if (/enores/){
	    for (0..15){
		$renta[$_] = 0;   # $data; la renta neta media de menores de 16 está calculada pero no tiene interpretación porque no COTIZAN. La ESO es obligatoria, etc.
	    }
	}elsif (/16/ && /29/){
	    for (16 .. 29){
		$renta[$_] = $data;
	    }
	}elsif (/30/ && /44/){
	    for (30 .. 44){
		$renta[$_] = $data;
	    }
	}elsif (/45/ && /64/){
	    for (45 .. 64){
		$renta[$_] = $data;
	    }
	}elsif (/65/){
	    for (65 .. $SUPUESTO_FIN_TRABAJO){
		$renta[$_] = $data;
	    }
	}else{
	    warn "DATO CHUNGO";
	}
    }
}
close IN2;

# ya aquí la cosa se ve que está mal de datos, pero seguimos

my $ingresos_reales_2020 =  157_562_000_000;
my $gastos_reales_2020 = 171_883_286_720;

my $SUPUESTA_PENSION_MEDIA = 1566 *14;
$SUPUESTA_COTIZACION = 0.279;
my ($res1, $res2) = (1566*14, 0.279);
my ($res, $resant);
$resant = 9e99;
my $resant1 = $res1;
my $resant2 = $res2;
my $sumaingresos;
my $sumagastos;
# my ($flag1,$flag2) = (1,1);
my $kount;
my ($SPM, $SC) = ($SUPUESTA_PENSION_MEDIA, $SUPUESTA_COTIZACION);     # OJO A LO DE SIEMPRE
my ($SPMANT, $SCANT) = ($SPM, $SC);
do{
    $kount++;
    
    $SUPUESTA_PENSION_MEDIA =  $SPM + (rand()-0.5) / 10_000; # if $flag2;
    
    $SUPUESTA_COTIZACION =  $SC + (rand()-0.5) /10_000; #  if $flag1;    # como mucho, ya que desempleo, contingencias, etc. no cuenta 

    $sumaingresos = $sumagastos = 0;
    for my $i (0 .. scalar(@pob)-1){
	if ($i <= $SUPUESTO_FIN_TRABAJO){
	    $sumaingresos += $pob[$i] * $renta[$i] * $SUPUESTA_COTIZACION;
	}else{
	    $sumagastos += $pob[$i] * $SUPUESTA_PENSION_MEDIA;
	}
    }

    $res1 = abs($sumaingresos-$ingresos_reales_2020);
    $res2 = abs($sumagastos-$gastos_reales_2020);
    
    $res = $res1+$res2;
    print "Resultado = $res\n" if $DEBUG;
    
#    if (($res1+$res2) <= ($resant1+$resant2)){

    if ($resant1 >= $res1 || $res <= $resant){
	print "Supuesta cotización efectiva = $SUPUESTA_COTIZACION\n";
	$resant1 = $res1;
	$resant = $res;
	$SC = $SUPUESTA_COTIZACION;
	$SCANT = $SC;
    }else{
	$SC = $SCANT;
    }

    if ($resant2 >= $res2 || $res <= $resant){
	print "Supuesta pensión media = $SUPUESTA_PENSION_MEDIA\n";
	$resant2 = $res2;
	$resant = $res;
	$SPM = $SUPUESTA_PENSION_MEDIA;              # OJO CON LA CIFRA QUE VA AL PAN
	$SPMANT = $SPM;
    }else{
	$SPM = $SPMANT;
    }  

}until ($res < 100);

print "Error = $res Iteraciones = $kount\n";

print "Primera estimación de ingresos = $sumaingresos\n";
print "Primera estimación de gastos = $sumagastos\n";

print "\tEstimación primera déficit = ", $sumaingresos-$sumagastos, "\n";
print "\tParte del ingreso fiscal = ", ($sumaingresos-$sumagastos)/(1e6*$INGRESO_FISCAL_2020), "\n";

if (0) {
    my $comment = <<EOT;
    El Presupuesto consolidado de la Seguridad Social para 2021 presenta un volumen de gasto
      no financiero de 171.840,8 millones de euros, unos ingresos de 157.547,6 millones y, por
      tanto, un déficit de -14.293,2 millones de euros, en torno a un -1,3 % del PIB estimado
      para ese año
EOT
      ;
}

my $CRECIMIENTO_VEGETATIVO = -153_167;   # dato a 2021
my $POBLACION_TOTAL = 47_394_223;

my $PORCENTAJE_PROYECTAR = $CRECIMIENTO_VEGETATIVO / $POBLACION_TOTAL;


print "SUPUESTO 1: crecimiento vegetativo negativo con subida IPC\n";

print "El porcentaje a proyectar 2020 es de $PORCENTAJE_PROYECTAR\n";

# $PORCENTAJE_PR0YECTAR = -0.0012;   # según Climent, pero lo obvio y sigo con lo mio porque es mayor (más desigual) y año reciente


my $SUPUESTO_INCREMENTO = 0.015;  # IPC SIMULADO

my $sumaingresos2 = 0;
my $sumagastos2 = 0;
my @pob2 = rollpob(5);
for my $i (0 .. scalar(@pob)-1){
#    $pob2[$i] = (rollpob(5))[$i];
    if ($i <= $SUPUESTO_FIN_TRABAJO){
	$sumaingresos2 += $pob2[$i] * $renta[$i] * $SUPUESTA_COTIZACION * (1 + $PORCENTAJE_PROYECTAR) ** 5;        # ejemplo 5 años  v3. interés compuesto
    }else{
	$sumagastos2 += $pob2[$i] * $SUPUESTA_PENSION_MEDIA * (1 + $SUPUESTO_INCREMENTO) ** 5;
    }
}
print "Segunda estimación de ingresos = $sumaingresos2\n";
print "Segunda estimación de gastos = $sumagastos2\n";

print "\tEstimación segunda de déficit = ", $sumaingresos2-$sumagastos2, "\n";

my $INGRESO_FISCAL_2025 = $INGRESO_FISCAL_2020 * (1 + 0.06) * (1 + 0.03) ** 4; # el primer año crece al 6% y el resto al 3% v3. interés compuesto

print "\tParte del ingreso fiscal = ", ($sumaingresos2-$sumagastos2)/(1e6*$INGRESO_FISCAL_2025), "\n";



print "SUPUESTO 2: crecimiento vegetativo nulo con subida IPC\n";

my $sumaingresos3 = 0;
my $sumagastos3 = 0;
my @pob3 = rollpob(5);
for my $i (0 .. scalar(@pob)-1){
#    $pob3[$i] = (rollpob(5))[$i];
    if ($i <= $SUPUESTO_FIN_TRABAJO){
	$sumaingresos3 += $pob3[$i] * $renta[$i] * $SUPUESTA_COTIZACION * (1 + 0) ** 5;        # ejemplo 5 años v3. interés compuesto
    }else{
	$sumagastos3 += $pob3[$i] * $SUPUESTA_PENSION_MEDIA * (1 + $SUPUESTO_INCREMENTO) ** 5;
    }
}
print "3.ª estimación de ingresos = $sumaingresos3\n";
print "3.ª estimación de gastos = $sumagastos3\n";

print "\tEstimación 3.ª de déficit = ", $sumaingresos3-$sumagastos3, "\n";
print "\tParte del ingreso fiscal = ", ($sumaingresos3-$sumagastos3)/(1e6*$INGRESO_FISCAL_2025), "\n";



print "SUPUESTO 3: crecimiento vegetativo negativo con congelación pensiones\n";

my $sumaingresos4 = 0;
my $sumagastos4 = 0;
my @pob4 = rollpob(5);
for my $i (0 .. scalar(@pob)-1){
#    $pob4[$i] = (rollpob(5))[$i];
    if ($i <= $SUPUESTO_FIN_TRABAJO){
	$sumaingresos4 += $pob4[$i] * $renta[$i] * $SUPUESTA_COTIZACION * (1 + $PORCENTAJE_PROYECTAR) ** 5;        # ejemplo 5 años v3. interés compuesto
    }else{
	$sumagastos4 += $pob4[$i] * $SUPUESTA_PENSION_MEDIA * (1 + 0) ** 5;    # CONGELACION NOMINAL DE PENSIONES
    }
}
print "4.ª estimación de ingresos = $sumaingresos4\n";
print "4.ª estimación de gastos = $sumagastos4\n";

print "\tEstimación 4.ª de déficit = ", $sumaingresos4-$sumagastos4, "\n";
print "\tParte del ingreso fiscal = ", ($sumaingresos4-$sumagastos4)/(1e6*$INGRESO_FISCAL_2025), "\n";


print "SUPUESTO 4: crecimiento vegetativo negativo con semicongelación pensiones\n";

$SUPUESTO_INCREMENTO = 0.008;   # tiene que ser bajo 

print "Supuesto de incremento = $SUPUESTO_INCREMENTO\n";

my $sumaingresos5 = 0;
my $sumagastos5 = 0;
my @pob5 = rollpob(5);
for my $i (0 .. scalar(@pob)-1){
#    $pob5[$i] = (rollpob(5))[$i];
    if ($i <= $SUPUESTO_FIN_TRABAJO){
	$sumaingresos5 += $pob5[$i] * $renta[$i] * $SUPUESTA_COTIZACION * (1 + $PORCENTAJE_PROYECTAR) ** 5;        # ejemplo 5 años v3. interés compuesto
    }else{
	$sumagastos5 += $pob5[$i] * $SUPUESTA_PENSION_MEDIA * (1 + $SUPUESTO_INCREMENTO) ** 5;    # SEMICONGELACION NOMINAL DE PENSIONES
    }
}
print "5.ª estimación de ingresos = $sumaingresos5\n";
print "5.ª estimación de gastos = $sumagastos5\n";

print "\tEstimación 5.ª de déficit = ", $sumaingresos5-$sumagastos5, "\n";
print "\tParte del ingreso fiscal = ", ($sumaingresos5-$sumagastos5)/(1e6*$INGRESO_FISCAL_2025), "\n";

exit 1;

sub rollpob {
    my $n = shift;
    return @pob if $n <= 1;
    my @roll;
    for my $i (0 .. scalar(@pob)-1){
	$roll[$i+$n] = $pob[$i];
    }
    for my $i (0 .. scalar(@roll)-1){
	if ($i > 67){               # EDAD CHUNGA
	    $roll[$i] *= (1-0.09);              # TASA DE MORTALIDAD 9% (MENOS 2018 y 2020, MÁS QUE 2019) 
	}elsif ($i <= $n-1){
	    $roll[$i] = $pob[1];   # arbitrary truncation on last data known 
	}
    }
    return @roll;
}



__END__ 
  Esta estimación se puede hacer más precisa, lo cual estaría bien,
  pero descubre que un método efectivo para contener el desbalanceo 
  de ingresos y gastos es la contención del incremento de las pensiones 
  (la última simulación) y otra manera es seguir el marco del pacto
  de Toledo, financiando parte de las pensiones con impuestos otros.

Esta estimación en concreto 'est5' es tan mala que la única solución 
  satisfactoria es la congelación inversa de todo el resto de partidas
  presupuestarias para salvar las pensiones y que los incrementos 
  nominales del presupuesto se vayan trasladando a las pensiones.

A continuación incluyo un par de ficheros de datos que se usan:

-----------------------------------------------pob_edad.csv
# Sexo;Edad;Provincias;Periodo;Total
# Total;47.394.223
0 a�os;338786
1 a�o;363341
2 a�os;382059
3 a�os;408287
4 a�os;431472
5 a�os;443170
6 a�os;450016
7 a�os;445284
8 a�os;469239
9 a�os;486596
10 a�os;491380
11 a�os;501426
12 a�os;527716
13 a�os;510018
14 a�os;507740
15 a�os;498747
16 a�os;495260
17 a�os;488794
18 a�os;475565
19 a�os;479813
20 a�os;484847
21 a�os;479183
22 a�os;470774
23 a�os;479101
24 a�os;478719
25 a�os;481718
26 a�os;490899
27 a�os;509736
28 a�os;525781
29 a�os;523378
30 a�os;531524
31 a�os;539644
32 a�os;549056
33 a�os;555864
34 a�os;569230
35 a�os;588611
36 a�os;611266
37 a�os;631276
38 a�os;667423
39 a�os;697388
40 a�os;731529
41 a�os;753631
42 a�os;787828
43 a�os;798467
44 a�os;812938
45 a�os;806728
46 a�os;800884
47 a�os;781296
48 a�os;777575
49 a�os;765164
50 a�os;753516
51 a�os;745521
52 a�os;738061
53 a�os;741354
54 a�os;721854
55 a�os;715168
56 a�os;722272
57 a�os;686102
58 a�os;661337
59 a�os;644202
60 a�os;644806
61 a�os;624692
62 a�os;607662
63 a�os;590685
64 a�os;549081
65 a�os;530169
66 a�os;503637
67 a�os;500877
68 a�os;495402
69 a�os;460630
70 a�os;446989
71 a�os;459503
72 a�os;475082
73 a�os;429268
74 a�os;404434
75 a�os;416011
76 a�os;388845
77 a�os;374888
78 a�os;317167
79 a�os;284752
80 a�os;338743
81 a�os;217288
82 a�os;236104
83 a�os;252791
84 a�os;262547
85 a�os;237765
86 a�os;218866
87 a�os;204540
88 a�os;182691
89 a�os;153104
90 a�os;134152
91 a�os;107330
92 a�os;88281
93 a�os;66962
94 a�os;53675
95 a�os;40002
96 a�os;29854
97 a�os;21407
98 a�os;15093
99 a�os;9851
100 y m�s a�os;13043
# 85 y m�s a�os;1.576.614
--------------------------------------------------renta_edad.csv
# Sexo;Edad;Renta anual neta media por persona y por unidad de consumo;Periodo;Total
# Total;Renta  ;2020;12292
# Total;Renta  ;2019;11680
# Total;Renta  ;2018;11412
# Total;Renta  ;2017;11074
# Total;Renta  ;2016;10708
# Total;Renta  ;2015;10419
# Total;Renta  ;2014;10391
# Total;Renta  ;2013;10531
# Total;Renta  ;2012;10795
# Total;Renta  ;2011;10858
# Total;Renta  ;2010;11284
# Total;Renta  ;2009;11318
# Total;Renta  ;2008;10737
+Menores de 16 a�os;Renta  ;2020;9025
Menores de 16 a�os;Renta  ;2019;8899
Menores de 16 a�os;Renta  ;2018;8919
Menores de 16 a�os;Renta  ;2017;8354
Menores de 16 a�os;Renta  ;2016;8078
Menores de 16 a�os;Renta  ;2015;7841
Menores de 16 a�os;Renta  ;2014;7726
Menores de 16 a�os;Renta  ;2013;7741
Menores de 16 a�os;Renta  ;2012;8151
Menores de 16 a�os;Renta  ;2011;8259
Menores de 16 a�os;Renta  ;2010;8355
Menores de 16 a�os;Renta  ;2009;8444
Menores de 16 a�os;Renta  ;2008;8005
+De 16 a 29 a�os;Renta  ;2020;11031
De 16 a 29 a�os;Renta  ;2019;10345
De 16 a 29 a�os;Renta  ;2018;10156
De 16 a 29 a�os;Renta  ;2017;9815
De 16 a 29 a�os;Renta  ;2016;9373
De 16 a 29 a�os;Renta  ;2015;8935
De 16 a 29 a�os;Renta  ;2014;9004
De 16 a 29 a�os;Renta  ;2013;9438
De 16 a 29 a�os;Renta  ;2012;9727
De 16 a 29 a�os;Renta  ;2011;10032
De 16 a 29 a�os;Renta  ;2010;10651
De 16 a 29 a�os;Renta  ;2009;10797
De 16 a 29 a�os;Renta  ;2008;10321
+De 30 a 44 a�os;Renta  ;2020;12053
De 30 a 44 a�os;Renta  ;2019;11422
De 30 a 44 a�os;Renta  ;2018;11397
De 30 a 44 a�os;Renta  ;2017;11040
De 30 a 44 a�os;Renta  ;2016;10536
De 30 a 44 a�os;Renta  ;2015;10353
De 30 a 44 a�os;Renta  ;2014;10348
De 30 a 44 a�os;Renta  ;2013;10543
De 30 a 44 a�os;Renta  ;2012;10912
De 30 a 44 a�os;Renta  ;2011;11106
De 30 a 44 a�os;Renta  ;2010;11480
De 30 a 44 a�os;Renta  ;2009;11687
De 30 a 44 a�os;Renta  ;2008;11024
+De 45 a 64 a�os;Renta  ;2020;13626
De 45 a 64 a�os;Renta  ;2019;12988
De 45 a 64 a�os;Renta  ;2018;12550
De 45 a 64 a�os;Renta  ;2017;12311
De 45 a 64 a�os;Renta  ;2016;11931
De 45 a 64 a�os;Renta  ;2015;11591
De 45 a 64 a�os;Renta  ;2014;11626
De 45 a 64 a�os;Renta  ;2013;11921
De 45 a 64 a�os;Renta  ;2012;12270
De 45 a 64 a�os;Renta  ;2011;12374
De 45 a 64 a�os;Renta  ;2010;13047
De 45 a 64 a�os;Renta  ;2009;12906
De 45 a 64 a�os;Renta  ;2008;12413
+65 y m�s a�os;Renta  ;2020;14132
65 y m�s a�os;Renta  ;2019;13315
65 y m�s a�os;Renta  ;2018;12758
65 y m�s a�os;Renta  ;2017;12528
65 y m�s a�os;Renta  ;2016;12416
65 y m�s a�os;Renta  ;2015;12206
65 y m�s a�os;Renta  ;2014;12183
65 y m�s a�os;Renta  ;2013;11960
65 y m�s a�os;Renta  ;2012;11867
65 y m�s a�os;Renta  ;2011;11463
65 y m�s a�os;Renta  ;2010;11783
65 y m�s a�os;Renta  ;2009;11685
65 y m�s a�os;Renta  ;2008;10893
# Menos de 18 a�os;Renta  ;2020;9110
# Menos de 18 a�os;Renta  ;2019;8944
# Menos de 18 a�os;Renta  ;2018;8909
# Menos de 18 a�os;Renta  ;2017;8479
# Menos de 18 a�os;Renta  ;2016;8153
# Menos de 18 a�os;Renta  ;2015;7824
# Menos de 18 a�os;Renta  ;2014;7702
# Menos de 18 a�os;Renta  ;2013;7710
# Menos de 18 a�os;Renta  ;2012;8155
# Menos de 18 a�os;Renta  ;2011;8259
# Menos de 18 a�os;Renta  ;2010;8354
# Menos de 18 a�os;Renta  ;2009;8468
# Menos de 18 a�os;Renta  ;2008;8024
# De 18 a 64 a�os;Renta  ;2020;12624
# De 18 a 64 a�os;Renta  ;2019;11961
# De 18 a 64 a�os;Renta  ;2018;11727
# De 18 a 64 a�os;Renta  ;2017;11390
# De 18 a 64 a�os;Renta  ;2016;10948
# De 18 a 64 a�os;Renta  ;2015;10650
# De 18 a 64 a�os;Renta  ;2014;10660
# De 18 a 64 a�os;Renta  ;2013;10942
# De 18 a 64 a�os;Renta  ;2012;11252
# De 18 a 64 a�os;Renta  ;2011;11415
# De 18 a 64 a�os;Renta  ;2010;11964
# De 18 a 64 a�os;Renta  ;2009;12002
# De 18 a 64 a�os;Renta  ;2008;11428