算法分析与设计——动态规划图像压缩C++

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
#include <stdio.h>
#include <iostream>
#include <fstream>
#include <string.h>
using namespace std;
 
const int N = 7;
 
int length(int i);
void Compress(int n,int p[],int s[],int l[],int b[]);
void Tracebace(int n,int& i,int s[],int l[]);
void Output(int s[],int l[],int b[],int n);
 
int main()
{
    int p[] = {0,10,12,15,255,1,2};//图像数组{p1,p2,p3,p4,p5,p6} 下标从1开始计数
    int s[N],l[N],b[N];
 
    cout<<"图像的序列为："<<endl;
 
    for(int i=1;i<N;i++)
        cout<<p[i]<<" ";
    cout<<endl;
 
    Compress(N-1,p,s,l,b);
    Output(s,l,b,N-1);
    return 0;
}
 
void Compress(int n,int p[],int s[],int l[],int b[])
{//n个像素，存在数组p[]中，
    int Lmax = 256,header = 11;
    s[0] = 0;
    for(int i=1; i<=n; i++)
    {
        b[i] = length(p[i]);//计算像素点p需要的存储位数
        int bmax = b[i];
        s[i] = s[i-1] + bmax;
        l[i] = 1;
 
        for(int j=2; j<=i && j<=Lmax;j++)//i=1时，不进入此循环
        {
            if(bmax<b[i-j+1])
            {
                bmax = b[i-j+1];
            }
 
            if(s[i]>s[i-j]+j*bmax)
            {
                s[i] = s[i-j] + j*bmax;
                l[i] = j;
            }
        }
        s[i] += header;
    }
}
 
int length(int i)//存储像素pi所需的位数
{
    int k=1;
    i = i/2;
    while(i>0)
    {
        k++;
        i=i/2;
    }
    return k;
}
 
void Traceback(int n,int& i,int s[],int l[])
{
    if(n==0)
        return;
    Traceback(n-l[n],i,s,l);//p1,p2,p3,...,p(n-l[n])像素序列，最后一段有l[n]个像素
    s[i++]=n-l[n];//重新为s[]数组赋值，用来存储分段位置，最终i为共分了多少段
}
 
void Output(int s[],int l[],int b[],int n)
{
    //在输出s[n]存储位数后，s[]数组则被重新赋值，用来存储分段的位置
    cout<<"图像压缩后的最小空间为："<<s[n]<<endl;
    int m = 0;
    Traceback(n,m,s,l);//s[0]:第一段从哪分，s[1]:第二段从哪分...，s[m-1]第m段从哪分
    s[m] = n;//此时m为总段数，设s[m]=n，分割位置为第n个像素
    cout<<"将原序列分成"<<m<<"段序列段"<<endl;
    for(int j=1; j<=m; j++)
    {
        l[j] = l[s[j]];
        b[j] = b[s[j]];	//此时s数组存放值为分割的位置 
    }
    for(int j=1; j<=m; j++)
    {
        cout<<"段长度："<<l[j]<<",所需存储位数:"<<b[j]<<endl;
    }
}