Замена функции Pow

@adel132 · Регистрация: 04.11.2012

Здравствуйте, можно ли как-то иначе заменить функцию Pow?

C

#include <stdio.h>
#include <math.h>
 
//Замена Pow
double Ext (int x, int n) {
    int ext = 1, i;
    for (i = 1; i <= n; i++) {
        ext *= x;
    }
    return ext;
}
 
double Sum(int n) {
    double sum = 0;
    for (int i = 2; i <= n; i++) {
        sum += Ext(-1, i) * (i/(Ext(i, 3) - i));
    }
    return sum;
}
 
int main(int argc, const char * argv[]) {
    int n;
    printf("n = "); scanf("%d", &n);
    printf("sum = %.10f\n", Sum(n));
    return 0;
}

@Pavel Kisliuk · 17.10.2015, 11:29

adel132, возведение в квадрат можно сделать как побитовый сдвиг влево.

@Vtulhu · 17.10.2015, 11:32

C

int power_recursive(int x, unsigned int n) {
    return n == 0 ? 1 : x * power_recursive(x, n - 1);
}
 
int power_effective(int x, unsigned int n) {
    int p = 1;
    while (n) {
        if (n & 1) { p *= x; }
        n >>= 1;
        x *= x;
    }
    return p;
}

@adel132 · 17.10.2015, 12:02 **[ТС]**

Еще фишка в том, что у нас проверяют большие числа. Например число 1 000 000 000 должно как можно быстро вычисляться в этой программе.

@Jewbacabra · 17.10.2015, 12:13

Сообщение от Pavel Kisliuk

возведение в квадрат можно сделать как побитовый сдвиг влево.

это только для числа 2 будет работать

Добавлено через 3 минуты

Сообщение от Vtulhu

power_effective

Я бы еще в начало добавил проверку, при x = -1

C

1
2
3

if (x == -1) {
    return n & 1 ? -1 : 1;
}

@Vtulhu · 17.10.2015, 12:30

Сообщение от adel132

Еще фишка в том, что у нас проверяют большие числа. Например число 1 000 000 000 должно как можно быстро вычисляться в этой программе.

Моя функция power_effective именно такова. Только тип int надо заменить на long или long long, или что вам нужнее.

Добавлено через 1 минуту

Сообщение от Jewbacabra

Я бы еще в начало добавил проверку, при x = -1

Все случаи от -2 до 2 являются специфическими. Можно их предусмотреть, да. Но в общем случае моя функция идеальна. Эффективнее только ассемблер, пожалуй.

@adel132 · 17.10.2015, 22:31 **[ТС]**

нашел еще один вариант, может кому понадобиться

C

long double Step(double x, double y)
{
    return exp(y * log(x));
}

@Pavel Kisliuk 25 / 26 / 16 Регистрация: 27.09.2015 Сообщений: 110
	17.10.2015, 11:29	2
	adel132, возведение в квадрат можно сделать как побитовый сдвиг влево. 0

@adel132 2 / 2 / 2 Регистрация: 04.11.2012 Сообщений: 70
	17.10.2015, 12:02 [ТС]	4
	Еще фишка в том, что у нас проверяют большие числа. Например число 1 000 000 000 должно как можно быстро вычисляться в этой программе. 0

@Vtulhu 429 / 383 / 200 Регистрация: 12.08.2011 Сообщений: 1,610
	17.10.2015, 12:30	6
	Сообщение от adel132 Еще фишка в том, что у нас проверяют большие числа. Например число 1 000 000 000 должно как можно быстро вычисляться в этой программе. Моя функция power_effective именно такова. Только тип int надо заменить на long или long long, или что вам нужнее. Добавлено через 1 минуту Сообщение от Jewbacabra Я бы еще в начало добавил проверку, при x = -1 Все случаи от -2 до 2 являются специфическими. Можно их предусмотреть, да. Но в общем случае моя функция идеальна. Эффективнее только ассемблер, пожалуй. 0

Опции темы

Замена функции Pow

Решение