C++层次分析法一致性检验

// ConsoleApplication3.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h"
#include<iostream>
#include<math.h>
#include<vector>
using namespace std;

int main()
{
    cout << "请输入行或列的长度" << endl;
    int row;
    cin >> row;
    cout << "请输入数字" << endl;
    vector<double> a;
    for (size_t i = 0; i < row*row; i++)
    {
        double x;
        cin >> x;
        a.push_back(x);
    }
    vector<double> w;
    double fenmu = 0, wzhi = 0;
    for (int i = 0; i < row; i++)
    {
        for (size_t j = 0; j < row; j++)
        {
            for (size_t k = 0; k < row; k++)
            {

                fenmu += a[k*row+j];
            }
            wzhi += a[i*row+j] / fenmu;
            fenmu = 0;
        }
        w.push_back(wzhi / row);
        wzhi = 0;
    }
    //以上为计算W值、
    vector<double> lanmude;
    double lanmudazhi = 0;
    for (size_t i = 0; i < row; i++)
    {
        for (size_t j = 0; j < row; j++)
        {
            lanmudazhi += a[i*row+j] * w[j];
        }
        lanmude.push_back(lanmudazhi);
        lanmudazhi = 0;
    }
    //以上是求lanmuda值
    vector<double> lmax;
    for (size_t i = 0; i < row; i++)
    {
        lmax.push_back(lanmude[i] / w[i]);
    }
    double max = 0;
    for (size_t i = 0; i < row; i++)
    {
        max += lmax[i];
    }
    max /= row;
    //以上是求max值
    double  ci = (max - row) / (row - 1);
    double  ri[5] = { 0,0,0.58,0.9,1.12 };
    double cr = ci / ri[row-1];
    //以上是求cr值
    cout << cr << endl;
    if (ci==0)
    {
        cout << "符合一致性检验";
    }
    else
    {
        if (cr<0.1)
        {
            cout << "符合一致性" << endl;;
        }
        else
        {
            cout << "不符合一致性应修改" << endl;

        }
    }
    system("pause");
    return 0;
}