init commit

contests/3contest/1.cpp

@@ -0,0 +1,81 @@
+#include <cmath>
+#include <string>
+
+namespace numbers {
+    class complex {
+        double c_re, c_im;
+    public:
+        complex(double r = 0, double i = 0) : c_re(r), c_im(i) {}
+        explicit complex(const std::string& s) {
+            sscanf(s.c_str(), "(%lf,%lf)", &c_re, &c_im);
+        }
+        double re() const {
+            return c_re;
+        }
+        double im() const {
+            return c_im;
+        }
+        double abs2() const {
+            return c_re * c_re + c_im * c_im;
+        }
+        double abs() const {
+            return sqrt(abs2());
+        }
+        std::string to_string() const {
+            char buf[100];
+            sprintf(buf, "(%.10g,%.10g)", c_re, c_im);
+
+            return std::string(buf);
+        }
+        complex& operator+=(const complex& other) {
+            c_re += other.c_re;
+            c_im += other.c_im;
+
+            return *this;
+        }
+        complex& operator-=(const complex& other) {
+            c_re -= other.c_re;
+            c_im -= other.c_im;
+
+            return *this;
+        }
+        complex& operator*=(const complex& other) {
+            double r = c_re * other.c_re - c_im * other.c_im;
+            double i = c_re * other.c_im + c_im * other.c_re;
+            c_re = r;
+            c_im = i;
+
+            return *this;
+        }
+        complex& operator/=(const complex& other) {
+            double r = (c_re * other.c_re + c_im * other.c_im) / other.abs2();
+            double i = (c_im * other.c_re - c_re * other.c_im) / other.abs2();
+            c_re = r;
+            c_im = i;
+
+            return *this;
+        }
+        friend complex operator+(const complex& C, const complex& other) {
+            return complex(C.c_re + other.c_re, C.c_im + other.c_im);
+        }
+        friend complex operator-(const complex& C, const complex& other) {
+            return complex(C.c_re - other.c_re, C.c_im - other.c_im);
+        }
+        friend complex operator*(const complex& C, const complex& other) {
+            return complex(C.c_re * other.c_re - C.c_im * other.c_im, C.c_re * other.c_im + C.c_im * other.c_re);
+        }
+        friend complex operator/(const complex& C, const complex& other) {
+            double r = (C.c_re * other.c_re + C.c_im * other.c_im) / other.abs2();
+            double i = (C.c_im * other.c_re - C.c_re * other.c_im) / other.abs2();
+
+            return complex(r, i);
+        }
+    };
+
+    complex operator-(const complex& C) {
+        return complex(-C.re(), -C.im());
+    }
+    complex operator~(const complex& C) {
+        return complex(C.re(), -C.im());
+    }
+}

contests/3contest/4.cpp

@@ -0,0 +1,33 @@
+#include <cmc_complex.h>
+#include <cmc_complex_stack.h>
+#include <cmc_complex_eval.h>
+#include <cmath>
+#include <iostream>
+#include <vector>
+#include <string>
+
+using namespace numbers;
+
+
+int main(int argc, char** argv) {
+    complex C(argv[1]);
+    double R = std::stod(argv[2]);
+    int N = std::stoi(argv[3]);
+    std::vector<std::string> record;
+    for (int i = 4; i < argc; i++) {
+        record.push_back(argv[i]);
+    }
+
+    complex I, z, next = C + R;
+    double h = 2 * M_PI;
+    double s = h / N;
+    for (double i = 0; i < h; i += s) {
+        z = next;
+        next = C + R * complex(std::cos(i), std::sin(i));
+        I += eval(record, (z + next) / 2) * (next - z);
+    }
+
+    std::cout << I.to_string() << std::endl;
+
+    return 0;
+}