package main import ( "bufio" "fmt" "os" "sort" "strconv" ) var sc = bufio.NewScanner(os.Stdin) var wr = bufio.NewWriter(os.Stdout) func out(x ...interface{}) { fmt.Fprintln(wr, x...) } func getI() int { sc.Scan() i, e := strconv.Atoi(sc.Text()) if e != nil { panic(e) } return i } func getF() float64 { sc.Scan() i, e := strconv.ParseFloat(sc.Text(), 64) if e != nil { panic(e) } return i } func getInts(N int) []int { ret := make([]int, N) for i := 0; i < N; i++ { ret[i] = getI() } return ret } func getS() string { sc.Scan() return sc.Text() } // min, max, asub, absなど基本関数 func max(a, b int) int { if a > b { return a } return b } func min(a, b int) int { if a < b { return a } return b } func asub(a, b int) int { if a > b { return a - b } return b - a } func abs(a int) int { if a >= 0 { return a } return -a } func lowerBound(a []int, x int) int { idx := sort.Search(len(a), func(i int) bool { return a[i] >= x }) return idx } func upperBound(a []int, x int) int { idx := sort.Search(len(a), func(i int) bool { return a[i] > x }) return idx } func gcd(a, b int) int { if b == 0 { return a } return gcd(b, a%b) } func solve() { N := getI() ans := 0 pre := 1 for k := 32; k > 0; k-- { l := pre r := N + 1 for l+1 != r { m := (l + r) / 2 t := 1 for j := 0; j < k; j++ { t *= m if t > N { break } } if t > N { r = m } else { l = m } } now := l tmp := 0 for i := 1; i <= k; i++ { for j := 1; j <= k; j++ { if gcd(i, j) > 1 { continue } t := N / max(i, j) tmp += t } } ans += tmp * (now - pre) pre = now } out(ans + N*N) } func main() { defer wr.Flush() sc.Split(bufio.ScanWords) sc.Buffer([]byte{}, 1000000) // this template is new version. // use getI(), getS(), getInts(), getF() T := getI() for i := 0; i < T; i++ { solve() } }