a/b % Mod = a % (b * Mod) /b % Mod;

Problem D: qwb与神奇的序列

Time Limit: 1 Sec   Memory Limit: 128 MB
Submit: 892   Solved: 103
[ Submit][ Status][ Web Board]

Description

qwb又遇到了一道题目：

有一个序列，初始时只有两个数x和y，之后每次操作时，在原序列的任意两个相邻数之间插入这两个数的和，得到新序列。举例说明：

初始：1 2
操作1次：1 3 2
操作2次：1 4 3 5 2
……
请问在操作n次之后，得到的序列的所有数之和是多少？

Input

多组测试数据，处理到文件结束（测试例数量<=50000）。

输入为一行三个整数x，y，n，相邻两个数之间用单个空格隔开。（0 <= x <= 1e10, 0 <= y <= 1e10, 1 < n <= 1e10）。

Output

对于每个测试例，输出一个整数，占一行，即最终序列中所有数之和。
如果和超过1e8，则输出低8位。（前导0不输出，直接理解成%1e8）

Sample Input

1 2 2

Sample Output

15

HINT

题目地址：http://115.231.222.240:8081/JudgeOnline/

推出公式：（3^n + 1）/2 * (x + y) % 1e8

思路：在这里2不能求逆元，所以采取这个公式
a/b % Mod = a % (b * Mod) /b % Mod;

这个公式适用于a|b的情况

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<vector>
#include<algorithm>
#include<map>
using namespace std;
const int maxn = 40000 + 10;
#define INF 0x3f3f3f3f
#define clr(x,y) memset(x,y,sizeof x )
typedef long long ll;
#define eps 10e-10
const ll Mod = 100000000;
typedef pair<ll, ll> P;
ll gcd(ll x,ll y)
{
    return y ? gcd(y,x % y) : x;
}
ll pow_mod(ll x,ll n,ll mod_val)
{
    ll ans = 1;
    ll t = x % mod_val;
    while(n)
    {
        if(n & 1)
        {
            ans = ans * t % mod_val;
        }
        n >>= 1;
        t = t * t % mod_val;
    }
    return ans;
}
void ex_gcd(ll a,ll b,ll &d,ll &x,ll &y)
{
    if(b == 0)
    {
        d = a;x = 1;y = 0;
        return ;
    }
    ex_gcd(b,a % b,d,y,x);
    y -= a / b * x;
}
int main()
{
    ll x,y,n;
    while( ~ scanf("%lld%lld%lld",&x,&y,&n))
    {
        ll sum = x + y;
        ll ans = ((pow_mod(3,n,2 * Mod) + 1)%(2 * Mod) /2 % Mod)* sum % Mod;
        printf("%lld\n",ans);
    }
    return 0;
}