herrDeng網內搜尋
自訂搜尋
熱門文章
-
11!=?
-
產生50個元素的陣列並排序
-
請輸出如下 * ** *** **** ***** ******
-
借錢500萬,年利18%，每月複利計算，請問連本帶利的金額 5, 10, 15, 20, 25, 30, 35, 40, 45, 50年的金額。並比較公式: A*exp(0.18n) A*(1+0.18/12)^(12*n), n=年數, A=本金 hint: math.h, p...
-
請問你所使的電腦，CPU的規格 Cache多大？其他硬體配置...... 另外， "INTEL® XEON PHI™ 處理器 7230F"的大致規格為？
-
寫一個是否上學期成績2/3的程式
-
計算你上學期加權平均
-
金鑰總數159*10^18 億=10^8 兆=10^12 京=10^16 垓=10^20 秭=10^24 穣=10^28 溝=10^32 澗=10^36 正=10^40 載=10^44 極=10^48 恆河沙=10^52 阿僧祇=10^56 那由...
-
用 C/C++算上學期加權平均
36 則留言:
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> X(70, 10);
double a[1000];
double sum = 0;
double sum2 = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
sum = sum / 1000.0;
cout << "平均數:" << sum << endl;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i]) - (sum * sum);
}
sum2 = sum2 / 1000;
double s = sqrt(sum2);
cout << "標準差:" << s << endl;
cout << (clock()-seed)/1000.0 <<"sec." << endl;
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> X(70, 10);
double a[1000];
double sum = 0;
double sum2 = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
sum = sum / 1000.0;
cout << "平均數:" << sum << endl;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i]) - (sum * sum);
}
sum2 = sum2 / 1000;
double s = sqrt(sum2);
cout << "標準差:" << s << endl;
cout << (clock()-seed)/1000.0 <<"sec." << endl;
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> X(70, 10);
double a[1000];
double sum = 0;
double sum2 = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
sum = sum / 1000.0;
cout << "平均數:" << sum << endl;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i]) - (sum * sum);
}
sum2 = sum2 / 1000;
double s = sqrt(sum2);
cout << "標準差:" << s << endl;
cout << (clock()-seed)/1000.0 <<"sec." << endl;
system("Pause");
return 0;
}
#include <iostream>
#include<cstdlib>
#include<random>
#include<ctime>
#include<cmath>
using namespace std;
int main(){
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> x(70, 10);
int a[1000];
double sum = 0;
double s = 0;
double sum2;
for (int i = 0; i < 1000; i++){
a[i] = (int)(round(x(gen)));
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9){
cout << endl;
}
}
sum2 = sum / 1000.0;
for (int i = 0; i < 1000; i++){
// s += (a[i] - sum2)*(a[i] - sum2);
s += (a[i] * a[i]) - (sum2*sum2);
}
s /= 1000;
s = sqrt(s);
cout << sum2;
cout << endl;
cout << s;
cout << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return (0);
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include<iostream>
#include<cstdlib>
#include<random>
#include<ctime>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> x(70, 10);
int a[1000];
double s = 0;
for (int i = 0; i < 1000; i++){
a[i] = x(gen);
cout << a[i] << "\t";
s = a[i]+s;
if (i % 10 == 9)cout << endl;
}
cout << (clock() - seed) / 10000.0 << "sec.\n";
double av = s / 1000.0;
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] - av)*(a[i] - av);
}
double var = sum2 / 1000.0;
double c = sqrt(var);
cout << c << "標準差";
cout << endl;
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> dis(70, 10);
double a[1000];
double sum = 0;
double sum2 = 0;
for (int i = 0; i < 1000; i++){
a[i] = dis(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum/1000.0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i]- av*av);
}
double var = sum2 / 1000.0;
cout << sum / 1000.0 << endl;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
cout << endl;
system("Pause");
return 0;
}
#include
#include
#include
#include
#include
using namespace std;
int main(){
int seed = clock();
mt19937_64 gen(seed);
normal_distribution x(70, 10);
int a[1000];
double sum = 0;
double s = 0;
double sum2;
for (int i = 0; i < 1000; i++){
a[i] = (int)(round(x(gen)));
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9){
cout << endl;
}
}
sum2 = sum / 1000.0;
for (int i = 0; i < 1000; i++){
s += (a[i] * a[i]) - (sum2*sum2);
}
s /= 1000;
s = sqrt(s);
cout << sum2;
cout << endl;
cout << s;
cout << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return (0);
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double>X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 100; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include <iostream>
#include<cstdlib>
#include<random>
#include<ctime>
#include<cmath>
using namespace std;
int main(){
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> x(70, 10);
int a[1000];
double sum = 0;
double s = 0;
double sum2;
for (int i = 0; i < 1000; i++){
a[i] = (int)(round(x(gen)));
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9){
cout << endl;
}
}
sum2 = sum / 1000.0;
for (int i = 0; i < 1000; i++){
s += (a[i] * a[i]) - (sum2*sum2);
}
s /= 1000;
s = sqrt(s);
cout << sum2;
cout << endl;
cout << s;
cout << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return (0);
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include <iostream> //程式如果用到cout或cin就必須引用iostream
#include <cstdlib>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;
int main(void) {
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> dis(70, 10);
int a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = (int)(round(dis(gen)));
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << "\n";
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double>X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 100; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include
#include
#include
#include
#include
using namespace std;
int main(){
int seed = clock();
mt19937_64 gen(seed);
normal_distributionX(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2+=(a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> dis(70, 10);
double a[100];
double sum = 0;
for (int i = 0; i < 100; i++){
a[i] = (int)(round(dis(gen)));
cout << a[i] << "\t";
sum+= a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 100.0;
cout <<"average= "<< av << "\n";
double sum2 = 0;
for (int i = 0; i < 100; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 100.0;
double s = sqrt(var);
cout <<"standard deviation= "<< s << endl;
cout << (clock()-seed/100.0)<<"sec. \n";
system("pause");
return 0;
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
B10233139 陳鈺璇
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
b10233137
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
B10233169
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
B10233139
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
B10233169
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9) cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include <iostream>
#include <cstdlib>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution<double> X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
double sum2 = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
sum = sum / 1000.0;
cout << "平均數:" << sum << endl;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i]) - (sum * sum);
}
sum2 = sum2 / 1000;
double s = sqrt(sum2);
cout << "標準差:" << s << endl;
cout << (clock()-seed)/1000.0 <<"sec." << endl;
system("Pause");
return 0;
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
double sum2 = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
sum = sum / 1000.0;
cout << "平均數:" << sum << endl;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i]) - (sum * sum);
}
sum2 = sum2 / 1000;
double s = sqrt(sum2);
cout << "標準差:" << s << endl;
cout << (clock()-seed)/1000.0 <<"sec." << endl;
system("Pause");
return 0;
}
#include
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
double av = sum / 1000.0;
cout << av << "\n";
double sum2 = 0;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i] - av*av);
}
double var = sum2 / 1000.0;
double s = sqrt(var);
cout << s << endl;
cout << (clock() - seed) / 1000.0 << "sec.\n";
system("Pause");
return 0;
}
ImageMatrix convert2green(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x + 0] = mm.term[y*m.rowsize + x + 2] = 0;
mm.term[y*m.rowsize + x + 1] = G;
}
return mm;
}
ImageMatrix convert2red(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x + 0] = mm.term[y*m.rowsize + x + 1] = 0;
mm.term[y*m.rowsize + x +2] = G;
}
return mm;
}
ImageMatrix convert2blue(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x +1] = mm.term[y*m.rowsize + x + 2] = 0;
mm.term[y*m.rowsize + x + 0] =G;
}
return mm;
}
ImageMatrix convert2half(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y<m.height; y++)
for (int x = 0; x<m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
BYTE gray = (BYTE)(0.2989*R + 0.5870*G + 0.1140*B + 0.5);
if (y<m.height/2)
mm.term[y*m.rowsize + x] = mm.term[y*m.rowsize + x + 1] =
mm.term[y*m.rowsize + x + 2] = gray;
else{
mm.term[y*m.rowsize + x + 0] = B;
mm.term[y*m.rowsize + x + 1] = G;
mm.term[y*m.rowsize + x + 2] = R;
}
}
return mm;
}
ImageMatrix convert2green(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x + 0] = mm.term[y*m.rowsize + x + 2] = 0;
mm.term[y*m.rowsize + x + 1] = G;
}
return mm;
}
ImageMatrix convert2red(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x + 0] = mm.term[y*m.rowsize + x + 1] = 0;
mm.term[y*m.rowsize + x + 2] = G;
}
return mm;
}
ImageMatrix convert2blue(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x + 1] = mm.term[y*m.rowsize + x + 2] = 0;
mm.term[y*m.rowsize + x + 0] = G;
}
return mm;
}
ImageMatrix convert2half(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y<m.height; y++)
for (int x = 0; x<m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
BYTE gray = (BYTE)(0.2989*R + 0.5870*G + 0.1140*B + 0.5);
if (y<m.height / 2)
mm.term[y*m.rowsize + x] = mm.term[y*m.rowsize + x + 1] =
mm.term[y*m.rowsize + x + 2] = gray;
else{
mm.term[y*m.rowsize + x + 0] = B;
mm.term[y*m.rowsize + x + 1] = G;
mm.term[y*m.rowsize + x + 2] = R;
}
}
return mm;
}
#include
#include
#include
#include
using namespace std;
int main()
{
int seed = clock();
mt19937_64 gen(seed);
normal_distribution X(70, 10);
double a[1000];
double sum = 0;
double sum2 = 0;
for (int i = 0; i < 1000; i++){
a[i] = X(gen);
cout << a[i] << "\t";
sum += a[i];
if (i % 10 == 9)cout << endl;
}
sum = sum / 1000.0;
cout << "平均數:" << sum << endl;
for (int i = 0; i < 1000; i++){
sum2 += (a[i] * a[i]) - (sum * sum);
}
sum2 = sum2 / 1000;
double s = sqrt(sum2);
cout << "標準差:" << s << endl;
cout << (clock()-seed)/1000.0 <<"sec." << endl;
system("Pause");
return 0;
}
ImageMatrix convert2green(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x + 0] = mm.term[y*m.rowsize + x + 2] = 0;
mm.term[y*m.rowsize + x + 1] = G;
}
return mm;
}
ImageMatrix convert2red(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x + 0] = mm.term[y*m.rowsize + x + 1] = 0;
mm.term[y*m.rowsize + x +2] = G;
}
return mm;
}
ImageMatrix convert2blue(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y < m.height; y++)
for (int x = 0; x < m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x + 0];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
mm.term[y*m.rowsize + x +1] = mm.term[y*m.rowsize + x + 2] = 0;
mm.term[y*m.rowsize + x + 0] =G;
}
return mm;
}
ImageMatrix convert2half(ImageMatrix m){
ImageMatrix mm = ImageMatrix(m.height, m.width);
for (int y = 0; y<m.height; y++)
for (int x = 0; x<m.rowsize; x += 3){
BYTE B = m.term[y*m.rowsize + x];
BYTE G = m.term[y*m.rowsize + x + 1];
BYTE R = m.term[y*m.rowsize + x + 2];
BYTE gray = (BYTE)(0.2989*R + 0.5870*G + 0.1140*B + 0.5);
if (y<m.height/2)
mm.term[y*m.rowsize + x] = mm.term[y*m.rowsize + x + 1] =
mm.term[y*m.rowsize + x + 2] = gray;
else{
mm.term[y*m.rowsize + x + 0] = B;
mm.term[y*m.rowsize + x + 1] = G;
mm.term[y*m.rowsize + x + 2] = R;
}
}
return mm;
}
張貼留言