結果

問題 No.681 Fractal Gravity Glue
コンテスト
ユーザー 👑 kakel-san
提出日時 2026-10-22 08:59:50
言語 C#
(.NET 10.0.402 + ACL)
コンパイル:
dotnet_c
実行:
/usr/bin/dotnet_wrap
結果
AC  
実行時間 42 ms / 2,000 ms
+ 434µs
コード長 2,354 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 14,052 ms
コンパイル使用メモリ 166,648 KB
実行使用メモリ 32,256 KB
最終ジャッジ日時 2026-10-22 09:00:46
合計ジャッジ時間 9,959 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge2_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 21
権限があれば一括ダウンロードができます
コンパイルメッセージ
  復元対象のプロジェクトを決定しています...
  /home/judge/data/code/main.csproj を復元しました (326 ミリ秒)。
  main -> /home/judge/data/code/bin/Release/net10.0/main.dll
  main -> /home/judge/data/code/bin/Release/net10.0/publish/

ソースコード

diff #
raw source code 

using System;
using static System.Console;
using System.Linq;
using System.Collections.Generic;

class Program
{
    static int NN => int.Parse(ReadLine());
    static int[] NList => ReadLine().Split().Select(int.Parse).ToArray();
    static int[][] NArr(long n) => Enumerable.Repeat(0, (int)n).Select(_ => NList).ToArray();
    public static void Main()
    {
        var n = NN;
        var c = NList;
        var (b, d) = (c[0], c[1]);
        var ans = V(b, d);
        var len = new long[31];
        len[1] = d;
        var size = new long[31];
        size[1] = d;
        for (var i = 2; i < 31; ++i)
        {
            len[i] = len[i - 1] * (d + 1) + d;
            size[i] = (size[i - 1] * (d + 1) + (long)d * i) % mod;
            if (n <= len[i])
            {
                var rest = (long)n;
                for (var j = i; j > 0; --j)
                {
                    if (rest == 0) break;
                    if (len[j] == n)
                    {
                        rest = 0;
                        ans = (ans - size[j] + mod) % mod;
                    }
                    else
                    {
                        var ok = 0;
                        var ng = d + 1;
                        while (ng - ok > 1)
                        {
                            var mid = (ok + ng) / 2;
                            if ((len[j - 1] + 1) * mid <= rest) ok = mid;
                            else ng = mid;
                        }
                        rest -= (len[j - 1] + 1) * ok;
                        ans = (ans - (size[j - 1] + j) * ok % mod + mod) % mod;
                    }
                }
                WriteLine(ans);
                return;
            }
        }
    }
    static int mod = 1_000_000_007;
    static long V(int b, int d)
    {
        return (Exp(d + 1, 2) * (Exp(d + 1, b - 1) - 1 + mod) % mod * Exp(d, mod - 2) % mod + d - b + 1 + mod) % mod;
    }
    static long Exp(long n, long p)
    {
        long _n = n % mod;
        var _p = p;
        var result = 1L;
        if ((_p & 1) == 1) result *= _n;
        while (_p > 0)
        {
            _n = _n * _n % mod;
            _p >>= 1;
            if ((_p & 1) == 1) result = result * _n % mod;
        }
        return result;
    }
}
0