#!/usr/bin/env python3

import sys


def main():
    mod = 998244353
    sys.setrecursionlimit(10**6)    # 1000 -> 1000000
    def input():  return sys.stdin.readline().rstrip()
    def mi():     return map(int, input().split())    

    n, k = mi()
    intervals = [tuple(mi()) for _ in range(k)]
    N = n
    
    # N: 処理する区間の長さ
    LV = (N-1).bit_length()
    N0 = 2**LV
    data = [0]*(2*N0)
    lazy = [0]*(2*N0)

    def gindex(l, r):
        L = (l + N0) >> 1; R = (r + N0) >> 1
        lc = 0 if l & 1 else (L & -L).bit_length()
        rc = 0 if r & 1 else (R & -R).bit_length()
        for i in range(LV):
            if rc <= i:
                yield R
            if L < R and lc <= i:
                yield L
            L >>= 1; R >>= 1

    # 遅延伝搬処理
    def propagates(*ids):
        for i in reversed(ids):
            v = lazy[i-1]
            if not v:
                continue
            if lazy[2*i-1] + v >= mod:
                lazy[2*i-1] = (lazy[2*i-1] + v) % mod
            else:
                lazy[2*i-1] += v
            # lazy[2*i-1] = (lazy[2*i-1] + v) % mod
            if lazy[2*i] + v >= mod:
                lazy[2*i] = (lazy[2*i] + v) % mod
            else:
                lazy[2*i] += v
            # lazy[2*i] = (lazy[2*i] + v) % mod 
            if data[2*i-1] + v >= mod:
                data[2*i-1] = (data[2*i-1] + v) % mod
            else:
                data[2*i-1] += v           
            # data[2*i-1] = (data[2*i-1] + v) % mod
            if data[2*i] + v >= mod:
                data[2*i] = (data[2*i] + v) % mod
            else:
                data[2*i] += v                  
            # data[2*i] = (data[2*i] + v) % mod
            lazy[i-1] = 0

    # 区間[l, r)にxを加算
    def update(l, r, x):
        *ids, = gindex(l, r)
        propagates(*ids)
        L = N0 + l; R = N0 + r
        while L < R:
            if R & 1:
                R -= 1
                lazy[R-1] = (lazy[R-1] + x) % mod
                data[R-1] = (data[R-1] + x) % mod
            if L & 1:
                lazy[L-1] = (lazy[L-1] + x) % mod
                data[L-1] = (data[L-1] + x) % mod
                L += 1
            L >>= 1; R >>= 1
    
    # [l] を求める
    def get(l):
        propagates(*gindex(l, l+1))
        return data[N0+l-1]



    update(0, 1, 1)
    for i in range(n):
        p = get(i)
        # print(p)
        for l, r in intervals:
            if i + l < n:
                update(i + l, min(i + r + 1, n), p)
    
    # print(query(n-1, n))
    print(get(n-1))


if __name__ == "__main__":
    main()
